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ABSTRACT 


A hybrid control system contains an analog plant and a hybrid {or 
sampled-data) compensator. In this thesis a new conic sector is determined 
which is constructive and can be used to (1) determine closed loop stability, 
(2) analyze robustness with respect to modelling uncertainties, (3) analyze 
steady state response to commands, and (4) select the sample rate. The use 
of conic sectors allows the designer to treat hybrid control systems as 
though they were analog control systems. The center of the conic sector 
can be used as a rigorous linear time invariant approximation of the hybrid 
control system, and the radius places a bound on the errors of this 
approximation. The hybrid feedback system can be multivariable, and the 
sampler is assumed to be synchronous. 

Algorithms to compute the conic sector are presented. Several ex- 
amples demonstrate how the conic sector analysis techniques are applied. 
Extensions to single loop multirate hybrid feedback systems are presented. 
Further extensions are proposed for multiloop multirate hybrid feedback 
system and for single rate systems with asynchronous sampling. 
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1. INTRODUCTION 
1.1 Motivation 

Digital computers are commonly used to control analog systems. 

Examples can be found across the spectrum of engineering disciplines - 
chemical process control, automobile engine control, aerospace systems, 
mass transit systems, electromechanical servomechanisms, and radar tracking 
systems. Much of the technology has been pushed by the aerospace industry 
for use in aircraft, helicopters, missiles, and satellites. Digital 
computers have gained wide acceptance because they are reliable, easily 
reprogrammed, and not subject to drifts in parameter values. Control 
algorithms can be implemented {such as gain scheduling) that are 
difficult or impossible to implement using analog techniques. Of 
primary concern to this thesis is the ability of digital computers to 
mimic (over a certain frequency range) the behavior of analog compensators, 
in spite of their "sampled-data" nature. 

The digital computer is embedded in a compensator that also contains 
a prefilter, sampler, and hold. This collection of analog and digital 
devices will be called a "hybrid compensator." A feedback system that has an 
analog plant controlled by a hybrid compensator is called a "hybrid 
feedback system." A block diagram of a hybrid feedback system ic shown 
in Figure 1.1. The multivariable analog components are modelled as 
Laplace transform matrices, and the multivariable digital computer is 
modelled as a z- transform matrix. 

Depending on where the feedback loop of Figure 1.1 is broken the 
signal is either an analog signal (points 1 and 2) or a discrete 


sequence (points 3 and 4) . At the input and output of the plant the 
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Hybrid Compensator Plant 


Figure 1.1; The hybrid feedback system. 
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physical signals are analog. It is only internal to the hybrid com- 
pensator that the physical signals are discrete sequences. 

If the feedback loop is broken at points 1 or 2 then the transfor- 
mation from these points to the output can be modelled by a continuous time 
linear time varying (CTLTV) operator. This operator is difficult to 
use for analysis techniques because it cannot be represented as a 
transfer function. 

A major simplification occurs if the feedback loop is broken at 
points 3 or 4. Here the hybrid feedback system can be modelled by a 
discrete time linear shift invariant {DTLSI) operator, which in turn 
can be represented by a z- transform. 

The use of z-transforms has led to the development of most of the 
analysis and design techniques in use today for hybrid feedback systems. 

It must be recognized, however, that these techniques can only be 
applied to points 3 and 4 in Figure 1.1. Whenever z-transform techniques 
are used the inherent (often unstated) assumption is that good feedback 
properties at points 3 and 4 (internal to the compensator) imply good 
feedback properties at points 1 and 2 (where the compensator interfaces 
with the real world) . Whether or not this implication is valid depends 
on the choice of prefilter, sample rate, and the nature of the hold 
device. Numerous ad-hoc ways to make these choices have been developed 
concurrently with the z-transform techniques in the sampled-data control 
literature. 

The difficult and important problem of analyzing a hybrid feedback 
system at points 1 and 2 of Figure 1.1 is the primary motivation of this 


thesis: 
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• Develop analysis techniques for hybrid feedback systems that 
can be applied, and used constructively, where the signals 
are analog. 

The analysis techniques are divided into three categories < 

■ closed loop stability 

• robustness with respect to plant uncertainties^" 

• steady state response to commands 

They should apply to multivariable as well as single input single 
output (SISO) hybrid feedback systems. 

The major difficulty with the analysis at points 1 and 2 of 
Figure 1.1 is that the differential equations that describe the hybrid 
feedback system at these points are time varying. One way to lessen 
this difficulty is to approximate the CTLTV system with a continuous 
time linear time invariant (CTLTI) system. Hence the second motivation: 

• Develop rigorous CTLTI approximations of the hybrid feedback 
system . 

The key word here is "rigorous". An approximation is rigorous if it 
applies to all possible inputs and disturbances. It is well known that 
for low frequency inputs a hybrid feedback system can be approximated by 
a CTLTI system. This approximation is not rigorous, however, unless 
it can be shown to be valid (in some sense) for all possible inputs. 

The research motivated by the above requires a branch of system 
theory that is general enough to be used for CTLTV systems and is able 
to make use of CTLTI approximations. Such a branch of system theory 
does exist ([5] to 19]). It was developed to analyze stability and 
robustness properties of nonlinear time varying feedback systems. 

One of the techniques is to represent a nonlinear time varying 
device by a conic sector. A conic sector has a center and a radius. 

feedback system is robust with respect to plant uncertainties if it 
remains stable for all perturbations in a defined set of tne plant. 
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The center is an approximation of the device, and the radius is a 

rigorous bound on the errors of this approximation (i.e. a bound that 

is valid for all possible input-output pairs) , 

These very general conic sector techniques were v nown to exist. 

The conjecture was made (by Gunter Stein) that they are useful for the 

analysis of hybrid feedback systems. Hence the third motivation: 

• Determine whether or not conic sector techniques are useful 
for the analysis of hybrid feedback systems, and if so, 
develop these techniques and demonstrate their usefulness. 

The full generality of these nonlinear conic sector techniques is not 

needed because the hybrid feedback system is linear. The center of the 

conic sector can be used as the CTLTI approximation of the CTLTV parts 


of the hybrid feedback system. 
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1.2 Literature Survey 

The background material for this thesis can be grouped into three 
categories: 1) conventional analysis of digital or sampled-data control 
systems, 2) conic sector analysis of nonlinear time varying feedback systems, 
and 3) conic sector analysis of hybrid feedback systems. There are many 
references in the first category, few in the second, and fewer in the 
third. The third category is the subject of this thesis. 

The conventional analysis of digital control systems uses zvtransform 
theory. The hybrid feedback system of Figure 1.1 is broken at points 
3 and 4, and z-transfoms are used to analyze stability, robustness, and 
performance. No attempt is made here to survey the extensive literature 
on digital control. Much of this literature has descended from journals 
into textbooks and is taught at the undergraduate level. A represen*- 
tative member of this set of textbooks (the one used the most for this 
research) is the text, by Franklin and Powell [4] . The frequency domain 
(z-transform) techniques presented in [4] can be generalized to the 
multivariable case by the use of singular value techniques, just as 
has been done for analog systems [1, 2, 3] . 

The stability of nonlinear time varying systems is a classic problem 
in system theory. Our interest in this problem dates back to the land- 
mark papers of Zames [5, 6] . He showed how conic sectors can be used to 
give sufficient conditions for closed loop stability of nonlinear time 
varying systems.^ 


The landmark papers by Zames are very readable, but for further reference 
see the textbook by Desoer and Vidyasagar [81 . 
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The starting point for Zames* work is the Small Gain Theorem / 
which states that a system is closed loop stable if the loop gain 
(appropriately defined) is less than one. Systems that satisfy this 
condition are of limited practical use because good command following, 
disturbance rejection, and insensitivity to modelling errors requires a 
loop gain much greater than one (most often at low frequencies) . A 
generalization which makes the Small Gain Theorem more useful is made 
possible by the Loop Transformation Theorem , which states that a system 
(call it System 1) can be transformed to another system (call it System 
2) in such a way that :he stability of System 2 implies the stability 
of System 1. Zames [5] was able to show that if System 1 satisfies 
conic sector conditions then it can be transformed to another system 
(System 2) which is stable by the Small Gain Theorem, thereby implying' 
stability of the original system (System 1) . 

Some important generalizations of Zames' work are due to Safonov 
[7] . He generalized conic sectors so that the centers and radii can 
be operators (instead of constant multipliers) . In this thesis conic 
sectors with centers and radii that are CTLTI operators are extensively 
used. 

Safonov [7] goes much further in his generalizations. He defines 
the feedback system of Figure 1.2 that has two "relations" K and G.^ 

The closed loop system is stable if K and -G* (the inverse relation) 
are "topologically separated." Conic sectors are one way to show this 
topological separation. The closed loop system is stable if a conic 

^A "relation" is any subset of a cross product space X x where X and 
y are extended normed linear spaces. 
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U.jf) < K (v) 


Figure 1.2: 


The feedback system defined by Safonov 17] 




sector can be found that K is strictly inside of and -G* is outside of . 

Conic sector stability theory and Liapunov stability theory are 
similar in that they are very general results which cannot easily be 
applied to specific pfoblems 1 . The usefulness of Liapunov stability 
theory depends on the determination of a specific Liapunov function. 
Similarly, the usefulness of conic sector techniques depends on the 
determination of a specific conic sector that is valid for the feedback 
system of interest. No general guidelines exist that help the engineer 
in his search for the "right” Liapunov function or conic sector. 

The conjecture (or hope) of Stein was that a useful conic sector 
could be found for hybrid feedback systems. The first conic sector was 
found by Kostovetsky [17, 18], who was one of Stein's S.M, students at 
MIT. Kostovetsky concentrated on a particular type of hybrid ccmpen- 

2 

sator, which he discovered and named the "optimal hybrid approximation." 
The center of the conic sector is the CTLTI compensator that the hybrid 
compensator is closest to, and the radius is a nondynamic operator 
(i.e. a constant multiplier) . Both the analog and hybrid compensator 
must be open loop stable. Kostovetsky demonstrated by an example that this 
conic sector gives conservative sufficient conditions for closed loop 
stability. Because this conic sector applies to only a very particular 
type of hybrid compensator it is not a "useful" conic sector. 

A second conic sector was found by Stein [18] , which is a 

*The similarity is no accident. Both conic sector and Liapunov stability 
theory are special cases of the general results of Safonov [7] . For a 
description of Liapunov stability theory see Willems [29] , 

2 

The "optimal hybrid approximation" is the hybrid compensator that is 
optimally closest to a CTLTI compensator. "Optimal" is defined as mini- 
mizing the mean square difference of the outputs of the CTLTI and hybrid 
compensators when the input is white noise. The optimal hybrid approx- 
imation turns out not to be practical because the prefilter, computer, 
and hold all contain copies of the CTLTI compensator. 
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generalization of the first conic sector and exists for any stable 
hybrid compensator. Any stable center can be used (poor choices result 
in large radii), and the radius is again a nondynamic operator. 

This second conic sector also proved not to be useful. The problem is 
that the radius is a nondynamic operator, and any such radius results 
in conservative sufficient conditions for closed loop stability. Hybrid 
compensators behave like CTLTI compensators for low frequency inputs but 
not for high frequency inputs. A nondynamic radius cannot distinguish 
between low and high frequency inputs. Because it must be valid for all 
inputs it will too large for low frequency inputs, which in turn results 
in conservative sufficient conditions for closed loop stability. 

Thus began the author’s search for a conic sector with a dynamic 
radius (i.e. a radius that is an CTLTI operator). The search lasted 
for the better part of a year, and ended with the results presented in 
this thesis. Preliminary versions of these result have been presented 


in [18] and [30] . 



-23- 


1.3 Contrlbu ions of Thesis 

Thir thesis brings together two disjoint areas of control theory - 
conventional analysis of digital or sampled-data feedback systems, and 
conic sector analysis of nonlinear time varying feedback, systems. The 
objective is to determine if conic sectors are useful for the analysis 
of hybrid feedback systems. 1 The conclusion is yes, with some restric- 
tions, they are useful. 

The major results are: 

* Proof of the existence of a new conic sector which contains 
a stable hybrid operator. 

* A modification of this new conic sector to create a conic 
sector that the inverse of a hybrid operator (stable or 
unstable) is outside of. 

• Development of analysis techniques based on these new conic 
sectors. Those pertaining to robustness form a unified 
approach to the analysis of robustness properties of hybrid 
feedback systems. 

• Development of algorithms to compute the center and radius of 
the new conic sectors, 

* Demonstration that the new conic sector analysis techniques 
are useful for the analysis of practical hybrid feedback 
systems, including the selection of the sampling rate, 

• Extension of the conic sector results to single loop multirate 
hybrid (SLMRH) feedback systems. 2 

The new conic sectors are more useful than previous versions (17, 18] 
because the radii of the new conic sectors are CTLTI operators and there- 
fore can be represented by Fourier transforms, thus avoiding the inherent 


The hybrid feedback system may be SISO or multivariable and has a single 
hybrid compensator with a synchronous sampler. 


2 

The SLMRH feedback system considered in this thesis (Chapter 6) is SISO 
and has two hybrid compensators in the same loop. The samplers are 
synchronized and have sampled rates that form an integer ratio. 
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conservativeness of conic sectors with radii that are constant 
multipliers . 

The results of this thesis are of interest to both control theoret- 
icians and practitioners. The theoreticians will probably be interested 
in the abstract properties of conic sectors (reviewed here) , the exis- 
tence of specific conic sectors for hybrid operators , the method of proving 
their validity, the gain of the hybrid operator, the signals that achieve 
the gain, the restrictions on the use of conic sectors, and so on. Care 
has been taken rigorously develop the mathematical background. The 
new results are stated as theorems and proved. 

Control practitioners are probably more interested in the new analysis 
techniques - what they do, when they can be used, how they are imple- 
mented, whether or not they are conservative, and whether or not they 
offer anything that can't already be done some other way. This thesis 
addresses each of these issues. The following results are highlighted 
as being of particular interest to practitioners! (1) the ability to 
rigorously approximate a hybrid feedback system by an analog feedback 
system, (2) robustness analysis techniques for single and inultlie.Le 
hybrid feedback systems that make direct use of the analog perturbation 
of the nominal plant (Theorem 3.7), and (3) the use of the conic sector 
radius to help select the sample rate. Both theoreticians and practitioners 
should be interested in the numerical examples, which complement and help 
to explain the theoretical results. 
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1.4 Summary of Thesis 

The material in this thesis is organized as follows i 

Chapter 2s Preliminary mathematics 

Chapter 3: New conic sectors for hybrid compensators 

Chapter 4; Algorithms to compute the radius of the new conic 
sectors 

Chapter 5t Examples of the use of conic sector analysis techniques 

Chapter 6: Multirate sampling issues 

Chapter 7 : Summary and extensions 

Each of these chapters is now summarised. 

Chapter 2 The mathematical review begins in Section 2,2 with 
analysis techniques for multivariable CTLTI feedback systems. The 
use of the loop transfer function and singular values is stressed. Plant 
uncertainties are modelled as additive or multiplicative perturbations 
that are bounded in magnitude. The relationship between multiplicative 
perturbations, phase margins, and gain margins is discussed. The goal 
(not fully realized) of conic sectors analysis techniques is to emulate 
these powerful multivariable techniques recently developed for CTLTI 
feedback systems. 

The emphasis of Section 2.3 is on the conventional analysis of 
hybrid feedback systems. The structure of the hybrid compensator is 
described, and it is shown that the CTLTV input-output transformation 
can be modelled by time and frequency domain methods. The z-transform 
analysis techniques are presented for multivariable systems using singu- 
lar values. Discretizing analog perturbations of the nominal plant is 
discussed. 


The difficult transition is then made to conic sector analysis 
techniques for general (i.e. nonlinear time varying) feedback systems. 
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Th« level of presentation in Section 2.4 Is more ma thematically precise. 
Definitions are given for relations, operators, gain, stability, and 
conic sectors, The sufficient conditions for closed loop stability are 
fully developed. Preliminary work by Stein {18] on using conic sectors 
to analyze command response is presented. 

As a breather before moving on to the new results, in Section 2.5 
it is shown how conic sectors are applied to multivariable CTLTI feedback 
systems. Results due to Safonov [9] are presented that prove the exis- 
tence of conic sectors fox analog plants which have bouided additive or 
multiplicative perturbations. 

Chapter 3 The major theoretical result of this thesis is the 
proof of existence of a new conic sector which contains a stable hybrid 
operator.* The most difficult step of this proof is a frequency domain 

2 

inequality (Lemma 3. A) which makes use of Lebesque Dominated Convergence 
and the Cauchy-Schwartz inequality. Theorem 3.2 is modified to show 
(1) ;he existence of a conic sector that contains the loop transfer 
operator and (2) the existence of a conic sector that the inverse hybrid 
operator is outside of (Theorem 3.6). 

Chapter 3 contains 8 theorems, which form the basis for the 
following conic sector analysis techniques! 

• An upperbound on the gain of a hybrid operator. For SISO 
hybrid operators the upperbound equals the gain (Theorem 3.1) 

• Sufficient conditions for closed loop stability (Theorem 3.3) 

*The hybrid compe nsato r is modelled mathematically as a hyo. id operator . 

Shankar Sastry and Marcel F, Coderch pointed out tne need to check the 
conditions of Lebesque Dominated Convergence and helped to do so. 
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• Sufficient conditions for robust closed loop stability, i.e. 
closed loop stability for all perturbations in a defined 
set of the nominal plant (Theorems 3.4, 3.5, and 3.7) 

• An upperbound for the steady state error of a hybrid feedback 
system (Theorem 3.8) 

A distinction is made as to whether the hybrid operator is inside 
of or outside of a conic sector. The inside conic sector techniques 
(Theorems 3.1 to 3.5) have the major restriction that the hybrid 
operator must be open loop stable. The outside conic sector techniques 
(Theorem 3.6 and 3.7) are less restrictive because they only require 
that the hybrid feedback system is closed loop stable (which can be 
determined by z-transform techniques) . 

When the hybrid operator is placed inside of a conic sector then 
Theorems 3.4 and 3.5 form a unified approach to the analysis of robustness. 
Both the hybrid compensator and the actual plant are approximated by 
CTLTI operators (centers of conic sectors) , and the respective errors 
of these approximations are modelled by bounded perturbations (radii 
of conic sectors) . The point of view taken is that a hybrid compensator 
is supposed to mimic a CTLTI compensator. The ex rent to which it does 
not is a sour ; of error. By comparing the sizes of the radii as a 
function of frequency the errors due to the use of a hybrid compensator 
can be compared to the errors due to uncertainties of the plant (unmodelled 
higher order dynamics, time delays, and so on). 

Chapter 4 Conic sectors would not be useful if it was not possible 
to compute the radius. 1 The most difficult parts of computing the radius 

Computing the center is no problem because any centei' can be used, 
although some centers are better them others. 
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are the sunning of infinite and double infinite series of Fourier 
transforms (each term shifted in frequency) . Three ways to compute 
the radius are 

* Sum a finite number of terms of single infinite series (the 
double infinite series can be broken down into several single 
inf ini ts series) . 

* Sum a finite number of terms of the double infinite series. 

* Find an exact analytical solution to the single infinite series. 
These approaches are discussed and compared. Bounds are placed on the 
remainders of the truncated infinite series. Several cases where exact 
analytical solutions can be found are presented. In the examples the 
second approach is used. 

Chapter 5 The conic sector analysis techniques are demonstrated 
in the examples of Chapter 5. They are shown to work well for SISO 
hybrid feedback systems but to be conservative for multivariable hybrid 
feedback systems. More work is needed to remove this conservativeness. 

Section 5.2 contains an extensive example of a SISO hybrid feedback 

nd 

system. A lead-lag compensator is designed for a stable 2 order plant 
using classical control techniques , and then the analog compensator is 
discretized to form the digital computer part of the hybrid compensator . 
Each of the 8 theorems of Chapter 3 are used to analyze the hybrid feedback 
system. In addition, the size of the radius (as a function of frequency) 
is used as a measure to 

* select sample rates 

* compare discretization techniques 

* compare different types of errors (due to the use of a hybrid 
compensator and due to time delays in the plant) 
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A control problem with an integral compensator is presented in 
Section 5.3. The hybrid operator cannot be placed inside of a conic sector, 
but the inverse hybrid operator can be placed outside of a conic sector, 
which is then used to analyze robustness. 

Section 5.4 contains a SISO and multivariable version of a hybrid 
feedback system for controlling motion in the pitch axis of a high 
performance aircraft. The analog compensators are designed using linear 
quadratic Guassian (LQG) techniques, and then the analog compensators 
are discretized to form hybrid compensators. The conic sector analysis 
techniques must be slightly modified when the loop transfer operator 
is placed inside of a conic sector to account for the fact that the 
plant is open loop unstable. The multivariable hybrid operator is 
placed outside of a cone, which is shown to result in a conservative 
robustness margin. 

Chapter 6 extends the conic sector results of Chapter 3 to the 
single loop multirate hybrid (SIMRH) feedback system of Figure 6.1. 

These extensions are useful for a limited class of feedback systems, 
and should be considered preliminary results for the more general and 
more important problem of analyzing multiloop multirate hybrid feedback 
systems . 

The new conic sectors are based on a frequency domain description 
of the input-output transformation of the cascaded hybrid compensators . 

The key observation is that the frequency domain description (6.3) has 
a similar structure to the description (2.51) of a single rate hybrid 
compensator operating at the slower of the sample rates. The detail? 
are presented in Section 6.2. An example is presented in Section 6.3. 

Chapter 7 contains a summary and suggestions for future research. 
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In the summary the entire thesis is described as revolving around the 
major result of Theorem 3.2. The suggested topics for future research 
are: 

* Removal of open loop stability restriction 

* Generalization of conic sector techniques to sector techniques 

* Less conservative multivariable robustness margins 

* Synthesis techniques 

* Multirate sampling issues 

* Asynchronous sampling is lues 


Finite word length issues 
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2. MATHEMATICAL PRELIMINARIES 
2.1 Introduction 

The background material for this thesis is reviewed in this chapter. 

•Hie starting point is in Section 2.2, where a set of results are presented 
about multivariable analog linear time invariant (LTI) feedback systems. 

The emphasis is on analysis techniques to determine closed loop stability, 
cocmand response, and robustness with respect to modelling uncertainties . 

The results presented in Section 2.3 are from the conventional 
analysis of multivariable hybrid feedback systems. The components of the 
hybrid compensator are described, and then it is shown how to use 
z-tranforms to determine closed loop stability, command response, and 
robustness. 

A much more general feedback theory is presented in Section 2.4. 

The components of the feedback system are modelled as "operators" and 
"relations." Closed loop stability, command response, and robustness 
are analyzed with "conic sectors." These analysis techniques are 
general enough to apply to both analog and hybrid feedback systems. 

Section 2.5 reviews how conic sectors are used to analyze analog feed- 
back systems. The contribution of this thesis is to use conic sectors to 
analyze hybrid feedback systems, which is done in Chapter 3. 

The general feedback theory is presented in a rigorous fashion. 

Starting from a set of definitions, the major results are presented in 
a series of lemmas. When it is felt to be necessary the proofs are 
included in the appendix to Chapter 2. An attempt is made to keep the 
presentation at the minimum required level of generality. This is 
consistent with the main thrust of this thesis, which is to apply a 
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general ft dback theory to a specific problem. 

The reader is assumed to be familiar with Laplace and z-tranaform 
theory, state apace techniques, and single-input single-output (STSO) 
analog and digital control* The use of the loop transfer function and 
the return difference equation for the analysis of feedback systems is, 
stressed. Students of control theory will probably hare had exposure 
to these subjects by the end of their undergraduate curriculum. 

The reader is assumed to be less familiar with the use of singular 
values for the analysis of multivariable control systems, and even less 
familiar with general feedback theory. Some exposure to the theory of 
Hilbert Spaces will be helpful for an understanding of the general 
feedback theory. 

There are numerous well-written journal articles and textbooks that 
can be used for further reference. The use of singular values for the 
analysis and design of multivariable analog feedback systems is explained 
in [1, 2, 3] . The analysis and design of SISO digital feedback systems 
is explained in [4] . Conic sector analysis of general feedback systems 
is explained in [51 to [9] . 

2.1.1 Notation 

if 1 , C n , 11™°°, C 0X311 » finite dimensional real and complex Euclidean 

spaces 

■ n-dimensional space of square integrable functions 

l 2 ° “ extended L 2 ° space 

* space of infinite dimensional vectors 



real numbers > o 
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A e C 0 ”” = matrix (underlined capital Foman letters) 

a, a € C° * vector (underlined small Roman and Greek letters) 

a, o e C° = scalar (small Roman and Greek letters) 

A 6 L 2 n e x L™ e = relation or operator (capital script letters) 

a 6 ij 1 - function (small Roman letters underlined by a tilda) 

— 2e 

A * - matrix inverse 
«»» 

T T 

A , a = matrix and vector transpose 
H H 

A , a = matrix and vector Hermitian (complex conjugate transpose) 
A 1 = inverse of relation or operator 


|| a || = Euclidean vector norm 

|| A || = matrix norm, induced by Euclidean vector norms 

|| a || . = L function norm 

2 2 

II a || - truncated function norm 

II A || . = operator norm, induced by functions norms 

2 

0 . [A] = minimum singular value of A 

nun - 

o [A] = maximum singular value of A 

max - - 


Associated with each LTI operator A (or with each a 6 L^) 
A(s) = Laplace transform 

A(jw) = wplace transform evaluated at s = jw 
A (t) = inverse Laplace transform, impulse response 
A (nT) = sampled version of A(t) 


A(z) = z-transform of A(nT) 
A* (s) = ~ I A (s-j ~ k) = 


A(z) 


z=e 


ST 
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The following Rotational abbreviations are usedt 

» F(jw - jto^k) (subscript k sometimes replaced by n) 

E<0 - ? (*> 

k k=-® 

I (•, = S (•) 

n^k n»-«° 

n/k 
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2.2 Analysis of Analog Feedback Systems 
2.2.1 The Analog Feedback System 

The analog feedback system is shown in Figure 2.1. The plant and 
compensator are multivariable LTI systems with Laplace transform matrices 
G(s) and K(s) . The output ^(s) and the error e(s) are related to the 
command input r(s) by 

y (s) ■ GK (I + GK)' 1 r (s) (2.1) 

e(s) = (I + GK) -1 r (s) (2.2) 

Analysis techniques repeatedly use the loop transfer function and 
the return difference equation, this due to the fact that closed loop 
properties can be determined by how signals propagate around the loop. 

For multivariable systems the loop transfer function differs depending 
on where the loop is broken. Natural places are where signals enter 
the loop, such as points (1) and (2) of Figure 2.1. Most of the atten- 
tion in this thesis is focussed on point (1) , where the loop transfer 
function is 

T x (s) = G (s) K(s) (2.3) 

and the return difference equation is I + T^(s). 

2.2.2 Plant Uncertainty 

The Laplace transform matrix G(s) is a mathematical model of a 
real system. The model will always be inexact, for many reasons 
including: 

• inaccurately measured or slowly time-varying parameters 

• unknown or purposely neglected high frequency dynamics (such as 
bending modes of mechanical systems) 
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• neglected time delays (far away right half plane zeros) 

* changes in the operating point about which a nonlinear model 
is linearized 

The actual plant is assumed to be 5(s), which is related to G(s) by 
the additive perturbation 

Sts) * G(s) + E (s) (2.4) 

It is assumed that G(s) and G(s) have the same number of unstable poles, 

which will definitely occur if E (s) is stable and may occur if E (s) 

"d “ft 

is unstable. 

The only other information known about E (s) is the following 

"3 

constraint on its magnitude: 

0 (E (jw)) < l (oj) for all to (2.5) 

max —a a 

The symbol "ct [E ( jw) ] " indicates the maximum singular value of 

IHctX “2 

E ( jo>) , as defined and discussed in the next subsection. Uncertainties 

*“ft 

that are bounded in magnitude, such as E (joi) , are called "unstructured" 

"ft 

uncertainties . 

It is also possible to characterize plant uncertainty by the 
multiplicative perturbation: 

G (s) = G(s) IT + E (s)] (2.6) 

-m 

It is assumed that G(s) and G(s) have the same number of open loop 

unstable poles. The only other information known about E (s) is 

-m 

0 IE < jco) 1 < l ((d) for all (D (2.7) 

max -m m 

The two types of perturbations are related by 

E (s) = G(s) E (s) 

-a _ -m 


( 2 . 8 ) 
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and their bounds are related by 

t«o) < J [G ( jco) ] l (co) (2.9) 

a max - m 

Whether the additive or multiplicative perturbation is used is often 
a matter f convenience, and results using one type of perturbation 
can usually be expressed using the other. The additive perturbation 
is more suited to conic sector analysis, whereas the multiplicative 
perturbation is more natural to use for an analysis of phase and gain 
margins. 

Analysis techniques usually asstime that the compensator Laplace 
transform matrix K(s) is exactly known- the reason being that it would 
be foolish to build an uncertain compensator. Nevertheless, if any 
perturbations exist they can be combined with the plant perturbations 
to create a perturbation of the loop transfer function. For example, 
if the actual plant and compensator are 

5(s) = G(s) + E (s) (2.10) 

" & 

K(s) = K\s) + E^fs) (2.11) 

then the actual loop transfer function is 

f^s) = 5 ( s) fc(s) ■ G(s) K(s) + 5 ta (s) (2.12) 

where 


E. (s) * G (s) E. (s) + E (s) K(s) + E (s) E. (s) (2.13) 

-ta - -ka -a - -a -ka 


a (E ( joi) ] < c [G(jo))l L (<o) + l (w)*o [K(jw) 
max -ta max - J ka a max - 


+ l (to) * l, (co) 
a ka 


(2.14) 
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More sophisticated (less conservative) techniques ft. r combining pertur- 
bations are given by Safonov [10] and Doyle [11] . 

2,2.3 Singular Values 

The singular values of a matrix A e C nxni are defined as: 

(? i (A) A |A (A H A)] 1/2 for i - 1, . , . ,m (2.15) 

The notation X^(*) indicates eigenvalue. Cf the many properties of 

singular values, the ones that used in this thesis are listed in Table 

2.1. For further reference see [12], 

Singular values are used to determine the gain of a matrix. Suppose 

the matrix A is multiplied by the vector x. The gain of A in the 

direction x is the ratio of the Euclidean vector norms^ |j Ax || _ / |j x |( . 

The maximum gain is 0^ (A) , which by property 2 of Table 1 is equal to 

the induced vector norm || A ||. The minimum gain is • If a square 

matrix A has a . (A) = o then the matrix is singular and cannot be 
min - 

inverted . 

Singular vclues can be used to give a quantitative measure of the 
2 

"size" of a matrix . The matrix A is "large" if <7 . (A) >> 1, and 

mm 


The Euclidean vector norm is defined 







1/2 


where x^ for i - l,...,n are the elements of the vector x. 


The number "1" is used for comparison because the return difference 

equation has the form I + A. Important properties (used for performance 

analysis) are that (I + A) s a . (A) if c , (A) » 1 and that 

1 mm - - nun „ 


a (i 

max - 


+ A) ; 1 if 


°max <*> 


« 1 . 


if P . (A) 
min - 



( 1 ) 

( 2 ) 

(3} 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

( 9 ) 

( 10 ', 
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CL is real and nonnegative, for i - l,...,n (the o^’s are ordered 

from maximum to minimum, with o, =0 and a • a . ) 

1 max n min 


0 (A) = A 

max - - 


11*5 ll E 

max 

x|i E <°° 11*11 E 


roax II Ax || 

x|l E -l 


lift* ll E 

0 . (A) “ min 

min ' i|x|L<» llx" 


max 


,„ E Hx|I e -1 


|ax|I e 


°min ( ^ Me - He ^ °max (A) Me for all il x |1 £ < - 


0 (A) = — — i > if A exists 


max - a . (a-1) 
min'- / 


-1 -1 

(A + E) exists if A exists and o (E) < 0 . (A) 
- - max - min - 


(Triangle inequality) t fT (A + E) < a (A) + o (E) 

max - - - max - max - 


(Fan’s Theorem [40]): o. (A + E) < a, (A) + o (E) for i = 1 ,n 

i - i - max - ' ‘ ' 

(Properties 9 and 10 are consequences of Paul’s Theorem) 
a. (A) - 1 < a. (A + I) < a ± (A) + 1 for i=l,...,n 

O . (A + E) > 0 . (A) - O (E) 

min - - mm - max - 


I 

i 

j 

* 



* 

1 


Table 2.1; Properties of singular values 
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"small" if a (A) << 1. 
max - 

Singular values can be computed by stable and reliable algorithms 
that are easily available well tested, and well documented [13] . 

Many of the analysis techniques for multivariable systems use 
singular values. Readers unfamiliar with singular values should note 
that for complex scalars*. 

Il»!l- W a) ’ < W a) - W (2 - 16) 

When working with SISO systems most of the singular value results can 
be replace 1 with absolute values. 

2.2.4 Command Response 

Specifications for control systems often include a constraint on 
how well the output y (t) follows certain types of inputs r(t) . Such 
specifications can be stated as constraints on the singular values of the 
loop transfer function. Specifications for disturbance rejection and 
sensor noise attenuation can also be stated as constraints on singular 
values. In this thesis, however, only the coimnand response will be 
analyzed. ^ 

Consider the set of commands that have energy only in the fre- 
quency range oo < oi < co^. The response to these commands will be 


Specifications for command response (which usually take the form of 
high gain at low frequency) are the same as specifications for the 
attenuation of disturbances that are added to the output. Therefore, 
results for command response can be restated as results for disturbance 
rejection bv replacing r(t) by r(t) - d(t). Specifications for sensor 
noise attenuation must be treated differently (low gain at high 
frequency) . For a good treatment of the various performance specifi- 
cations see [1] . 
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good if over this frequency range : 


^(jw) as r(jcu) 

<= e(j<o) = o 

<= [I + GK(jO))] _1 ~o [by (2.2)] 

<= a [I + GK(jOJ)]’' 1 « 1 

max - -- 

<= O . [I + GK(jtu)] >> 1 (property 5, Table 2.1) 

mm - ’ — 

<-= a m i n (GK(jW) ] >> 1 (property 9 , Table 2.1) (2.17) 

Hence, a specification for command response can be stated 


o . [GK(ja>)] > p (co) for co < w < w. (2.18) 

nun — - o 1 

For example, if the command is the step function r(t) = a for t > o 
then the steady state error is 

e (t) = II + GK(jo)]" 1 a = b (2.19) 

The relative error, expressed as a ratio of vector norms, has the 
upper bound 


*W E -1 

<0 [(I + GK) ] 

H max 

E 


1 = 1 

a . [I + gk] ~ p (o) 
nan - 


( 2 . 20 ) 


where it is assumed that p(o) » 1. 

The command response specification must be satisfied not just for 
the nominal plant, but also for every possible perturbation of the 
nominal plant. This is guaranteed by increasing the lower bound of 
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equation (2.18) : 


min — - l - L 


P to) 


to) /a EG ( ja>) ] 

a max - 


( 2 . 21 ) 


Therefore, if the additive perturbation to) is 10% the size of 

a [G(jai)] then command response specification must be increased by 
max ~ 

about 10%. 

2.2.5 Stability 

It is important to kr >./ whether or not the closed loop system is 
stable. This is relatively easy to check for analog feedback systems. 

In fact some desig:- methods, such as the linear quadratic Gaussian 
method [Athans, 14] , guarantee closed loop stability (under mild 
assumption;-) . 

The closed loop system is stable if all of the closed loop poles 

are in the left half plane (i.e. Rets^ < o for all closed loop poles 

s^.) . It is easiest to check this condition if a state space description 

is known for G(s) K(s) . In this case the closed loop poles are the 

eigenvalues of the closed loop system matrix. 

It is not always feasible to compute the closed loop poles. This 

is the case when G(s) K(s) contains polynominals of very high order or 

when they contain infinite dimensional terms such as the pure time 
-sT 

delay e . For these cases a better way to check closed loop stability 
is to use the multivariable Nyquist theorem. This is due to Rosenbrock 
[15] . A nice statement of the theorem is given by Lehtomaki [3] . 
Basically, this is applied by plotting the imaginary versus the real 


It is assumed that i { ) < o tG(jco)] and that o . [GK(jco)] » 1. 
See [1]. a max " min ” 
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part of -1 + detll + GK(jtu)] (as 0) varies) and counting the msnber of 
encirclements of the -1 point. 

2.2.6 Robustness 

It is not enough just to know whether or not the closed loop system 
is stable. If it is stable it must also be known how close it is to 
becoming unstable. This is the subject of stability-robustness (often 
referred to just as robustness) . The closed loop system must be stable 
for all possible perturbations of the nominal plant G(s), in other 
words for all possible 5(s) of equation (2.6). Sufficient conditions 
for this to be true have been derived by Doyle, Stein, Lehtomaki, and 
others 11, 2, 3]. 

The sufficient conditions for robust stability have two parts. 

The first is that the nominal plant G(s) must result in closed loop 
stability. This condition is called "nominal closed loop stability." 

The second condition can be expressed as a constraint on the maximum 
singular value of the multiplicative perturbation: 

0 [E GK (I + GK) _1 (jOJ)] < 1 for all 0) 
max -m — - J 

<— 0 [E (j< 0 )] < 0 . (I + (GK) ^ ( jOJ) ] for all 0) (2.22) 

max -m min - 

Most people find the nominal closed loop stability condition to be entirely 
reasonable and the singular value inequality to be somewhat less that in- 
tuitive. The basic idea of the second condition is that it guarantees that 
the return difference equation I + GK(jw) remains invertible for all poss- 
ible G's in the set defined by E , 

— m 

There are a plethora of other ways to express the robustness 
results. For example, if the additive perturbation is used, the 
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singular value inequality changes toi 

-1 

o [E K (I + GK) (jco)l < 1 for all co 
max -a- 

<= 0 IE <ja»] < O . [K 1 ( jaj) + G(jco)] for all U (2.23) 

max -a min - - 

It is not always obvious which of the various results is best to use. 

For a method to derive these results and for some guidance on which to 
use see Lehtomaki [3] . 

2,2.7 Phase and Gain Margins 

These terms have been handed down to us from classical control 
theory and are used to characterize the uncertainty of SISO plants. 

Let the nominal and actual SISO plants be related by 

g (s) = y(s) e(s) (2.24) 

— 

To verify a phase margin specification let e(s) = e . If the closed 
loop system is stable for all g(s) such that |4>j < 45° then the closed 
loop system has a phase margin of 45°. Similarly, to verify a gain 
margin specification let e(s) = Z. If the closed loop system is stable 
for all g(s) such that Z^ < Z < Z^ then it has a gain margin of Z^i • 

The gain margin is usually expressed in terms of decibels: [20 log^g(£^) • 
20 log 10 (£ 2 )]. 

Phase and gain margins can be expressed as multiplicative per- 
turbations. For SISO plants the multiplicative perturbation is 

g (s) - g (s) II + e (s) ] (2.25) 

m 

which is related to e(s) by 

e (s) = e (s) - 1 
m 


(2.26) 
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A phase margin specification of 45° corresponds to the following bound 
on the multiplicative perturbation: 


| e (ja» | < |e j45 * -l| - l (ai) - .77 
m m 


(2.27) 


Thus, if the closed loop system is nominally stable and if 


1 + (gk) 1 ( jco) ( > .77 for all w 


(2.28) 


then the closed loop system has at least a 45® phase margin. 
The correspondence can be turned around. Suppose that 

|l + (gk)" (jo) | > l (w) = a for all o 

HI 

Then the closed loop system has the following guarantees 1 : 

2 

guaranteed phase istrgin = ± arccos (1 — - ) 


guaranteed gain margin = [- 


-] 


(2.29) 


(2.30) 


l 1-kx ' 1-a 

These are guarantees - the actual margins may be better. 

The phase and gain margins cannot be simultaneously achieved. They 
indicate robustness with respect to pure phase or pure gain perturbations. 
It is possible, however, to analyze robustness with respect to simultaneous 
phase and gain changes [16). Let e(s) = £e^, a.id then find regions in 
the 20 log^g (t) x 4> space where |£e^ -l| < a. These regions are plotted 
in Figure 2.2. Suppose that a ■ .8. Then any combination of t and <J> 
inside of the ellipse marked a = ,8 will not affect closed loop stability. 


1 i jd> | , 

This is shown by finding (ft such that |e -1| < a and by finding L 
such that |£-1 I < a. 



phase angle 
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gain, dB 

Simultaneous phase and gain changes that correspond 
to a constant multiplicative perturbation of ^ (W) 


Figure 2.2: 
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The guaranteed phase and gain margins lie on the boundary of the ellipse. 

Phase and gain margins can be generalized to multivariable systems 
[7, 3]. This is done by inserting the diagonal perturbation 

£*-?) = diag (e^s)} (2.31) 

In which case, if 

°min + < a (2.32) 

then the e^'s can simultaneously undergo phase and gain changes as 
indicated in Figure 2.2. 

The generalizations of phase and gain margins are sometimes subject 
to misinterpretation. The danger lies in restricting attention to diagonal 
perturbations. The closed loop system may be very robust with respect to 
diagonal perturbations, but sensitive to off diagonal perturbations. 

When this happens then a of (2.32) will be small, leading to conservative 
guarantees for perturbations that are restricted to be diagonal. 
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2.3 Conventional Analysis of Hybrid Feedback System* 

2.3.1 The Hybrid Compensator 

A hybrid compensator consists of a pre filter, .^mpler, digital 
computer, and hold. A block diagram is shown in Figure 2.3. The word 
"hybrid" emphasizes that the compensator has both analog and digital 
parts. Both the input and output are analog signals, so from an input- 
output point of view the hybrid compensator is an analog device. In- 
ternally the signals are represented by discrete sequences, so from 
this point of view the hybrid compensator is a digital device. 

It is the sampler that converts an analog signal into a digital 
sequence, and it is the sampler that complicates the analysis of hybrid 
feedback systems. Associated with the sampler are the signals e^(t), 
e^lnT) , e d *(t), and their respective transforms e^s), 63 (z), ^ 
e^*(s). We show below how they are related. 

The input to the sampler is the analog signal e^tt). The sampler 
is periodic and outputs a sample every T seconds, so the output is the 
discrete sequence e^(nT). Another way to represent the output is the 
input multiplied by a train of impulses: 

e* (t) = £ e,(t) 6 (t-nT) (2.33) 

— d — d 

n 

where 6 (t) is an impulse at time t=0. The following identities are 
well known, and are derived in [4] : 

^•(s) = - £ e^(s-ju) s n) , where (jJ g = — (2.34) 

n 


®d* (s) 


ejU) 


z=e 


sT 


(2.35) 
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Prefilter Sampler Digital Hold 

Computer 


'■■■» analog signal 
digital seguence 


Figure 2.3: The hybrid compensator 
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The star operation may be considered a mathematical operation 
independent of its above use with samplers. Its definition for Laplace 
transform matrices is 

AMs) A ^ l A(s-jO) n) (2.36) 

L n 8 

Two of the properties of the star operation are 


AMs-ju> n) = AMs) for any integer n (2.37) 

s 

[AMs) B (s) ] * = AMs) BMs) (2.38) 


Every z-transform matrix A(z) has associated with it on AMs) defined by 


AMs) A A(z) 


sT 

z=e 


(2.39) 


An AMs) matrix defined by (2.39) obeys properties (2.37) and (2.38), 
and can be used interchangeably with an AMs) defined by (2.36). In 
this thesis the star notation is preferred to the z-transform rotation 
(e.g. Figure 2.3). 

The first part of the hybrid con^ensator is the prefilter. It is 
a linear time invariant (LTI) system which has the Laplace transform 
matrix F(s) . Its main purpose is to low pass filter the input and 
reduce aliasing. If F e(juj) is nonzero for ju>| > Ti/T then aliasing 
will occur, as readily seen from (2.34).* 


Aliasing occurs if J(Fe)Mjoo)| > ] Fes ( jco) 1 for 0 < Cu < tt/T. Aliasing is 
undesirable phenomena for control systems. One of the effects of 
a-Liasinq is that high frequency inputs are interpreted as low frequency 
inputs, which is particularly true of high frequency noise. 
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The digital computer is a linear shift invariant 1 device with the 
z- transform d(z) . How the z-transform is derived is a synthesis prob- 
lem and is not the main concern of this thesis. The two basic approaches 
are (1) to discretize an analog compensator and (2} to discretize the 
plant and synthesize D(z) by a direct method. 

In the examples of chapter 5 the Tustin with prewarplng method 
[4, p. 55] is used to discretize K(s) . The z-transform for the computer 
is set equal to 


D(z) 



a 


z-1 

z+1 


where 



(2.40) 


Tustin with prewarping has the property that (respectively) : 


D(z) 


z=l, 


jOJjT 

e 


K(s) 


s*0. 


(2.41) 


The analog compensator is perfectly matched at s=0 and s=jOJ^. 

The prewarped frequency should in some respect be a "special" 
frequency, such as the natural frequency of a pole or zero, the 
frequency of maximum phase lag or lead, or the closed loop bandwidth. 


^Linear shift invariance is the discrete equivalent of linear time in- 
variance. It the input sequence is shifted an integer number of sample 
periods then the output is the same except for being shifted the same 
amber of sample periods. Only linear shift invariant operators have 
z-transforms defined for them. 



» ■ 4 * 


-bJ- 


^SSs 


Tustin with prewarping has the advantage of being easy to compute. 
Suppose that K(s) has the state space realization 


X_ 

1 


* Ax + B u 

*= r x + Du 




( 2 . 42 ) 


Then the state space realization of D(z) is 


x(n+l) 

= F x(n) 

+ G u. 
d 

*d (n) 

= H x(n) 

+ K u, 
d 


(n) 

(n) 


where 



( 2 . 43 ) 


( 2 . 44 ) 


The hold device transforms a digital sequence into an analog signal. 
The hold is modelled as a LTI sv^cem with the Laplace transform H(s) . 

The input to the hold is a train of impulses u^Mt) with Laplac - trans- 
form u *ls), and the output of the hold has the Laplace transform 
— d 
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The most consnon 
response and Laplace 


type of hold is the zero-order- hold, 
transform (for the S1S0 case) are 


Its impulse 


h ,t) 


1 0 < t < T 

0 elsewhere 


h(s) = j (1 - e' sT ) 


A first order analog approximation of h(s) is 


(2.45) 


(2.46) 


h (s) 

a 


2 

s + 2/T 


(2.47) 


Magnitude and phase Bode plots of h(ju>) and h (ju>) are shown in Figure 
2.4 (for T = .6283- w = 10). Below j o> rad/sec the approximation 

S 4 S 

is very good. 

2.3.2 The Hybrid Operator 

An operator transforms input signals (belonging to a set of al- 
lowable input signals) to unique output signals. 1 An operator is a 
mathematical moael, as opposed to a physical system. The hybrid 
operator is a mathematical model of a hybrid compensator. It is given 
the symbol K, and the transformation from an input signal e to an out- 
put signal u is represented by 


u = K e 


A more precise definition is given in Section 2.4. 


(2.48) 
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The hybrid operator is a linear time varying (LTV) operator . The time 
variations are due to the sampler. The same input signal shifted in 
time by a fraction of sample period results in a different sequence 
of samples. 

A time domain description of the hybrid operator is given by the 
convolution 

(0)d0 (2.49) 

where K(t,0) is an impulse response matrix. It is periodic in the 
sense that 

K(t+Tn, 9+nT) = K(t,0) for any integer n (2.50) 

The same input shifted by an integer number of sample periods results 
in the same output shifted by the same integer number of sample periods. 

It is this time domain description that was used by Kostovetsky 
{17] (see also [18]) to investigate some properties of hybrid operators. 
He was able to show that the gain of the hybrid operator is unbounded 
as the prefilter approaches an impulse (F(s) I). Also, he showed 

how to select a prefilter, digital computer, and hold such that K(t,0) 
is optimally close to a specified LTI impulse matrix K(t) . 

This thesis uses the following frequency domain description of the 
hybrid operator. If the input signal has the Laplace transform e(s) 
then the output signal has the Laplace transform 

u(s) = 


u(t) 


-r* 


K(t ,0)e 


H (s) D* (s) (F(s) e (s) ] * 


(2.51) 
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This equation is crucial to the conic sector results of Chapter 3. 

Note that (2 . 51) does not define a transfer function. Only LTI operators 
can be represented by transfer functions. 

One of the properties of the hybrid compensator that distinguishes 
it from LTI compensators is that the hybrid compensator spreads out the 
power spectrum of the input. This property is due to the sampler, which 
when viewed in the frequency domain shifts and adds the power spectrum 
of the sampled signal. An example is shown in Figure 2.5. The plots are 
magnitude versus frequency for signals at different points in the hybrid 
compensator. The input signal is bandlimi^ed, but the output signal 
has energy outside of this bandlimited region. A LTI compensator would 
have energy only in the same bandlimited region as the input. 

2.3.3 The Hybrid Feedback System 

The hybrid compensator is one part of the hybrid feedback system 
of Figure 2.6a. The plant G(s) is the same as in the analog feedback 
system in Section 2.2.1 and Figure 2.1. As with the analog feedback 
system, closed loop properties are determined by how signals pass 
around the loop. The difference is that loop transfer operators must 
be used instead of loop transfer functions . 

Consider the loop broken at point (1) in Figure 2.6a. Inject the 
input signal e^ n and let e Qut be the signal that returns after passing 
around the loop. Their Laplace transforms are related by 

e ( s) = G (s) H (s) D* (s) [F(s) e. (s) 3 * (2.52) 

—out — — — — —in 

This trans format ion can be represented by the loop transfer operator 
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Figure 2.6c; The digital feedback system 
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Sout 


e, 

~xn 


(2.53) 


but it is not possible to define a loop transfer function. The same 
is true for the loop broken at point (2) in Figure 2.6a. 

The conventional analysis of hybrid feedback systems avoids the 
problems of dealing with LTV operators by analyzing the system at points 
(3) and (4) of Figure 2.6a. Here the loop transfer operators T and 
are linear shift invariant and can be represented by the z-transforms 
T 3 (z) and T^z). This simplifies the analysis, but at the cost of only 
examining the system at the sample times. 

The conventional analysis proceeds by transforming the hybrid 
feedback system to a discrete feedback system. An intermediate step 
is shown in Figure 2.6b. The block diagram maniupulation used to derive 
Figure 2.6b is to pass the prefilter and sampler backwards across the 
sum. The prefilter and hold can be grouped with the plant to form the 
discretized plant 


G^ls) = [F (s) G(s) H(s)]* 


(2.54) 


which can also be represented by the z-transform G^(z) . 

The discrete portion of Figure 2.6b is extracted to form the 

digital feedback system of Figure 2.6c. The output y,(z) and the 

— d 

error e^Cz) are related to the command input £^(z) by 


-1 


y d (z) = G^z) D(z) U+£^(z)D(z)] r^ (z) 


(2.55) 


e^Cz) = U + (^(z) D(z) ) 1 r^z) 


(2.56) 
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It is possible to confute the analog output y(s) given a command 
input r(s) : 

y(s) = G(s) H (s) D^(s) [F(s)r(s)]» (2.57) 

where D*^(s) = D*(s) [I + GjMs) D*(s)] _1 (2.58) 

The closed loop operator is a linear time varying operator with the 
same structure as the hybrid operator defined by (2.51). 

2.3.4 Command Response, Stability, and Robustness 

The discrete loop transfer function and the discrete return dif- 
ference equation can be used to analyze the command response, stability, 
and robustness of the digital feedback system of Figure 2.6c. These 
results are analogous to those for the analog feedback system of 
Figure 2.1. One of the differences is that the results in this sub- 
section apply to discrete sequences - not to the analog signals - that 
appear in the hybrid feedback system. 

The command response is considered to be good if the discrete out- 
put y^(nT) tracks the discrete input r^(nT) with small error over some 
frequency range. As seen from (2.56), this will be the case if the 
return difference equation (similarly the loop transfer function) is 
large over the frequency range where the discrete input has significant 
energy. A command response specification can be stated 

°min D* ( jui) ] >_ p(cd) for u o < u < o) 1 


(2.59) 
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The loop transfer function is periodic, and the specification should 
only be given over a portion of the frequency range 0 to tt/T rad/sec. 

A discrete sequence r,(nT) is said to have significant energy over some 

“U 

frequency range if r^(z) evaluated at z * e^ WT has large magnitude for 
0 ) in the specified frequency range. 

The digital closed loop system is stable if all of the digital 
closed loop poles have magnitude less than unity. The digital closed 
loop poles can be found by various frequency domain and state space 
techniques (4] . 

Another way to determine stability is to use the discrete version 
of the multivariable Nyquist theorem. It will be stated here, using 
the same notation as in [31. Let N(ft, f(s), C) denote the "number of 
clockwise encirclements of the point ft by the locus f(s) as s traverses 
the closed contour C in the complex plane in a clockwise sense" [3, p. 76], 
Let the contour V be the unit circle, with small expansions to include 
the open loop poles of G^fzjDfz) with unit magnitude (these are con- 
sidered open loop stable) . Let P be the number of open loop unstable 
poles of GjfzjDtz). The multivariable Nyquist theorem states that the 
discrete closed loop system is stable (has no poles with magnitude > 1) . 
if and only if 

N (-1, -1 + det [I + G,*(s)D*(s)l , V) = - P (2.60) 

— a — 

An important point to note is that the hybrid feedback system is 
stable if and only if the discrete feedback system is stable. Hence, 
the two stability tests just mentioned are useful for the hybrid feed- 


back system. 
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The model of the plant is not exact, and closed loop stability 
must be preserved for all possible perturbations of the plant. Robust- 
ness results similar to those for analog systems can be derived. The 
starting point is a z-transform description of the additive perturbation: 

G d (z) = G^U) + E^z) (2.61) 

It is assumed that G^z) and G^tz) have the same number of open loop 
unstable poles, and that the magnitude of (z) is bounded by 

< £, a (U> for 0 < U) < W (2.62) 

-max —da — da — — s 

The digital closed loop system is stable if it is nominally closed loop 
stable and if 

C [E* D* (I + G*D*) -1 ( jo)) J <1 for 0 < w < U) 

max —da — — — d — — — s 

<= a [E*(j(i>)l < o ID* ( jui) + G*(ju>)] for 0 < (D < w 
max —da nun — — u — — s 

(2.63) 

The singular value inequality need only be checked over the fundamental 
frequency range. 

A digital robustness analysis can also be performed using the 
discrete multiplicative perturbation: 


(2.64) 
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It is assumed that G , (z) and G,(z) have the same number of open loop 

— u — 

unstable poles, and it is assumed that 

a [E*(j<o)l < £. (to) for 0 < a) < a) (2.65) 

max —am — cun — — s 

The closed loop system is stable if the nominal system is closed loop 
stable and if 

a [E* G* D* (I + G* D*)~ 1 (jU))] < 1 for 0 < W < 0 ) 

max —dm — d — — — d — — — s 

(2.66) 

<*= a (E . * ( jto) ] < a . [I + (G,* D*) -1 (jw)l for 0 < to < w 
max —dm min — — d — — — s 

(2.67) 

A digital robustness analysis starts with a discrete perturbation. 
This is not a natural place to start, however, for a robustness analysis 
of a hybrid system. It is the analog plant G(s) that is uncertain, 
and its uncertainty is expressed by an analog perturbation. The analog 
perturbation must be discretized in order to obtain a discrete pertur- 
bation. 

Consider when the actual plant is G(s) = G(s) + E (s) . The actual 

& 

discretized plant is 

GjMs) = G^Ms) + Ej a *(s) (2.68) 

where E. Ms) = (F E H(s}]* (2.69) 

—da — — a - 

The additive perturbation is discretized the same way as tne analog plant. 
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When multiplicative perturbations are used then the actual plant 

is G(s) * G(s) [I + E (s) ] and its discretized version is 
— — — — m 

G,*(s) - Gj*(s) + [FG^Hts)]* (2.70) 

This cannot easily be expressed as a discrete multiplicative pertur- 
bation except when E is a constant that commutes with H(s) 1 , in which 

-in — 

case 


GjMs) . 


'3a* (s) [I + E. ] , where E, 
-a — —ant 


(2.71) 


The analog and discrete multiplicative perturbations are equal. 

Constant E^'s can be used to find guarenteed phase and gain margins 
(se^ Subsection 2.2.7). By the argument of the above paragraph the phase 
and gain margins of the digital system are also phase and gain margins 
of the hybrid system. 

Discretizing analog perturbations is one way to analyze the ro- 
bustness of closed loop hybrid control systems. This is not, however, 
the approach that is persued in this thesis. Rather, the approach is 
to approximate the hybrid operator by a LT1 operator and then to use 
analog techniques to analyze robustness (Theorems 3.4, 3.5, and 3.7). 

The "analog techniques" are conic sector techniques. They are 
now reviwed, using a precise mathematical format. 


For SISO systems any constant E = e satisfies this condition. For 
multivariable systems the easiest example i*? when both E and H(s) are 
diagonal matrices. 
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2.4 Results from General Feedback Theory 

The hybrid compensator is modelled with a LTV operator which is 
called a hybrid operator. The analysis techniques described so far 
make no direct use of this hybrid operator. In order to do so, 
this thesis makes use of a feedback theory that is more general than the 
analog or digital feedback theory of the previous v.’o sections. Components 
of a feedback system are modelled with mathematical entities called 
"relations. " Hybrid compensators and the LTI plants are special cases 
of relations, so analysis techniques for general feedback systems can 
be applied to hybrid feedback systems. 

The description given here of the general feedback system follows 
the work of Zames, Safonov, Athans, Desoer, and Vidyasagar [5] to [9] . 

The most general part of the description is contained in the subsections 
on relations and conic sectors. The subsections on the feedback systems 
ars less general because they assume that the feedback system is causal 
and well-posed.^ Even with this restriction the components can be 
nonlinear and time varying. It is not until specific conic sectors are 
developed that further restrictions are needed. 

2.4.1 Extended Normed Linear Spaces 

The analog signals in a feedback system are members of an extended 
normed linear space. A "linear space" is a basic concept of analysis. 
Definitions and properties can be found in many textbooks [e.g. 19], A 
"normed linear space" is a linear space with a norm defined for it. 

The norm introduces the c _>ncepts of "size" and "distance”. Elements 

* A feedback system is causal if its output at time t is independent 
of its input after time t. It is well-posed if fer every possible 
input there exists an output. 
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of a normed linaar space must have finite norm. This restriction can 
be removed by extending the normed linear space. The result is called 
an "extended normed linear space." 

In this thesis the only normed linear space that will be extensively 
used is *- 2 °' the space of square integrable n-dimensional functions. 
Elements of L™ are functions x: R + ■* R n (from the set of real numbers 
> o to the set of n-dimensional vectors) that have finite norm. The 
norm is defined 





A 


r 

1 o 



[ 2 . 12 ) 


L has engineering significance as being the set of signals with finite 

energy. The square of the norm, || jc||? is proportional to the energy of 

L 2 

the signal x(t). 

The extension of L,, 11 is the extended normed linear space L 

2 2e 

Elements are functions x: R *► R n that have finite truncated norm II x |i 

~ + " ■» " c 

for all T 6 R + , where 

T 2 1 1//2 

||x(t) || E dt I (2.73) 

Elements of L ^ a re automatically in the extension and have the 

property that the limit as T -*• 00 of || x |! " || x || . . Examples of 

T L 2 

functions that are in L _ but not L „ are x(t) = t and x(t) = exp(t) . 

ze 2 

Examples of functions in neither space are x(t) = tan(t) and x(t) = 
i/d-t 2 ) . 
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2.4.2 Relations , Operators, Gain, and Stability 

A relation H is any subset of the product space x L 2 ™. 
range and domain of H are defined by 

Ra(H) A (y| (x,y) € H for some x 6 L_ n } (2.74) 

— ~ - - - 2e 

Do(K) A { x | (x,y) € H for some y € L_ m } (2.75) 

The inverse of the relation H is a set with the elements x and y arranged 
in the reverse order. This inverse relation always exists and is defined 
by 

H 1 £ {(y,x) 6 L™ x L 2 " j (x,y) € H} (2.76) 

The composition product HK and the sum H + K are relations defined by 

UK A { (x,z) € L n x L r | there e: cs a y € L™ such that 

— - ~ 2e 2e ~ 2e 

(x,y) 6 K and (y,z) G H) (2.77) 

H + K A { (x,y) 6 L n x | x 6 L_ n and y = y + y for some 
= ~ i 2e 2e 2e ~ ~1 i2 

y L G Ra(tf) and y 2 G Ra(fO) (2.78) 

An operator H is a relation which satisfies two conditions: 1) 

Do (H) =L„ n , and ?' for each x € L „ n there exists a unique y € L™ such 

2e - 2e i 2e 

that (x,y) 6 H. The same notation is used for both relations and 
operators. It is usually not important to distinguish between them, 
and when it is, it is usually apparent. For the thesis, it would not 
have been necessary to define relations, except that the inverse of 
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operators are not necessarily operators. 

The relation H cam be considered to be a transformation from the 
input space to the output space Lj™, in which case the notation 
y - Hx indicates that (x,y) € H. This transformation is somewhat 
dangerous to use, however, because there may be none, one, or several 
y’s for which (x,y) e H. If H is an operator then the transformation 
is well defined (it exists for every x e L. n and is unique) . 

The gain (or norm) of the relation H is defined by 


I|h|| l a 
2 


sup 


IIMI X 

!M1 t 


(2.79) 


where the suprenum is taken over all nonzero x e Do (H) , all corres- 
ponding Hx € Ra(tf), and all T e R + , In other words, for all possible 
input-output pairs and for all possible truncations. Note that the 
rat'o is always finite, and only the suprenum can be infinite. 

The relation H is defined to be L^-stable if |jff ||^ < “, in 

which case there exists a constant k such that 


IlHx |I T < k || X || T (2.80) 

for all x e Do(tf), all corresponding Hx 6 Ra(H) , and all T e R + . This 
type of stability is usually called "bounded input bounded output" 
stabi i icy. 

It may not be immediately apparent why an extended spaces are 
needed. The reason is that unstable : nations cannot be defined on 

the unex tended product space L ^ x A relation that is L ^-stable 

maps L n into L™ , but an unstable relation maps into l 2 ™. There- 
fore, to consistently define L .^-stable and unstable relations it is 
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necessary to use extended normed linear spaces. 

2.4.3 Conic Sectors 

In this section the concept of conic sectors is defined and dis- 
cussed. Necessary and sufficient conditions for relations to belong to 
a conic sector are presented. The reason for using conic sectors will 
become apparent later - when sufficient conditions for closed loop 
stability are sta ced in terms of conic sectors. 

Let H be a relation, and let C and R be operators. If 

II y - c x|l T < IIM! t - e ||x|l T (2.8D 

for all (x,y) fc H, T e R , and for some £ > o then H is said to be 
"strictly inside the conic sector with center C and radius R; M which 
is equivalently stated "st.ict'.y inside cone (C,R)." If (2.81) is 
true for some e > o then H is "inside cone <C,R) . " 

Now turn around the inequality sign. If 

i! y - c ? II ; * II II* + G IU W 2 X ( 2 . 82 ) 

for all (x,y) € H, T 6 R + , and for sane £ > o then H is "strictly outside 
cone (C,R) and if (2.82) is true for seme £ > o then H is "outside 
cone (C, R)." 

The easiest visual example of conic sectors is obtained from 
relations H which are memoryless nonlinear operators y = h(x). For 
example, consider the function y = h(x) plotted in Figure 2.7. 
function is bounded by lines with slopes c-r and c+r, which can be used 
to show 


y - cx 


(2.83) 
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Figure 2.7: A memoryless nonlinearity inside conic sector (c,r) 
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This inequality implies (2.81) with £=0, and therefore H is inside of 
cone (C,R) , where C and R are memoryless linear operators defined by 
y = cx and y = rx . 

One way to interpret H being inside of a cone (C,R) is to think 
of the center C as an approximation of H and to think of the radius R 
as a bound on the errors due to this approximation. Presumably if the 
approximation is to be of any use then R must in some sense be small. 

At the very least it should be the case that for a certain class of 
inputs x 6 S that ||Rx |j^ << J] Cx || . 

Another way to interpret conic sectors is to think of them as a 
bound on the energy of various signals. Consider the case* - when for 
all x e L 2 

ll(ff-C) X II , < !|Rx|L 2 (2.84) 

2 ~ L 2 

The L 2 norm is a measure of energy, so this inequality states that the 
energy of (tf-C)x is less than or equal to the energy of Rx. Let 
x € S be a set of signals that are in some way special. For example, 
inputs that have > 99% of their energy below 10 Hz. Then the center 
is a good approximation of H if the energy of (H-C)x is small for all 
x € S. One possible way to be quantitative about the approximation is 
to require for all x 6 S that Rx has < 1% of the energy of Cx. 

Two lemmas are now presented that give necessary and sufficient 
conditions for an operator to be either inside or outside of a cone. The 


The assumptions for this case are that W,C, and R are operators, 
H-C is stable, R is stable, and H is inside cone (C,R) . 


1 
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conic sector inequalities (2.81' and (2.82) are satisfied if and only if 
certain composite operators have gain < 1. Safonov 1 states the results 
[9, p. 42Cj , which to him fall into the "it can be shown" category. Here 
the attitude is not so cavalier, and proofs are included in the Appendix 
to Chapter 2. 

Lemma 2,1 Let H, C, P, and R 1 be operators such that H-C, R, and 
R 1 are (^“Stable. Then the following are equivalent: 

(a) H is strictly inside cone (C,R) 

▼ 2 

(b) || (H-C) R 1 II, < 1 - C for some G > o ■ (2.85) 

2 

Lenina 2,2 Let H, C, R, R 1 , and (I + CH) 1 be operators such that 
R, R 1 , and H(I + CH) 1 are L^-stable. Then the following are equivalent 

(a) -H 1 is outside cone (C,R) 

(b) ||RH (I + CM) 1 II ^ <1 • (2.86) 

2 

The radius R is a very special type of operator. In these two 

lemmas, and everywhere else in this thesis, it is assumed that R and R 1 

are (^“stable operators. The assumption that they are operators is 

enough to imply a one-to-one mapping between functions in the domain 

and range of R. The additional assumption that R and R 1 are L„ -stable 

2e 

implies a one-to-one mapping between the finite norm parts of the 
domain and range. 

What all of this means for LTI operators is that the Laplace 


1 


See also Zames (5] and Desoer and Vidyasagar [8] . 
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transform matrix r(s) is proper* and have no poles or zeros with real 
parts > o, and in addition, R *(s) is proper and has no poles or zeros 
with real parts > o. Examples are r(s) = 1, r(s) = {s+l)/(s + 10), 
and r (s) = (1 - .5e" sT )/(l - .9e~ ST ). 

One of the subtle interpretations of Lemma 2.2 is that being outside 

2 I 

of a cone is an inherently closed loop property . If ~H is outside 
of cone (C,R) then it is not useful to think of C as an approximation 
of -W . It is better to think of C as any operator that stabilizes 
the feedback system 


y = Hx v 

l (2.87) 

x = u - Cy \ 

The assumption that (I + CH ) 1 is an operator guarantees that the 
feedback system is well-posed (see Willems [20]). The additional 
ass um ption that H( I +CH)* is stable is another way of saying that the 
feedback system of (2.87) is closed loop stable. 

2.4.4 The General Feedback System and the Small Gain Theorem 

The general feedback system is shown in F’ jure 2.8. It will be 
referred to as System 1. The equations that define System 1 and the 
assumptions that the components satisfy are 


u = Ke 


§ = I - Gu 
E. e . u e L * 

G, K, and (I +GK ) 1 are causal operators / 


1 R (s) is proper if o max [R( 00 )] < ®, ie # poles 
not roll-off or grow as oj -*■ 00 . 


( 2 . 88 ) 


= # zeros, ie R(jo>) does 


In contrast to outside conic sectors being closed loop properties, 
inside conic sectors are inherently open loop properties. 
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Figure 2.8: System 1, the general feedback system 
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Any mathematical model that represents a physical system should be 
well-posed and causal. The assumptions in (2.88) guarantee this. The 
stability theory of Zames 1 5 , 6 ] and Safonov \1 , 9] does not require 
well-posedness and causality. These assumptions are made here because 
the extra level of generality is not needed to analyze hybrid control 
systems. 

The closed loop operators of Sysiatn 1 are E and U, where 

E A {(r,e) [ (r,e) € L r x L_ r and there exists u such that 
- ~ ~ ~ ~ 2e 2e ~ 

(2.88) is satisfied} (2.89) 

and where U is defined in a similar way. The closed loop system is 

stable if E and U are I^-stable. Hence, closed loop stability is proved 

by showing that | { E j[ . and |jU |j . are bounded. 

L 2 L 2 

The Small Gain Theorem is used to show that a certain cl^ss of 
systems is closed loop stable. The proof differs little from those of 
[5, Theorem 1] and [8, Theorem III. 2.1]. 

Lemma 2.3 (The Small Gain Theorem) Consider System 1. If 

II* II, • II G II i < 1 
2 2 

then E and U are -stable. • 

2e 

The Small Gain Theorem has a nice interpretation using the SISO 

Nyquist Theorem. Assume GK is SISO, LTI , L_ -stable, and that 

2e 

II GK ||. = max |gk(j«)| < 1 (2.90) 

L 2 6J 

The sufficient condition of the Small Gain Theorem is satisfied and 


therefore the closed loop system is stable. The Nyquist plot is inside 
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the unit circle and cannot encircle the (-1, 0) point. Hence, by the 
Nyquist Theorem the same conclusion 1 b reached that the closed loop 
system is stable. 

Unfortunately, any system that satisfied (2.90) is of little 
practical use, because good command response (and other performance 
specifications) generally require that |gk(jaj) | » 1 over some frequency 
range. The Small Gain Theorem would likewise be of little practical 
use if it were not for the Loop Transformation Theorem, which is 
discussed next. 

2.4.5 Sufficient Conditions for Closed Loop Stability 

If the operators G and K of System 1 satisfy certain conic sector 
conditions then the system is closed loop stable. This result is an 
extension of the Small Gain Theorem. 

The extension ir made possible by the Loop Transformation Theorem. 

A system is transformed to another system in such a way that the stability 
of the transformed system implies the stability of the original system. If 
G and K satisfy certain conic sector conditions then by the Small Gain Theorem 
the transformed system is stable, and by the Loop Transformation Theory 
the original system is stable. 

This approach is due to Zames [5, 6] . His results are less general 
because the cone center and radius are constant multipliers. The 
generalization presented here, which allows the center and radius to 
be operators, is due to Safonov [7] . He goes even further in his 
generalizations, and the main results presented in this preliminary 
chapter (Leninas 2.5, 2.6, and 2.7) are special cases of [7, Theorem 2.1]. 
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The transformations of System 1 are shown in Figure 2.9a. After 
a few block diagram manipulations the result in System 2 of Figure 2.9b. 
The defining equations and the assumptions for System 2 are: 


u = K 

e* + 

Cr 


~2 2 

~2 

•V 


e_ <= Rr 

- <?, 



-2 

2 

-2 


r, e. 6 

l r . 

u^ 

6 L “ 

~ -2 

2e 

~2 

2e 

g 2 . K 2 , 

and 

(I + 

<W 


(2.91) 


The transformations from System 1 are described algebraically by: 


k 2 = (K-o r 


G 2 = R G (T +CG) 


u 2 = (I + CG) u 


(2.92) 


?2 = 


The closed loop operators E_ and U respectively relate r to e. and u . 


Lemma 2.4 (Loop Transformation Theorem) 

(a) r, u, and e are solutions of System 1 if and only if r, 
and e 2 are solutions of System 2. 

(b) Furthermore if (I + CG) 1 and R 1 are E. -stable then the 

2e 

stability of System 2 implies the stability of System 1. • 

The proof of part (a) differs 1 ittle from those of (5, Theorem 2i] 
and [8, Theorem 111.6.3]. Part (b) follows after a few steps from (2.92). 
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Figure 2.9b: System 2 
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Lemma 2.5 (Sufficient Conditions for Closed Loop Stability, K 

Inside of a Cone). Consider System 1. Let C ana R be operators such 

that (I + CO) 1 , G(I + CG) 1 , K-C , ft, and ft* are L_ -stable. Then 

2e 

System 1 is closed loop stable if 

(i) K is strictly inside cone (C,ft) 

(ii) -G 1 is outside cone (C,ft) • 

The proof is now sketched out. If the conic sector conditions 
(i) and (ii) are true then (by Lemmas 2.1 and 2.2) it follows that 

II*, Ilf * IlG J|. < 1-e for some e > o (2.93) 

L 2 L 2 

By the Snail Gain Theorem (Lemma 2.3) it follows that System 2 is closed 

loop stable, and by the Loop Transformation Theorem (Leimna 2.4) it 

follows that System 1 is closed loop stable. 

The major restriction of Lemma 2.5 is that K-C must be -stable. 

2e 

Whatever is placed inside of the conic sector must be bounded by a stable 
radius . 

Leimna 2.5 is a robustness as well as a stability result. The 
closed loop system will be stable not just for a particular K and G but 
for any K that satisfies (i) and any G that satisfies (ii) . 

Condition (i) of Leimna 2.5 stages that the compensator K is inside 
of a cone. Sometimes it is more convenient to think of the plant G as 
being inside of a cone. Sufficient conditions for closed loop stability 
can also be derived when this is the case. A different transformation 
of System 1 is required, as shown in Figures 2.10a and 2.10b; and a 
different version of the Loop '.ransformation Theorem is required, so 
that stability of System 3 in Figure 2,10b implies the stability of 
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Figure 2.10a: Alternate transformation of system 1 



Figure 2.10b: System 3 
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System 1. The final result is shown below. 

Lemma 2.6 (Sufficient Conditions for Closed Loop Stability, “K 1 
Outside of a Cone) . Consider System 1. Let C and R be operators such 
that (I + CK) 1 , K(1 + CK)* , G**C, R, and R 1 are L 2 e -stable. Then System 
2 is closed loop stable if 

(i) -K 1 is outside cone (C, R) 

(ii) G is strictly inside cone (C,R) • 

The major restriction of Lemma 2.6 is that G-C must be L^-stable. 
Just as is the case for Lemma 2.5, Lemma 2.6 is a robustness as well as 
a stability result. 

Some of the references for t V *> eneral feedback theory set up a 
version of System 1 that has symmetric inputs, thereby avoiding separate 
derivations for inside and outside conic sector conditions for closed 
loop stability. Symmetric inputs are not usual for hybrid control 
systems, however, and therefore for this thesis it was decided to include 
the two derivations. 

One more set of sufficient conditions for closed loop stability 
is included. Tl.e difference here is that the loop transfer operator 
T = GK is placed inside of a cone. The general feedback system can 
still be considered to have two elements - one is and the other is I. 

Lemma 2 . 7 (Sufficient conditions for closed loop stability. 
inside of a cone). Consider System. 1. Let C and R be operators such 
that (I + C) 1 , T - C, R, and R 1 are L 2e -s«-able. System 1 is closed 
loop stable if 

(i) T is strictly inside cone (C,R) 

(ii) I! R (I + C ) 1 li . < 1 , 

2 
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If condition (ii) is true then -I is outside of cone (C,R). The 

major restriction of Lemma 2,7 is that T, - C must be L_ -stable. 

1 2e 

2.4.6 Conic Sector Analysis of Command Response 

Conic sectors can be used to analyze the response to signals that 

enter the feedback loop. The only signal that enters the general feed- 

back system of Figure 2.8 is the command input r, which enters just before 
the prefilter. 

In this thesis we will be interested in the steady state response to 
commands that belong to a set of command signals S C L^. In particular, 
we will be interested in the steady state response of low frequency 
sinusoidal inputs. Transient errors (overshoot, risetime, settling time, 
etc.) are not analyzed in this thesis. 

The error signal is 

e = (I +GK) 1 r (2.94) 

and the steady state command response specification is 

li e!| T 

lim < q for all r C S (2.95) 

T-** II r|| T “ 

The constant "q" is called the "quality measure”.. By letting the transient 
time T approach infinity the transient errors become insign if icant , and 
therefore the quality measure is a measure of the steady state error. It 
is in general difficult to compute the quality measure. As will be shown, 
conic sectors can be used to find an upperbound for q. 

There are two approaches. The first is to place the loop transfer 
operator GK inside of a cone, and the second is to place the closed loop 
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operator (T + GK) 1 inside of a cone. The first approach is due to Stein 
118] 1 , and the second approach is new. 

Placing GK Inside of a Cone 

Assume that the operators GK and C are L_ ^-stable, and assume that 

2e 

GK is inside of the cone (C, R^) . Expand the error equation (2.94) as 
follows : 


PaG'L 

QUALITY 


e = (T + GK) 1 r 
- (I + C + GK - C) 1 r 
= (I + C) 1 (I + (GK - C) (I + C) 1 ! 1 r 

= (I + C) 1 ^ - (I + (GK - C) (I + C) 1 ] 1 (GK - C) (I + C) 1 ) r 

(2.96) 

Take the truncated norm of each side of (2.96), and then use singular 
value and conic sector inequalities to obtain 


e|l T < li (I + r|| T + 


II a + c^Wl 


1 - II \(T + C) ‘;j L 


2 — II ^(T + C) 1 r |, T 


for all r 6 L J and all T e R (2.97) 

~ 2e + 


Stein [18] '“fines the quality measure by 


~ T 


INI 


q for all r € S and for all t 6 R + 


This quality measure is valid for transient as well as steady state error 
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An upperbound for the quality measure is 


q < q (1 + — r ) for all r 6 S 
- o m* 1 


(2.98) 


where 


q < lim 

° - T -K» 


11(1 + C) 1 r || 


for all r 8 S 




m* 


= 1 - R, (I + C) 


s* = || (I + C) 


(2.99) 


R, (I +C) A r 


r,q < lim 

1 ° " T -*=o 


~ T 


for all r € S 


T 


The constant "q is the "nominal quality measure/ 1 which is obtained 
o 

s * 

when the operator GK is replaced by C. The constant (1 + — r, ) is the 

m* l 

amount by which q^ must be raised to obtain an upperbound for q. If the 
set 5 is sinewaves with frequencies less then co^, and if C and R^ are 
LTI operators, then the consta ts in (2.99) are 


q = max 0 [(I + C) (jco)] 

o max 

O<C0<C0 

o 

r.* = 1 - max 0 [R (I + C) ^ ( joi) ] 
max 1 
CO 


( 2 . 100 ) 


s* = max a [ (I + C) (jco)] 
max 

CO 


r q = a (I + C)" 1 v ] 

1 c o«o<co max 1 
o 


-87- 


Two restrictions to the use of (2.98) are (1) that GK must be open 
loop stable and (2) that the cone around GK must meet the sufficient 
conditions for closed loop stability.^ in addition to these restrictions 
there is also the problem that (2.98) may give a conservative upperbound 
for q. By "conservative" it is meant that the upperbound for q is higher 
than it needs to be. If the cone around GK ii not very tight, then m* will 
be small and r^ will be large - both of which lead to conservativeness. 

Two more sources of conservativen^ss are that both m* and s* use operator 
norms, and hence must be t ue for all input signals and not just y 6 S 

Despite the problems mentioned above, equation (2,98) is a natural 
generalization of LTI results - specifically (2.21) and [1, equation (19)1. 
It is a first attempt at using conic sectors to analyze command response 
(as opposed to the more common uses of conic sectors to analyze closed 
loop stability and robustness) . One advantage of working with GK instead 
of (I + GK) 1 is that the perturbations of G and K are easier _o separate. 
For example, a multiplicative perturbation of G is also a multiplicative 
perturbation of GK, but not of (7 + GK) 1 . 

Placing (I + GK) 1 Inside of a Cone 

Assume that the operators (T + GK) 1 and (I + C) 1 are t- 2 e -sta k^ e ' 

I 12 

and assume that (I + GK) is inside of the cone [ (I + C) , R^J • Expand 


^It must be true that j] R (I + C ) 1 1 j f < 1, 

1 L 2 

condition for closed loop stability given 
equality is not satisfied then m* < o. 


which is the sufficient 
in Lenna 2,7. If this in- 


2 

The operator G is an approximation of GK. The open loop stability 
assumptions of C and GK are replaced by the less restrictive assumptions 
of closed loop stability. 
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the error equation (2.94) by adding and then subtracting (I + C) 1 rj 
e * (I + GK) 1 r 

= (I + C) 1 r + [(I + GJO 1 - (I + C)] r (2.101) 

An upperbound for the quality measure is 

q < q (1 + r.) for all r € S (2.102) 

~ o z. ~ 


where 


i (I + C) r I 

q * lim for all r G S 

°-T~ ||r|l 

•V i 


(2.103) 


(1+ GK) 1 - (I + 


r 2 q o < lim 

J-+G0 


C ) 1 r|| T 


for all r € S 


(2.104) 


~ T 


As before, the constant q is the nominal quality measure. The constant 

(1 + r ) is the amount that q must be raised to obtain an upperbound tor 
2 o 

q* 

One way to compute the constant r 2 is to construct a cone [(I + C) 1 , 

R 2 ] that contains (I + GK) 1 . It follows that 

II R 2~^ t r 

r_q < — — — for all r 6 l. and all T 6 R (2.105) 

2° - 2e 

Equation (2.105) may result in a conservative (i.e. large) estimate for r 2 q o 
because (1) the conic sector does not take advantage of the restricted 
set r € S and (2) the conic sector must be valid for all possible trunca- 


tions of the input signal. 
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An alternate form of (2.104) is 


[ GK (I + GK) 1 - C (I + C) 1 ] r 


r„q < lim 
2 o - _ 

J-KO 


r 

~ T 


for all r e S 


(2.106) 


which is derived from (2.104) by use of the identity 


(I + A) 1 - I - A (I + A) 1 


(2.107) 


This alternate form is easier to use for hybrid operators. It is 
interesting to note that conic sectors that contain (I + GK) 1 and 
GK (I + GK) 1 have the same radius. 

In this subsection two different approaches have been discussed for 
using conic sectors to analyze the steady state response to commands. 
Either the loop transfer operator GK or the closed loop operator (I + GK) 
can be placed inside of a cone. 
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2.5 Conic Sector Analysis of Analog Feedback Systems 

The general feedback theory is just that - general. It is only 
when specific assumptions are made about the feedback system that the 
results from the general feedback theory can be applied. 

In this section the conic sector analysis techniques are applied 
(for better or for worse) to the analog feedback system of Section 2.2. 

The feedback system is briefly described, sufficient conditions for 
LTI operators to be inside or outside of a cone are presented, and then 
sufficient conditions for robust closed loop stability are presented. 

The conic sector analysis techniques are not as general as the analysis 
techniques of Section 2-2. The distinction is that conic sectors cannot 
be used for unstable perturbations of the analog plant. 

The compensator and nominal plant are modelled, respectively, by 
K and G. The uncertainty of the nominal plant is modelled by the 
additive perturbation E , where 

G = G + E (2.108) 

a 

and where E is a L -stable LTI operator . The only other information 

ci Z 6 

known about E^ is a bound on the magnitude of its Fourier transform 
matrix : 

0 [E (ju>)] < l (to) for all to (2.109) 

max -a - a 

The uncertainty of the nominal plant can also be modelled by the 

multiplicative perturbation E , where 

m 

G = G(I + E ) (2.110) 

n» 

0 [E <jto)l < l (to) for all to 
max -m m 


< 2 . 111 ) 
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and where E = GE is L, -stable.^ 
a m 2e 

Three Lemmas are now presented: Lemma 2.8 gives sufficient condi- 
tions for G to be strictly inside of a cone, Lemma 2.9 gives sufficient 
conditions for - K 1 to be outside of a cone, and Lemma 2.10 combines 

the previous two and gives sufficient conditions for robust closed loop 
2 

stability. 

^ I ‘X 

Lemma 2.8 Let G, G, R, and R be LT1 operators such that G-G, R, 
and R 1 are L^-stable. G is strictly inside cone (G , R) if 


a . [R ( joo) ] > 

mm - - 


(l-e) 


1/2 max - 


[G(jCJ) - G(jw)] 


( 2 . 112 >) 


for all gj and some £ > o 

Lenuna 2 . 9 Let K, G, R, and R 1 be LTI operators such that 

K(J + GK) 1 , R, and R T are -stable. - K 1 is outside cone (G,R) if 

2e 


O [RK (I + OK) ~ 1 (iCJ)] < 1 for all GJ 
max ™ " 


■( 2 . 112 ) 


Lemma 2.10 Consider the analog feedback system containing the LTI 
operators K, G, and G. Let R and R 1 be i- 2 e _sta bl e LIT operators. The 
closed loop system is stable for all possible G's if there exists an 
R(jco) such that 


The requirement that E is (^-stable means that the exact number and 
the exact location of unstable poles of G(s) must be known. The less 
restrictive requirement in Section 2.2 is that G(s) and G(s) must have 
the same number of uns .able poles, which means that only the exact 
number of unstable poles must be known. 


Lemmas 2.8, 2.9, and 2.10 (for analog systems) correspond, respectively, 
to Theorems 3.?, 3.C. , and 3.7 (for hybrid systems). Lemmas 2.8 and 2.9 
are due to Safonov [9, Lemas A4 and A2). He gives necessary and 
sufficient conditions to be conic. Only the sufficient conditions are 
presented here, to facilitate comparison with the results for hybrid 
systems . 



92- 


(i) 

K(J + 

GK) 1 is -stable 

2e 



(ii) 

^max 

[RK (I + GK)" 1 (jOJ)] < 1 for all o> 


(2.114) 

(iii) 

°min 

[R(jo»] > t (oj) for all 03 

- a 

■ 

(2.115) 


The proof for Lemma 2,10 can be quickly sketched out. Conditions 
(i)_ and (ii) guarantee that -K* is outside of cone (G , R) . Condition 
(iii) guarantees that G is strictly inside cone (G , R) . Hence, by 
Lemma 2.6, the analog system is closed loop stable , 

The first step in applying conic sectors to analog feedback systems 
is to construct a cone (G, R) such that all possible G's are strictly 
inside of it. By Lemma 2.8, this is true if 


o . 
min 


[R(j0)) > \~py 

(1 -c) 1 '* 


l Co» 

3 


for all 0 ) and some e > o 


(2.11S) 


The e term is a nuisance. It suffices to repace (2.115) with the strict 
inequality i 


C [R( jo)) ] > t ( 03 ) uniformly in 03 (2.117) 

If the uncertainty is modelled by a multiplicative perturbation then 
(2.117) is replaced by 

O . (R ( jo)) ] > t (oj) • a [G(j 03 )j uniformly in 03 (2.118) 

Note that the additive perturbation E^ cannot in general be used 

as the radius of the cone. A radius must have the property that both 

it and its inverse are t -stable. This is not assumed about E - only 

2e a 

that E is stable. The problem is that E ( jos) may roll off, whereas 

cl rft 

F(j03) must eventually flatten out (or become periodic). This distinction 
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is not cf practical importance, however, because the radius can be 

arbitrarily close to bound t (oj) out to arbitrarily high frequencies. 

A 
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Appendix to Chapter 2 

This appendix contains proof of Lemmas 2.1 and 2.2. 

Proof of Lanina 2.1 Here i'_ is shown that H is strictly inside of 
cone ( C, R) if and only if the operator (tf - C) R 1 has gain < 1 - e for 
some £ > o. The operators are defined on the cross product spaces: 


H, C C L * x l " 


R C L 2> ^2e 


W-C1P 1 C X L* e 


(2.A.1) 


Each step in the proof is equivalent, 
parenthesis. 

H is strictly inside cone (C,R) 


Explanations are enclosed in 

(2.A.2) 

(2. A. 3) 

and some E > o 


<-*■ II y - C;||* < ||Rx||* -e ||*||* 

for all (x,y) € H, all T € R + , 
(definition of strictly inside) 


o*> 


li(H-C) x[|* < 

for all x € 
(because H, 


IMI 2 -e II *11 2 

L„ n , all T 6 R , and sane £ > o 
2e + 

C, and R are operators with the same 


domains) 


(2. A. 4) 


<-> II (h-cjr 1 ell 2 < lie II 2 - c Hr 1 ! II 2 

for all E 6 U n , and T G R , and some £ > o 
£ 2e f 

where K £ Rx 


(2. A. 5) 


(because R and R 1 are operators) 



<■> 


<v»> 
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|i (H-C)R 1 g|j T 2 < i| £ || 2 a-e<) 


(2.A.6) 


for dll 


£ € all T e R + , and some e* > o 


where e l = a e and 


II^IL 

O < for all £ | ft o 

llcll T - T 


(a < “ because R is -stable) 

2e 


J y ll 

||(H-C)R || A -,u P < 1 ~ e' 

II ? ll T 


(2.A.7) 


suprenum taken over all £ e L with 


£ || / o and all T € R 

n T + 


(only possible if H-C and R are L^-stable) 


This complete the proof. 
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Proof of Lemma 2.2 Here it is shown that -H* is outside of a cone 
if and only if the operator RH (I + CH) 1 has gain < 1. The operators are 
defined on the following cross product spaces: 


-H 1 . CC L m x L n 
2e 2e 

R C L ® x L ® 

2e 2e 

RH(I + CH) 1 c L n x L® 
2e 


(2.A.8) 


Note that the center and radius differ from those of Lemma 2.1. Each 
step in the following is equivalent; 


-H 1 is outside cone (C,R) 


<~> j| x - Cy|| 2 > ||Ry || 


(2.A.9) 

(2. A. 10) 


for 


all (y,x) e -H and all T e R + 


(by definition) 


<=> || (I + CH) x || ^ > ||RHx|| 2 (2. A. 11) 

for all x 8 L n and all T e R 
2e + 

(because H is an operator) 

<=> li § li 2 >11 RH (I +CH) 1 | II 2 (2. A. 12) 

for all £ £ L_ n and all T 6 R 
-w 2e + 

(because (I + CH) and (T + CH) 1 are operators) 
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<-> pH (I + CH) T II 2 


sup 


|| RH(I + cm 1 5 1|| 
II 5 III 


(2. A. 13) 


n 2 

suprenum taken over all £ € L with || 5 II T ^ o 
and all T 6 R . 

T 

(only possible if R and H( I + CH ) 1 are L^-stable) 


This completes the proof. 
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3- CONIC SECTORS FOR HYBRID OPERATORS 
3 . 1 Introduction 


The major theoretical results of this thesis are presented in this 
chapter. It is here that the conic sector analysis techniques developed 
for general feedback systems are applied to hybrid feedback systems. 

The results are condensed into eight theorems, which are now briefly 
described: 


Theorem 3.1: An upper bound for the gain of the hybrid operator. 

Theorem 3.2: Sufficient conditions for a hybrid operator to be 

inside of a cone. 


Theorem 3.3: Sufficient conditions for closed loop stability 

when the hybrid operator is inside of a cone. 

Theorem 3.4: First approach to sufficient conditions for robust 

closed loop stability when the hybrid operator is 
inside of a cone. 


Theorem 3.5: Second approach to sufficient conditions for robust 

closed loop stability when the hybrid operator is 
inside of a cone. 

Theorem 3.6: Sufficient conditions for a hybrid operator to be 

outside of a cone. 


Theorem 3.7: Sufficient conditions for robust closed loop sta- 

bility when the hybrid operator is outside of a cone. 

Theorem 3.8: Use of the closed loop hybrid operator to analyze 

command response. 


The theorems about closed loop stability and command response 
(3.3, 3.4, 3.5, 3.7, and 3.8) all begin with the statement: "Consider 

the hybrid feedback system." This refers to the hybrid feedback system 
of Section 2.3, which is a special case of System 1 of Section 2.4. 


For the hybrid feedback system: 

• The compensator is modelled with the casual hybrid operator K, 
which relates the Laplace transforms of its input and output 
by u(s) = H (s) D* (s) (F(s)e(s)]*. 
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. The nominal and actual plants are modelled by the causual LTI 
operators G and G, which are related by either the additive 
perturbation E or the multiplicative perturbation E^, as 
described by equations (2.108) to (2.111). 

. Associated with the hybrid feedback system are the LTI operators 
C, R, and R 1 ; which form the center and radius of various conic 
sectors. It is always assumed that R and R* are l 2 e -stable. This, 
and any other assumptions of open loop or nominal closed loop 
stability, are explicitly stated in the Theorems. 

, The hybrid feedback system is assumed to be well-posed and causal, 
in other words it is assumed that (I+GK)* is a casual operator. 

For Theorems 3.3, 3.4, and 3.5 it is also assumed that (T+CG)* is a 
casual operator. 

A distinction is made between "closed loop stability" (Theorem 3,3) 
and "robust closed loop stability" (Theorems 3.4, 3.5, and 3.7). The 
former applies only to the nominal plant G, and the latter to all pos- 
sible plants G. A distinction is also made as to whether the hybrid 
operator is inside of a cone (Theorems 3.2 to 3.5) or outside of a cone 
(Theorems 3.6 and 3.7). When the hybrid operator is inside of a cone, 
then the center C can be viewed as a LTI approximation of K. On the 
other hand, when the hybrid operator is outside of a cone, the center 
C cam be viewed as an approximation of the plamt, and is usually set 
equal to the nominal plant G. 
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3.2 The Gain of the Hybrid Operator 

One of the properties of an operator is its gain [see (2.79)]. 

It is the maximum ratio of output norm to input norm. An upperbound 
for the gain of the hybrid operator K is presented in Theorem 3 . 1 . 
When K is SISO then an input signal can be constructed which achieves 
the upperbound, hence the upperbound actually is the gain. 

Theorem 3.1 Let K be a stable hybrid operator. 

(a) An upperbound for the gain of K is 



max T r 

°± W ± i ( T Ik 



He* II 



( 3 . 1 ) 


(b) Furthermore, if K is SISO then (3.1) is true with equality.* 
The proof is presented in the Appendix to Chapter 3. Included in 
the proof of part (b) is a signal which achxeves the upperbound. At 
the end of the proof a conjecture is made about the actual gain in the 
multivariable case. 

After some thought (possibly after considerable thought) , Theorem 
3.1 should be as intuitively reasonable as the corresponding result for 
a LTI operator H, which is 


W\ L " T II H ( ju»|| (3.2) 

2 

The input signal that achieves the maximum gain for LTT. operators is a 

sinewave at the frequency u) o that maximizes (3.2). The corresponding 

input signal for hybrid operators is an infinite series of sinewaves. 

The fundamental frequency is the (jJ q that maximizes (3.1) , and the other 

frequencies are shifted away from co by integer multiples of u> . 

o s 
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The gain of the hybrid compensator depends on how fast the pre- 
filter rolls off. Consider when the impulse response of the SISO pre- 
filter approaches an impulse, i.e. when f(ju>) 1. The infinite sum 
2 

E IjfJJ approaches infinity, which i\iplies that the gain of the hybrid 
n 

operator approaches infinity. Spurious inputs such as noise will be 
greatly amplified. Hence, Theorem 3.1 can be used to justify the need 
for adequate prefiltering. This duplicates a result of Kostovetsky [17] . 

The gain also depends on how fast the hold rolls off. In fact, the 
hold is treated exactly the same way as the prefilter. This results in 
the intuitively pleasing symmetry of equation (3.1). This symmetry is 
important when the hybrid compensator is included in a feedback system, 
in which case it is not obvious which parts of the loop should be in- 
cluded in the prefilter and which parts should be included in the hold. 
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3.3 A Cone that Contains the Hybrid Operator 
3.3.1 E xistence 

For any stable hybrid operator K there e :ists a cone (C,R) such 
that K is strictly inside of the cone. This result is presented in 
Theorem 3.2. 

Theorem 3.2 Let K be a stable hybrid operator, let C be any LTI 
/.^-stable operator, and let R and R 1 be LTI L^-stable operators. 

(a) K is strictly inside cone (C,R) if 


a . 

min 


R(jo))| j> 


(1-e) 


1/2 


r i " r 2 (U)) +r 3 (a)) 


1/2 


for all U) and some e > 0 


(3.3) 


where r^w) = ll*i k l| 2 ] * llE *!| 2 2 ] ( 3 - 4 > 

r 2 (u) = -L Z IlHjjD'rJI 2 (3.5) 

T k 

r 3 (w> = l II* - ^ll 2 (3.6) 

k 


(b) Furthermore, the optimal center 

C(s) = i H(s) D* (s) F (s) (3.7) 

minimizes the lower bound for C . (R( jw) ) . • 

min — 

The proof of this result is in the Appendix to Chapter 3. A few 
remarks are now made to highlight the theorem. A conic sector exists 
for any stable hybrid operator. The stability requirement is the major 
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limitation of this theorem, in fact it is the major limitation of all 
the inside conic sector results. Unstable hybrid operators (which in- 
cludes hybrid operators with digital integrators) cannot be bounded by 
finite radii. 

The conic sector is not unique - and not always useful. The choice 
of center determines to a large extent how useful the conic sector is. 

The center should be a good approximation of the low frequency behavior 
of the hybrid operator . 

Usually, but not always, the optimal center (3.7) is a good choice. 
It has the distinct advantage of minimizing the radius. Also, since 
this choice makes ^(w) - 0, there is one less term to compute. The 
optimal center is an infinite dimensional LTI operator. If a low order 
rational polynominal center is desired then use 

C(s) = h H (s) K(s > F < s > (3.8) 

X u “ ' 

where h (*) : s an approximation of a zero-order-hold, such as (2.47). if 
the objective is to compare different discretization techniques then use 

C(s) = i H (s) K(s) F (s) (3.9) 

The magnitude of the radius will rise or fall depending on how well D*(s) 
approximates K(s) at low frequencies. 

An alternate equation for the radius uses a double summation: 

°min [— * ^ | - [ r 4 <*» « 3 «■» ] ^ (3.10) 
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where 

r 4 (u>) = \ 1 1 ll^£*rjl 2 (3.11) 

T k nyflc 

and wher. (ca) is given by (3.6). For SISO systems (oj) (ui) (<*)) , 
and for multivariable systems r 4 (w) <r ]L (w) -r 2 (co) . Whether to use (3.3) 
or (3.10) to compute the radius is primarily a computational issue, 
as discussed in Section 4.3. 

Theorem 3.2 only gives a lower bound for the radius. It is always 
best to choose R(ju>) so that it actually equals the lower bound, because 
this choice gives the least conservative stability and robustness re- 
sults. If R(jw) is multivariable then all of its singular values should 
be set equal to the lower bound. This is done by making R(jto) a diagonal 
matrix with equal diagonal entries. 

An important point to make about the radius is that it is computable. 
In the examples of Chapter 5 the radius is computed for many different 
values of w and then plotted on a magnitude Bode plot. The radius is 

periodic with period w , and need only be computed over the frequency 

s 

range 0 to U) g /2. 

The radius R(ju)) corresponds to an additive perturbation, whereas 
the quantity R(ju))C ^"(jw) corresponds to a multiplicative perturbation. 
The quantity R C ^ ( ju>) is called the "multiplicative radius." It con- 
tains information about the size of the radix’s relative to the center, 
which is usually more important than the absolute size of the radius. 

The center C(jw) is a good approximation of the hybrid operator K over 

the frequency range where 0 |R C ^(ju))| << 1. 

max [ — — j 

If the hybrid operator is SISO and if the optimal center is used 
for the cone then the multiplicative radius is independent of the 
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computer transform d*(s). In this case the multiplicative radius de- 
pends only on the prefilter and hold. All of the closed loop stability, 
robustness, and performance results do depend (of course) on the com- 
puter z-transform. 

The last comment is that Theorem 3.1 is a special case of Theorem 3.2. 
If the choice of center is C(s)=Q, then r 2 (w) = r 3 (to) , which means that 
-r 2 (to) + r (lo) = 0* What remains of (3.3) can be used in place of (3.1) 
of Theorem 3.1. 

3.3.2 Closed Loop Stability 

The hybrid operator K is one part of a hybrid feedback system. 

Sufficient conditions for closed loop stability are now presented. 

Theorem 3.3 Consider the hybrid feedback system (see Section 3.1) . 

Assume K, C, R, and R 1 are L. -stable. The hybrid feedback system is 

2e 

closed loop stable if a C and R exist such that 

(i) K is strictly inside cone (C,R) 

(ii) Gtl+CG) 1 is L„ -stable 

2e 

(iii) o (r G(I + C G) 1 (ju))| < 1 for all to ■ (3.12) 

ntdx l "*■ “ ■ j 

The proof can be quickly sketched out. If condition (i) is true 
then K is strictly inside cone (C, R) . Such a cone always exists (by 
Theorem 3.2), because it is assumed that K is (.^-stable. If conditions 
(ii) and (iii) are true then -G 1 is outside of the same cone (by Lemma 
2.8). The fact that K is strictly inside and -G 1 is outside of the same 
cone implies that the hybrid system is closed loop stable (by Lemma 2.5) . 

Theorem 3.3 is applied by first constructing a C(jw) and R(jto) such 

that condition (i) is true, and then by using C(jco) and R(joo) to check 

conditions (ii) and (iii). The major restriction ot this theorem is 
that the hybrid operator must be open loop stable. This guarantees 
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the existence of a C(joo) and jl(jw) such that condition (i) is true. 

If the hybrid operator is not open loop stable then Theorem 3.3 cannot 
be used to determine closed loop stability. 

Condition (ii) is a check for nominal closed loop stability. The 
nominal system is an analog LTI system, and therefore standard techniques 
can be used to determine closed loop stability (see Subsection 2.2.5). 

The singular value inequality (3.12) can be graphically checked on 
a magnitude Bode plot. This is done in the examples of Chapter 5. The 
inequality (3.12) is implied by either of the following: 1 

O ( R( jo))l • O Ig (I + C G) _1 ( joj) 1 < 1 for all u> (3.13) 

max [- ] max[— — J — 

a [r C -1 (j(jL))l < 0 . 1 1 + (C G) _1 ( jt*>) I for all U) (3.14) 

It is (3.13) that is used in the examples. The problem with (3.14) is 
that it cannot be used if C(jw) = 0^ for some w. This happens whenever 
the hold H(juj) is a zero-order-hold and the center that is chosen in- 
cludes H(jw) . The only way that (3.14) can be used is if the center 
includes an approximation to H(jw), such as H (jw) of (2.47). 

— u 

3.3.3 Robust Closed Loop Stability 

One of the problems with Theorem 3.3 is that plant uncertainties 
are not explicitly included in the sufficient conditions for closed 
loop stability. The closed system should be stable not just for the 
nominal plant G but for all possible plants G in a defined set. 

1 Note the similarity of (3.13) and (3,14) to the robustness conditions 
(2.23) and (2.22) for analog feedback systems. 
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Two robustness results are now derived. The first result uses 
the fact that conditions (ii) and (iii) are true for a set of linear 
operators, of which G is one element. This set will contain all pos- 
sible G's if the condition of Theorem 3.4 is satisfied. 

Theorem 3.4 Consider the hybrid feedback system. Assume that all 
of the assumptions and conditions of Theorem 3.3 are satisfied. The 
hybrid feedback system is closed loop stable for all possible G's if 


l(u) < 

m 


a . [I + G C ( jw) ] - 0 [R G(jU3)] 
nun — ~ — m ax 

o v lG + a [R G(jw)] 

max max — — 


for all u 


{3.15} 


The proof of Theorem 3.4 is included in the Appendix to Chapter 3. 
Theorem 3.4 is set up to be used when the uncertainty of the nominal 
plant G is characterized by a multiplicative perturbation. The hybrid 
operator K is inside of the cone (C,R), and if the conditions of Theorem 3.3 
are satisfied then -G 1 is outside of this cone. If in addition the 
condition of Theorem 3.4 is satisfied then all possible -o's are out- 
side of this cone. 

The next robustness result uses a different approach. The per- 
turbation of the nominal plant is used to increase tna size of the 
radiur {the radius of the cone that the hybrid operator is inside of) . 

This is accomplished by grouping (£ + E^) , the multiplicative pertur- 
bation, with the hold.^ 

Theorem 3.5 Consider the hybrid feedback system. Assume that 

K. G, E , C, R, and R 1 are L_ -stable. Choose the center 
m 2e 

^"Similar results are possible by grouping the multiplicative pertur- 
bation with the prefilter. 
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rv 


C(s) - | H (s) DM») F(s) 


(3.16) 


(a) The composite operation (I+E m )K is strictly inside cone 


(C,R) if 


i 7^172 K ( “' * V U) * r 7 ( “>r /2 


(3.17) 


for all a) and some £ > 0 


where 

V"* - ^■[^ 1+£ rok )2 • IImJI 2 j • lloii 2 • 1 2 i|fJ| 2 j 0.18) 

r 6 (u) = \ £(1 + • IK D* F k || 2 (3.19) 

T k 

r 7 ( “> - l3 - 20 > 

T k 

(b) .he hybrid system io closed .loop stable for all possible 


G* s if 





(i) 

(I+E )K is 
m 

strictly inside cone 

(C,R) 


(ii) 

G(l + CG ) 1 

is -stable 

2e 



(iii) 

o Irg 

max l 

(I + C G)" 1 (ju))j < 1 

for all u) 

II (3.21) 


The proof of Theorem 3.5 is a composition of the proofs of Theorems 
3.2 and 3.3. The proof of part (a) differs from the proof of Theorem 
3.2 only in that H(s) is replaced by (1_ + E^)H(s), and the following 
inequality is inserted in the appropriate place in the proof: 
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1 ll<r%)|| • ||h( ju>) 1) 

1 (l + llEjl) “ llH(jw)l) 

1 (l + A*) * l|H(j“) II (3. 22) 

The proof of (b) differs from the proof of Theorem 3.3 only in that K 

is replaced by (I+E )K. The hybrid operator (I+E^JK is strictly inside 

of cone (C,R) if condition (bi) is true, and -G 1 is outside of the same 

cone if conditions (bii) and (biii) are true. Hence, the closed loop 

system is stable for all possible G's in the set defined by E . 

m 

An alternate expression for the radius (3.17) uses a double sum- 
mation : 


0 mn[£ l H > iV 1 *” + r 7 ( “>] 


1/2 


(3.23) 


for all u) and some e > 0 


where 


r g (u>) = r 5 (w) - r 6 (a)) 


3 l * (1 + v 2 -llvtji 

T k n?*k 


(3.24) 


and where (w) is defined by (3.20) 

Constant values of £ m (w) are useful for phase and gain margin 
analysis (see Subsection 2.2.’' Suppose 


£ ( w) -■= 

m 


a 


for all u) 


(3.25) 
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Then the various components of the radii of Theorems 3.3 and 3.5 are 
related by 


r 5 (w) 

= (1+a) 2 r 1 (a)) 

r 6 (u» 

= (1+a) 2 r 2 (io) 

r ? (w) 

= a 2 r 2 (w) 

r 8 (o» 

= (1+a) 2 r 4 (w) 


(3.26) 


The significance of (3.26) is that the same computations for Theorem 
3.3 suffice for Theorem 3.5. 

It is possible to systematically find the maximum value of a such 

that the robustness condition (3.21) is satisfied. Let £(<»>) = a and 

m 

substitute parts of (3.26) into (3.21) to obtain: 


r. 


,2 2 , 2 2 

L d**> r 4 + a r 3 


.1/2 


g ( 1+cg ) 1 < 1 


for all 0 ) 


(3.27) 


Change the inequality to an equality, and manipulate (3.27) to obtain 


a + 'xi + c = 0 


(3.28) 


where b and c are functions of r^, c, and g. If the stability con- 
dition (3.12) is true (i.e. if (3.21) is true when a=0) then c < 0 
and (3.28) has one negative and one positive real root. Let the posi- 
tive real root be a^tu), which is a function of u). The maximum value 
of a that satisfies (3.21) is 
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a = min a. (co) (3.29) 

a) 1 

Note that a^(u>) is minimized over w 

As just shown. Theorem 3.5, equations (3.28) and (3.29), can be 
used to find a guaranteed phase and gain margin. Minimizing over u 
the positive real root of the quadratic equation (3.28) is not diffi- 
cult, but it is nevertheless easier to use Theorem 3.4, equation (3.15) , 
to find a guaranteed phase and gain margin. As will be seen later. 
Theorem 3.7, equation (3.40) can also be used to find a guaranteed phase 
and gain margin. In general, Theorems 3.4, 3.5, and 3.7 will give 
different margins. 
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3.4 A Cone that Contains the Loop Transfer Operator 

The entire loop transfer operator can be placed inside of a cone . 

This will usually result in less conservative sufficient conditions for 
closed loop stability and robustness. The size of the radius depends 
on the amount of aliasing of both the prefilter and the hold. In- 
cluding the plant with the compensator will usually help because the 
rolloff of the plaint reduces the amount of aliasing of the prefilter 
and/or the hold. 

There are some pitfalls, however, the most important of which is 
where to include the plant. There are three choices: 

Case 1: with the hold 

Case 2: with the prefilter 

Case 3: with a combination of the hold and prefilter 

The three cases are illustrated in Figure 3.1. 

Case 1 corresponds to breaking the feedback loop before the pre- 
filter. The loop transfer operator at this point cam be used to amalyze 
the response to signals injected there (such as r) . 

Case 2 corresponds to breaking the loop before the plant. The 
loop transfer operator at this point can be used to amalyze the response 
to signals injected at this point. 

Case 3 does not correspond to breaking the loop at amy physical 
point. Both the prefilter amd the hold are replaced by (F^ G H)*^ 2 . 

2 

Case 3 is the best to use to analyze closed loop stability and robustness. 

1 Whatever is placed inside of a cone must be open loop stable. If the 
plant is open loop unstable it can be split into a stable part amd a 
unit magnitude unstable part, and then the stable part can be included 
in the cone. See Section 5.4 for an example. 

■> 

“Nowhere does the square root actually have to be computed, because 
|| (F G H) 1/2 || 2 = || F G H || . 
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Case 1 : Plant included with the hold 



Case 2 . Plant included with the prefilter 



Case 3: Plant included with a combination of the hold and prefilter 


Figure 3.1: Three cases of including the plant with the hybrid 

compensator. 
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Theorems 3.1 to 3.5 are used when the hybrid operator is inside 
of a cone. With appropriate modifications these theorems can also be 
used when the loop transfer operator is inside of a cone. Tie modifica- 
tions are described in the following paragraph. 

In each of the theorems the hybrid operator K is replaced by the 
appropriate loop transfer T. Depending on where the plant is included: 
Case i: Replace H by G H 

Case 2 : Replace F by F G 

1/2 

Case 3: Replace H and F by (F G H) 

In Theorem 3.3 replace G by I, and the three conditions for closed loop 
stability become 


(i) 

T is strictly inside cone 

(C,R) 


(ii) 

(I+C) 1 is L -stable 
2e 



(iii) 

<W 1 1 

for all w 

(3.30) 


In Theorem 3.4 replace the condition for robust closed loop stability 
by 


£«*)) < 
in 


o . [I + C(joj)] - a [r ( jto) ] 
nun — — max — 

a [c(jai)J + a [R(jo))] 
m ax — max — 


(3.31) 


In Theorem 3.5 replace G by I and replace the three conditions for robust 


closed loop 

stability by 


(i) 

(I+E )T is strictly inside cone (C,R) 
m 


(ii) 

(I+C) 1 is L -stable 
2e 


(iii) 

a |R(I+C) 1 (ju)) 1 < 1 for all w 

max j— ■ — J — 

(3.32) 
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The significance of case 3 is that for SISO systems the radius is 
smaller (for each w) than the radii for cases 1 and 2. Hence, case 3 
gives the least conservative sufficient conditions for closed loop 
stability and robustness. The Cauchy-Schwartz inequality is used to 
show this. 

The equation for the radius in Theorem 3.3 is 1 

r(w) = Jr^w) - r 2 (w) + r 3 (w)| 1//2 (3.33) 


For SISO systems r 2 (ui) and r^ (u) do not change for the three cases. 

Therefore, whatever differences exist in r(tt) cure due to differences in 
2 

r 1 (jw). The equations for r^U)) are: 


Case 1: r^lw) 

Case 2: r^lw) 


Case 3: r^w) 


s\i |d *' 2 •[; |f » |2 ] 

^M 1 ] •l*‘l a -l£l^l 2 ] 

|£ A \!] 2 • l d *i 2 


(3.34) 


By the Cauchy-Schwartz inequality: 


Without loss of generalty the e term has been ignored. 


For multivariable systems both r 2 (w) and r 3 (w) differ for the three 
cases, because matrices do not commute. 
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[H^i] 2 i(^K\i 2 ]-[Mg 2 ] 

> (3.35) 

[? l f A\l] 2 1 \l Kl 2 ] • [* l £ n 9 n' 2 ] 1 


Therefore r^ui) for case 3 is less (for each w) then r^((ii) for the other 
cases. It follows that r(o>) for case 3 is less (for each u>) then r(ti)) 
for cases 1 and 2. 

Without further information about h(s) and f(s) it is rot possible 
to tell which of the radii for cases 1 and 2 is t mailer. Some indica- 
tion of which radii is smaller can be gleaned by comparing the relative 

2 2 

magnitudes of — ijh, | and E |f [ , but the best way is just to graphically 
T 2 k * n n 

compare the radii. In general, neither will be smaller for all w. 



-117- 


ORIGINAL. PAvi- 
OF POOR QUAU 


3.5 a Cone that the Inverse Hybrid Operator Is Outside of 
3.5.1 Existence 

Just as there exists a conic sector that K is inside of, there 

also exists a conic sector that - K 1 is outside of. This result is 

presented in Theorem 3.6. It is proved in the Appendix to Chapter 3. 

Theorem 3.6 Let K be a hybrid operator, let C be any LTI operator 

such that Kd+CK) 1 is L.„ -stable, and let R and R 1 be LTI -stable 

2e 2e 

operators. -K 1 is outside cone (C,R) if 


a 

max 



/r g (w) 


for all to 


(3.36) 



(3.37) 


D^(s) 


= D* (s) I + [F C H (s) ] * 


D*(s)] 


-1 


■ (3.38) 


The outside conic sector conditions do not require open loop sta- 
bility of the hybrid compensator. This removes the major restriction 
of the inside conic sector conditions. Taking its place is a require- 
ment that KlI+CK) 1 is L -stable. The composite operator Kd+CK) 1 is a 

ZQ 

hybrid operator (see the proof) . It is stable if all of the poles of 
D^d) have magnitude < 1 (see Subsection 2.3.4). 

3.5.2 Robust Closed Loop Stability 

Sufficient conditions are presented for the hybrid feedback system 
to be closed loop stable for all possible plants G. The robustness re- 
sult follows much more naturally than in Subsection 3.3,3. He .e the 
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hybrid operator -K* is outside of a cone and the plant operator G is 
inside of a cone. This is more natural because the nominal plant G 
and its additive or multiplicative perturbation can be used to place 

•V 

G inside of a cone. 

Theorem 3.7 Consider the hybrid feedback system. Assume R and R^ 
are L^-stable. The hybrid feedback system is closed loop stable for all 
possible 6's if an Rtju) exists such that 


(i) 

K(1+GK )* is L -stable 
2e 


(ii) 

- K 1 is outside cone (G,R) 


(iii) 

0 . [R( jto) ) > £ (u>) for all u). 

nun — a 

■ (3.39) 


The proof is sketched out here. If condition (i) is true then a 
R(jto) can be constructed such that condition (ii) is true (by Theorem 
3.6) . If condition (iii) is true then G is strictly inside of cone 
(G,R) . Hence, the hybrid system is closed loop stable for all possible 
G's (by Lemma 2.6). 

If the uncertainty of the nominal plant G is modelled by a multi- 
plicative perturbation then (3.39) has the alternate form: 

£ m (w) < o min [R(joj) )/a max [G( jw) ] for all go (3.40) 

If £ m (w) * a then the maximum value of a such that (3.40) is true can 
be used to find guarenteed phase and gain margins. 
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3-6 Command Response of Hybrid Feedback Systems 

Conic sectors can be used to analyze the steady state oonmand 
response of a hybrid feedback system. The objective is to find an upper- 
bound for the quality measure "q", which is defined by equation U.95) of 
Subsection 2.4.6. Either the loop transfer operator GK or the closed 
loop operator (I + GK) 1 can be placed inside of a cone. 

The first approach uses Theorem 3.2 or Theorem 3.5 to place the 
loop transfer operator GK inside of a cone.^ The loop is broken where 
the command signal enters the loop, which is just before the prefilter, 
and the plant is grouped with the hold (case 1 of Section 3.4). The cone 
is placed around either GK (if the nominal plant is used) or GK (if the 
perturbed plant is used) . The center and radius are used via (2.98) and 
(2.99) to compute an upperbound for q. 

The second approach to computing am upperbound for q is to place 

the closed loop operator GK( 1 + GIO 1 inside of a cone. The center of the 

I 2 

cone is (Ml + C) , where C is a LTI approximation of GK. The closed 

I 

loop operator GK( I + GfO has the same structure as a hybrid operator. 

The input-output transformation of GK( 2 + Glfl 1 is given by (2.57), which 
is repeated here: 

y(s) - G(s) H (s) D*^(s) [F (s) r(s)]* (3.41) 

where 

D* „(s) = D* (s) [I + G .* DMs)]" 1 (3.42) 

- CL - -Q 

^Theorems 3.2 and 3.5 must be modified as discussed in Section 3.4. 

^This is equivalent to placing (I + GK) 1 inside of a cone with center 
(I + C) 1 . See (2.106) and (2.107). 
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G *(s) « [HGF (s) ] * (3.43) 

*-q 

Theorem 3.2 (appropriately modified) can be used to construct a cone that 
contains GK (T + GK) 1 . The radius of this cone can then be used, via 
(2.105), to compute an upperbound for q. 

Significant improvements can be made to the second approach by 
taking advantage of the fact that the command signals are restricted to 
a set S C L Rather than meet the conic sector inequality 

|| [GK(I + GK ) 1 - C(I + C) I ]r|| < || R.r|| T 

for all rcL* and T e (3.44) 

- 2e + 

it is only necessary that 

||[GK(I + GK ) 1 - C(I + C) I ] r || ||R,r|| 

lim < lim - 

T"* l|r|| T T-*» ||r|| x 

for all r € S 

By virtue of the fact that S C L * it follows that 

2e 

IIMI T II ML 

lim < lim — for all r e S 

t-*° II ^ || T-*» || r || 

~ l •»* T 

r 

A particular set S C L„ is now defined: 
a 4c 

S a = (sinewaves with frequency < “} 


(3.45) 


(3.46) 


(3.47) 
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An important property of signals in this set is that they do not cause 
any aliasing when they are sampled. 

The following theorem can be used to analyze the steady state res- 
ponse to conmands that are members of $ . The proof of Theorem 3.8 is in 
the Appendix to Chapter 3. 

Theorem 3.8 Consider the hybrid feedback system. Assume that 

GK( I + GK) 1 and C ( I + C) 1 are L -stable. Then 

2e 


lim 

t-k» 


tGK(I + GK ) 1 - C(I + C) 1 ]r|] T 



< 


lim 



for all r € S 

a 

if 


(3.48) 


a min [ - (jW)1 - [r io (W) + r il (C0) 1 1/2 for < ? (3.49) 

where 

r io ( “> ■ k J 0 IlSk^Scllll 2 

r ( «> = il ^ - -d - (jW> " ?<I + 2 )-1( ^> II 2 ■ (3-50) 

Theorem 3.8 can be used to find an upperbound for the quality measure 

q. Let the set of command signals be sir.ewaves with frequency less than 

to (which is a subset cf S ) . Then it follows that 
o a 


q < q (1 + r ) 
- o 3 


(3.51) 



where 
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q • max a I (I + C)’*(jw)l (3.52) 

o .. max 
o< 0K(0 
- - o 

r q « max a (R_(jw)l (3.53) 

jo . max — j 
cKaXco 
- - o 

This completes the section on using conic sectors to analyze the 

steady state response of the hybrid feedback system to command signals. 

The suggested approaches are (1) use Theorem 3.2 to construct a cone 

that contains GK, (2) use Theorem 3.2 to construct a cone that contains 

GK( I + GK) 1 , and (3) use Theorem 3.8 to construct a cone that contains 

GK( I + GK) 1 but is only valid for command signals r 6 5 . In the example 

— a 

of Subsection 5.2.11 the third approach is shown to give the least 
conservative (lowest) upperbound for the quality measure q. 
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Appendix to Chapter 3 

This appendix contains proofs for Theorems 3.1, 3.2, 3.4, 3.6, 
and 3.8. The proofs of thu remaining theorems of Chapter 3 are 
sketched out in the text. 

This appendix begins by stating and proving a lemma which is an 
intermediate step used in the proofs cf Theorems 3.1, 3.2, and 3.6. 
The lemma is a frequency domain inequality that is used for signals 
that are input; of hybrid compensators. The proof uses the Cauchy- 
Schwartz inequality (Id, p. 30] and Lebosque Dominated Convergence 
121, p. 44). 

Lemma 3 . A Let lj^(s) be Laplace transform matrices for all inte- 
gers n, and let e(s) be the Laplace transform of o t L . Furthermore 

» *■ 

(to assure Leber.que Dominated Convergence) assume for ju'j suf- 
ficiently large that 


'■ IlK/juO ir 1 — fo1 SOmf * a '^ > 0 ( 3 . A . 1 ) 

n j u' J 


Then 


1 

2n 


II i: ( iu))e _( ju>- !«0 n)|| ‘ da; <•_ 

n 


1 


r"' 

\ ■ 1 1 

— < \ 

• l! 

y n 


't 


i.la' 


(3. A. 2) 


Proof of 


inequalities. 



-CO 
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Lemma 3. A The proof is a long series of equalities and 
Justification for each step is included in parenthesis. 


1 2 K (ju>)e(ju>-ju> n))r do 
n 




(because ||a+b|y < H a.|y + ||b|^> 


da) 


(3.A.3) 


< - 1 - |°° |e ||K n (jw) || *|l e(ja>-ja) g n) || £ J 2 du 


(3.A.4) 


(because | A b J <_ |a|| • ||b|| ) 


'- 2 K j IIJS n ^ w )H 2 J *| Z .,e(ju>-ju) s k) || 2 j da) (3.A.5) 

2 

(By the Cauchy-Schwartz inequality. Let = || ( jto) || 

2 

and b = lle(ju)-jo) n)||_ . For each w, a and b are non- 
n " — s E n n 

negative real numbers. Let a and b be vectors (of 

infinite length) containing the a ' s and b 's. The 

n n 


Cauchy-Schwartz inequality states that 
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|£(jU)-jw s k) II E da) 


(3. A. 6) 


(move sum over k outside integral, which is valid 
by Lebesque Dominated Convergence) . 




K (juH-ju) k) 
— n s 



^ 2 

|£(ju>) || da) 


(3. A. 7) 


(for each k change the variable of integration 

to o)=ai-u) k) . 
s 


2tt 



K (jca+jo) k) 
~n s 



• ||e(3'1)) ||^ da) 


(3. A. 8) 


(move sum over k back inside the integral) . 


This completes the proof. 
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Proof of Theorem 3 . 1 The objective is to show that equation (3.1) 
is an upperbound for the gain of a hybrid operator. The proof essentially 
boils down to an application of the frequency domain inequality of 
Lemma 3. A. Some extra work is required, however, both before and 
after this lemma cam be applied. The operator gain is defined in the 
time domain using truncated norms [see (2.79)]. A truncated function 
is defined so that function norms can be used in place of the trun- 
cated norms. function norms axe needed so that Parseval’s Theorem 
cam be used to switch from the time to the frequency domain. While 
in the frequency domain Lenina 3. A is applied. Parseval's Theorem is 
then used to switch back to the time domain, aind finally the truncated 
function is used to switch back to truncated norms. 

The truncated function is defined by 

j e(t) t £ T 

e^t)^ ) (3.A.9) 

| 0 t > T 

The proof uses the notation: 

K (jio)a H( jw)D* ( jw) F(jw-jw n) (3. A. 10) 

n r ~ “ s 

and the following: 

I E ll^ju-ju k> || 2 -EE IlH^F, J | 2 

k n k n 


(3. A. 11) 
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i £ E [II4II 2 • ll°*H 2 • ll? n+ ) ( l| 2 j (3.A.1J) 

- [2 IlHfclt 2 ] -l|D*|| 2 •[rllr.ll 2 ] (3.A.13) 


2 

The last step, (3.A.13) , is possible because I ||f || is periodic with 

n 

period w s . 

The main part of the proof is now presented. Fo. all e £ L and 

2e 

all i e R +: 


2 2 

j | #Ce 1 1 ^ = j | ^ e T 1 1 T (because K is causal) (3. A. 14) 


(3.A.15) 


(3. A. 16) 


00 

2 

|| (Ke ) ( jw) || da> (Parseral ' s Theorem) 

N l L 

OO 

(3. A. 17) 

00 

||z K^juOe^joj-ju^n) ||^ dw (3. A. 18) 

a, n 

[where K is defined by (3. A. 10)] 




£ ll^ e T ll^ (norm increases as T -+• °°) 


OO 

|| (^e T ) (t) j| ^ dt (by definition) 

—OO 
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2TT 


[ Z Z ||K i (jw-jw s k) |1 2 1 
Ik n ■* 


|£ r (jw)|| do) 


(3. A. 19) 


(by Lemma 3. A) 


<-k \ * It d* || 2 - [e ||F n || 2 J -II II : 


du> 


[by (3. A. 13) ] 


(3. A. 20) 


i r 

2TT J 


^.(jw)!! du) 


(3. A. 21) 


where g = 


max 

0<w<J 

T 


iK i|2 l 


1/2 r ,iV 2 | 

•||D*|| •[l||P n l| 2 j 


(restricted frequency range used because 
the bracketed term is periodic with period 


and is symmetric about w= 0) 


(3. A. 22) 


g 2 IMl 


(Parseval’s theorem) 


(3.A.23) 


= g 2 II s II 2 


[e^(t) = 0 for t > T] 


(3.A.24) 


By definition of the gain of an operator: 


K|j . = sup 


Ke 






(3. A. 25) 


This completes the proof of part (a) . 
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Part (b) of Theore 3.1 is proved by showing that the following 
input (for the SISO case) achieves the upperbound: 

e(t) = £ la I • cos [(a) -U) n)t + Arg (a )] (3. A. 26) 

1 n 1 (os n J 

n 


where a = d*(-jw ) f(-jw + jw n) 
n o o s 

io q = frequency the maximizes (3. A. 22) 


The norm of this function, expressed as a limit, is: 


lim _1 

-f-KD f 




(3. A. 27) 


The output of the hybrid compensator is 


u(t) = ~ l |b k | • cos |(o) o — u> s k) t + Arg(b k )J 


(3. A. 28) 


where b R = h( ja^-juMc) • j l |aj 

n 


and the norm of this output signal is 


lim 

T~MX> 


tWt 


2l 1/2 

| h (ja) o - j u) s k)| J 


\ l l\l 

l n 


(3. A. 29) 
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The gain for this Input signal is 


lim 

x-h» 



<3 


(3.A.30) 


This completes the proof. 

Remark The input (3. A. 26) and the output (3. A. 28) have sinewaves 
at the same frequencies. The a n 's in (3. A. 26) are chosen so that the 
Cauchy-Schwartz inequality is met with equality [see (3.A.5)]. 

Conjecture {Multivariable version of Theorem 3.1, part (b) ) . 

The gain of the multivariable hybrid operator is 


|K| 


max 

TT 

0<W<- 
T 



(3.A.31) 
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Proof of Theorem 3.2 The objective is to show that tha hybrid 
operator K is strictly inside cone (C,R). This proof is similar to the 
proof of Theorem 3.1, and again essentially boils down to an applica- 
tion of Lemma 3. A. Parseval's Theorem is used to pass between the time 
and frequency domain, and the following truncated input is used to 
pass between the truncated norm and the L 2 function norm: 

! (%) (t) t > T 

(3. A. 32) 

0 t < T 

The input is truncated only after it is convolved with the radius. 

Rather than show that || (K-C)e||^ < ||Re|[^. , the main part of the proof 

shows that j|(K-C)R I e || , < j|e || , . Yet another complicating factor 

- T 2 ~ T L 2 

is that the e term must be manipulated to show strict inequality. This 

is possible because R is L_ -stable. The middle part of this proof 

2e 

[steps (3. A. 45) to (3. A. 50) , which includes the application of Lemma 3. A) 
is the only part that differs significantly from [9, Theorem A. 4]. 

It is convenient to introduce the notation 


K 
— n 



-• H(s) DMs) F(s) - C(S) 


n=0 


— H (s) D* (s) F(s-ju) n) 
T — — — s 


nj<0 


(3. A. 33) 


which is used to express the Laplace transform of (K-C)e as 
u(s) - c(s)£(s) = Z (s) e (s- ju) g n) 


n 


(3. A. 34) 
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Midway through the proof the following Inequality is used: 


PAGE is 
quality 


Z Z IlK^jw-jw JO II 2 - Z Z ||k (ju>-jwk)|| 2 + 2 11^(30)- 
k n k n/0 k 

m I 1 3 Ke*VJ 2 

k n T 

-f 3 Kfi’U 2 

k T 

*1 Il?« 1 c5*^-c k ll 2 

k 

1 7 lk 11 ^"1 ■ |l - 112 'In 115 " 112 ] 'V"’ 

= r 1 (ca) - r 2 (w) + r 3 (u) 

= r (oi) 

Part (a) is now proved. For all e e l 2e and all T e R + : 

|| (K— C) e || ^ = II ( IC-C) it ^ (because R 1 exists) 

= |! (K-C)R X e T ||^ (because R is causal) 

< || (K-C)R I e T ||? (norm increases as T _HX> ) 
L 2 


II 2 

(3.A.35) 


(3.A.36) 

+ r 3 (0J) 

(3. A. 37) 
(3. A. 38) 
(3. A. 39) 

(3. A. 40) 
(3.A.41) 
(3. A. 42) 
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CftlGiiVii. r * 

0F p 00 « Quality 


||(K-C) R*e (t) ||p dt 

n 


(3. A. 43) 


(by definition of the i- 2 norm, the integral 

exists because K,C, and R 1 are L _ -stable) 

2 e 


1 f°° 

= 2TT ||(R-C)R I e_ T (ju) |j 2 du (Parseval's) (3. A. 44) 


r°° 

■ 5 

* n 


( jtu) || E dw [by (3. A. 33)] 


(3. A. 45) 


U Z III, II 1 2 

' Ik n 

KJO 


! ] • He" 1 : 


'^.(jw) II E dw 


(by Lemma 3. A) 


(3. A. 46) 


-OC 

1 ~ 2 ~ r(u)) * ||R -1 e T (juj) || E dw [by (3. A. 39)] (3.A.47) 


- 2n r(w) *ll£ 1( 3W)|| 2 * ||e T (jo)) ||| dw 


(3. A. 48) 


f 

j * °~J n [R( 3 W)] • ||e r (jto)|| 2 dw 


(3. A. 49) 


< -L 

- 2tt 


(1-0 * ||e x (jco) II ‘ du 


[by (3.3)] 


(3. A. 50) 
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' U-e) ||e |jf (3.A.51) 

L 2 

- (1-E) ||e T ||J <3. A. 52) 

- (1-e) {] Re || ^ (3.A.53) 

1 IMI t - s ’ll;ll? (3.A.54) 


where e' » e ||R||^ 


It has been shown that K is strictly inside cone (C,R) . This completes 
the proof of part (a) . 

The objective of the proof of part (b) is to show that the optimal 

C(s) of (3.7) minimizes the lower bound of a . £R( jco) ] . Note that 

" mm — 

C(juj) enters only in r^(u }) ; and that r^ (w) , r 2 (w) , and r 3 (w) are all 
nonnegative real numbers. If the optimal center is used men r , (to) -0, 
which minimizes r^tco) - r 2 (<u) + r 3 (w) . This completes the proof of 
Theorem 3.2. 


-135- 


Proof of Theorem 3.4 This theorem gives sufficient conditions for 

robust closed loop stability when the hybrid operator is inside of a 

conic sector. The objective of the proof is to show that conditions 

(ii) and (iii) of Theorem 3.3 are true for all possible G's. 

The first part of this proof uses the LTI robustness result (2.22) 

Conditions (ii) of Theorem 3.3 is that Gd+CG ) 1 is -stable for all 

2e 

possible G's. The assumption is made that Gd+CG ) 1 is L^-stable. It 
follows from (2.22) that Gd+CG ) 1 is L 2e ~stable for all possible G's if 

o . [I + G C(jU))] 

V“> < ~ iGc I iS ) - — £ ° r al1 “ (3 - A - 55) 

max J 

This inequality is implied by (3.15) of Theorem 3.4, Hence, if (3.15) 
is true then condition (ii) of Theorem 3.3 is true for all possible G's 
The second part of the proof uses the following string of in- 
equalities: 


0 (R G (I + C G)" 1 ] < 

max — — 


< 


0 

max 


a 

max 


[R Gl/O . [I + C G] (3. A. 56) 

min — 

IR G] (1 +IJ/0 , [I + C G) 

m min — 


(3.A.57) 


(by the triangle inequality) 
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1 V tB 01 /{«,!„ H + ££! I££l 

(3. A. 58) 

(by property 10 of Table 2.1) 

< 1 (if (3.15) is true) (3.A.59) 

Hence, if (3.15) is true then conditions (iii) of Theorem 3.4 is true 
for all possible G. This completes the proof. 
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Proof of Theorem 3.6 The objective ia to show that -K 1 ia outside 

of cone (C,R). This will be true (by Lemma 2.2) if and only if 

|| RF (I+CK) 1 1| . <_ l. The proof of Theorem 3.6 ia basically an applica- 

L 2 

tion of Lemma 3. A to show that the composite operator has gain <_ 1. 

Condition (i) of Theorem 3.6 assumes that the following feedback 
system is closed iooj. stable: 

e = - Cu \ 

u - Kc > (3. A. 60) 

*= Ru / 

A block diagram of this feedback system 1 is shown in Figure 3.A.I. The 
closed loop operator is RMI+CK) 1 . 

Given the Laplace transform e.^®) **** * n P ut * the Laplace trans- 

form of the output is 

U-j (s) 15 jUs) H(s) D^(s) ~ X F (s- jtu^n) (s- ju^n) (3. A. 61) 

n 

An alternate expression for (3. A. 61) is obtained by defining 

K ( s ) ^ ~ R(s) H ( s) C*o(s) F(s-ju) n) for all n (3. A. 62) 

n t 'Cv s 


^Tho subscript "3" is used to be consistent with Figure 2.10b of 
Subsection 2.4. r 5. 
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and then by substituting (3. A. 62) into (3. A. 61) to obtain 


u_(s) = E K (s)e(s-joJ n) (3. A. 63) 

— J — n — s 

n 

The closed loop system has the same structure as a hybrid compensator, 
except that D*(s) is replaced by its closed loop counterpart D^(s). 
Define the truncated input 


e (t) 
— T 


e(t) t < T 

0 t > T 


(3. A. 64) 


The main part of the proof now follows. For all e e L^ e and all T e R + 


|RK(I+CK) I e|| T = |iRK(I+CK) I e T || T 2 


(3. A. 65) 


< HRKd+CK ) 1 ^!! 2 


(3. A. 66) 


_1_ 

2tt 


I K^( ju))e^ (jw-jw g n) || E do) V 3.A.67) 
n 


■ *ri“ 


|K n (ja)-ju) g k) 


e^tjw) I! I dw 


(3. A. 68) 


(by Lemma 3. A) 
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« 5T f ^[Ml£ll 2 -llH k ll 2 j-ll%ll 2 -[^l|r n l! 2 ]-lls T (3«.) 


,du) 


(R^ = R because R is periodic) 


(3.A.69) 


= | I|r || 2 • r 9 (w) *||e r (ju)t ! 2 dw 


(3. A. 70) 


tr g (o)) defined in (3.37)] 


< — 
- 27t 


r« 


e (j0))|| dii) 
1 E 


(by ( .3*0] 


(3. A. 71) 


.2 
e . 
~t 1 L 


(3. A. 72) 


- «;»? 


(3. A. 73) 


It has been shown that the operator RKd+CK ) 1 has gain < 1, and there- 


/I , 


fore that -K is outside of cone (C f R) . This completes the proof. 


i 
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Proof of Theorem 3,8 The objective is to show that the inequality 
(3.46) is satisfied for all r 6 S . The property of signals r 6 S that 

a y a 

is exploited in this proof is that these signals are not aliased when 
they are sampled. This .Lack of aliasing eliminates the need to shift and 
add the prefilter F(j(o). 

The following is valid as T-*» and for all r € 5 a : 


|| [GK(I + GK ) 1 - C (I + C) 1 ]^* 


= 1|CGK (I + GK) 1 - C(I + C) I ]r||. 

~ 2 

•00 

= ~ 1 1 GHD ! ^ Z F r - C(I + C) _1 r || 2 d £0 

2ir 11 cz T n -n-rt - - - -"e 

J— 00 


(3. A. 74) 


(3. A. 75) 


2tt 


(2ktl)^ 


(2k-l) 


7T 


f 5 ?n- r n ' + 9> _1 rll E 2d “ 


(3. A. 76) 


-±- I 
2ir k 


■n/T 


-7T/T 


G,H,D ! ^ E F r - C, (I + C, ) _1 r, ||l dw 
-k-k-ct T n -n+k-n+k -k - -k -k"E 

(3. A. 77) 


2tt 


>7T/T 
J -TT/T 






2 

E 


+ 



^ Fr - C(I + C) 



(because there is no aliasing) 


d to (3. A. 78) 


< 

- 2t 


i r 7T / T r 

l r 

W k/ 

n/T L 


* 1 


k^o II S&Scl T 5 


* 1 


O / -P . s-1 \ 


-In 2 


n 


ii 2 



I A 
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This completes the proof. 


(3.49)] (3. A. 80) 

(3. A. 81) 
(3 .A. 82) 
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4. NUMERICAL CONSIDERATIONS FOR THE COMPUTATION OF THE CONE RADIUS 
4.1 Introduction 

The previous chapter contains many different cones that the hybrid 
operator is either inside of or outside of. For each of these cones the 
radius must be computed. In this chapter it is shown how to compute 
the different radii. 

Equations for the radius are given in Theorems 3.2, 3.5, and 3.6. 
The equations contain several parts that are combined in straightforward 
ways. What is not straightfoward is that each of the parts contain 
inf inite series that must be summed. The infinite series take the form 



(4.1) 

(4.2) 


The "A" matrix is either the prefilter, the hold, or any of several 
different combinations of the prefilter, hold, and plant. Equation 

(4.2) has appeared earlier as equation (3.11). 

The computational results of this chapter are grouped into three 
sections. In Section 4.2 the infinite sum (4.1) is approximated by summing 
a finite number of terms. Sufficient conditions are given for the 
remainder to be bounded. In Section 4.3 the double infinite sum (4.2) 
is approximated by summing a finite number of terms. Analysis of the 
double infinite sum is aided by splitting it into several single in- 
finite sums. The double infinite sum is finite if the single in- 
finite sums are finite. In Section 4.4 analytical solutions for 
s (oj) of (4.1) are presented. Analytical solutions have been found 
when A(s) = a (s ) is single input single output (SISO) . The major 
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results of this chapter are summarized in three theorems: 

Theorem 4.1: Upperbound for the remainder of the truncated 

infinite series for s (o» . 

Theorem 4.2: Analytical solution for s(oi) when a(s) has a state 

space realization. 

Theorem 4.3: Analytical solution for s(u>) when a (s) = h(s) a(s), 

where h(s) is a zero-order-hold and a(s) has a state 
space realization. 

It is helpful to understand some of the qualitative behavior of the 

infinite series for s(co). The sum iseven[i.e. s (to) = s(-to)] and periodic 

with period to li.-e. s(to) * s(to-kto ) for any interger k] . Because of 
s s 

these two facts it is only necessary to compute s(to) for points between 

0 and 01/2= TT/T (the foldover frequency). If A(jto) is small above the 
s - 

1/2 

foldover frequency then the aliasing is insignificant and s (co) will 
be approximately equal to || A( jto) || for |to| < tt/T. 

An example will help to clarify this qualitative behavior. Consider 
the first order lag 


a(s) 


s + 1 


(4.3) 


A magnitude Bode plot of a (jto) is shown in Figure 4.1. The same Bode 

1/2 

plot contains several different plots of s(0)) , each for a different 

sasqple interval T. 


The break point of a (jto) is at w = 1. If 1 << to then s(to) 


1/2 s 


|a{jU))j for o < to < to /2. On the other hand, if to << 1 then s(to) s 
” ” s s 

1/2 

(T/2) for all to. At high frequencies (to << cu) it is always true 
that |a(jw)| << s(co) 1/2 , because |a(jw) | rolls off and sfoo)^ 2 is 


periodic. 



Magnitude 



Frequency 


Figure 4.1: Comparison of j a ( joj) | and s( 0 )) 1//2 
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4.2 Truncating an Infinite Series 

One way to compute s(co) is to sum a finite number of terms of the 
infinite series (4.1). Whether or not this is feasible depends on the 
behavior of the remainder. 

In this section sufficient conditions are presented for the re- 
mainder to be bounded. One example of these sufficient conditions is 
that the remainder is bounded if A(jco) has an one pole rolloff. These 
sufficient conditions show that the remainder approaches zero as 
the number of terms of the truncated series increases. This information 
can be used to determine where to truncate the infinite series so that 
the remainder is less than a specified amount. 

Some new notation is needed. Let the truncated series and its 
remainder be defined by 



r„(co) A s(co) - s (co) (4.5) 

N = N 


Two assumptions sure made about A(j<o) . It is assumed that above 
some frequency C 0 q that || * (jco) || lies below a straight line asymptote: 

|| A (jco) || < b(co) for co > co q (4.6) 

/co \ P 

where b(co) - f (4.7) 


This assumption is illustrated in Figure 4.2. On a log-log scale the 
straight line asymptote has a slope of -p. If p = 1 then A(ico) is said 
to have at least an one pole rolloff. 



log magnitude 
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Figure 4.2: Assumption that || A(jco) || lies below a straight 

line asymptote 
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The second assumption about A(jco) is that 


* 


|| A(jco) || < 00 for all 0 ) 


(4.8) 


If this assumption is true then the boundedness of r.,(b>) implies the 
boundedness of s (a>) . 

The main result of this section is now presented. The proof is in 
the Appendix to Chapter 4. 

Theorem 4.1 Consider the Laplace transform matrix A(s). 

(a) If (4.6) is true for p > j then 


-=(*)©(*) 


for o < co < cu and 
s 


(b) Furthermore, if (4.8) is true then 


(4.9) 


8 (co) < g for all co and some g < 00 ■ (4.10) 

The upperbound on r , (co) is only valid over the fundamental frequency 

N 

range. Only at the cost of considerable notation can the result be 
extended to other frequency ranges (see the remark in the Appendix to Chapter 
4). The lower bound on N, which is one of the conditions of (4.9), assures 
that all of the terms of the remainder are bounded by b(co) of (4.7). 

When A(jco) has an one- pole rolloff then (4.6) is true for p * 1. 

The upperbound on the remainder is given by 

V“> i 2 (s*) (4 - u > 

The remainder is (approximately) proportional to which is characteristic 
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of an incredibly slow rate of convergence. Each extra decimal point of 

accuracy requires 10 times as many terms in s^(<0) • 

Theorem 4 . 1 can be used to determine the value of N such that the 

2 

upperbound of r N (oo) is less than a specified amount. If this is im- 
portant then it should be recognized that the upperbound (4.9) can be 
tightened, and may result in a value of N larger than necessary. 

The best use of Theorem 4.1 is to help understand the qualitative 
behavior of the remainder (such as when it is bounded and, if bounded, 
what is its rate of convergence) . 


*Slow compared to exponential convergence. 

'< — 2 

For example, if p - 1 then N > (2/e) (co^/u^) guarantees that r^fw) 


for o < co < a. . 

s 
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4.3 Truncating a Double Infinite Series 

The double infinite aeries (4.2) can be decomposed into several 
single infinite series: 


r . (( ) ■ -r ^ E 
4 -7 k n/k 


9*5*? J 


(4.12) 


J H II ?*5* 1 -, il 2 ' IT J H 9*9*?* II 

T * n t K 


(4.13) 




• lie* II 


s II ?„ II 


- \ ? ii 9 k s*? k ir 


T k 


(4.14) 


The inequality in (4.14) can be replaced by an equality if the matrices 
are replaced by scalars. 

The double infinite series (4.12) is finite for each (ti if I |j H |j ^ , 

k 

|| D* || , and E !|f || are all finite for each (d. Each term of (4.12) 

n ‘ " n 

is nonnegative, *.• '4.12) is lower bounded by zero. The subtraction 
in (4.14) cannot result in a negative number. 

A truncated version of (4.12) is defined: 


r 4N l "» ‘-7 * “ Il9*9*? n ll 2 (4 - 15) 

T k“-N n=-N 
n^k 

The remainder r, (oi) - r.„ (o>) can be analyzed by using the remainders 
4 4N 

of the single infinite series of (4.14) . Tiis analysis is cumbersome 
and will not be attempted here. The endpoint of this analysis is 
clear - if H(jo>) and F(ju) both have at least a one pole rolloff then 


1 


Equation (4.15) with N 


20 was used in the examples of Chapter 5. 
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the remainder r 4 (<o) - r 4n (to) is bounded and approaches rero as N + ®. 

Truncating the double infinite series is often prefereable to trun- 
cating the various single infinite series of (4.14). For multi variable 
systems the double infinite series has the distinct advantage of being 
less than or equal to (4.14). 

The double infinite series has another advantage in that it avoids 

a common numerical problem. The double infinite series (which subs to 

X 2 

r (to)] is converted to (4.14) by adding and substracting — E * 

T k 

For some values of to this terxr. may be orders of magnitude larger than 

r. (to). This introduces numerical problems when finite precision 

4 

arithmetic is used, because (4.14) involves the substraction of two 
large numbers, which results in a loss of accuracy. An example of this 
loss of accuracy is shown in Subsection 5.2.6. 


m 
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4.4 Analytic Solutions 

In some cases an analytic solution can be found for s(o>) of (4.1). 

The Laplace transform matrix Ate) must be a scalar, in other words 
A(s) = a(s) . 

Equation (4.19) is used to find analytic solutions for s(o>). This 
equation is now derived. Let b(s) be a Laplace transform, b(t) the 
inverse Laplace transform, and b(z) the z -transform of the samples b(kT) 1 . 
It is well known (4, pp. 77-80] that 

Define 

b(s) = a(-s) a(s) (4.17) 

and then it follows that 

b(jw) = a(-jto) a(jui) = |a(jd))| (4.18) 


s ( 0 )) 


= E 
k 




b, = T*b (z) 


jtoT 

z=e 


(4.19) 


Whether or not (4.19) is useful depends on how easy it is to derive 
an expression for b(z). Given 2 m expression for a(s) the way to find 
b(z) is to respectively find a(-s), bis), b(t), b(nT), and then b(z). 

This procedure will be followed for the following cases: 


1) Zero-order-hold 

2) First order lag 

3) Cascaded zero-order-hold and first order lag 

4) ate) given by a state space representation 


^Any region of convergence (ROC) can be used. If the ROC ^f b(s) is s ^ 

s < Re(s) < s, then the corresponding ROC of b(z) is e”° < |zj < e 1 . 

o 1 
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5) a(s) given by a state space representation cascaded with a 
zero-order-hold . 

The last two cases are sufficiently general to present them as theorems . 

Case 1 If the hold device is a zero-order-hold then a^(s) = h(s) 
has the following Laplace transform:* 


a^s) 



(4.20) 


It is straightforward but tedious to compute b^ (s) = a^ (-s) a^ (s) and 
then to compute the inverse Laplace transform to obtain b^ (t) . A 
shortcut is to recognize that the inverse Laplace transform of a^(-s) 
is a^(-t), and that b^(t) = a^(t) * a^(-t). This convolution is shown 
in Figure 4.3. There is only one nonzero sample of b^tt), which is 
b^o) = T. The z-transform of b (nT) is 

b x (z) = T (4.21) 


Therefore the analytical solution for s(oj) is 


s (to) = Z j a | 2 = T 2 (4.22) 

n In 

which is a remarkable simplification! 

Case 2 If the prefilter is a first order lag then a.^(s) = f(s) 

has the following Laplace transform: 
a 


V s) 


o 

s + a 


(4.23) 


o 

Unlike the previous case there is no shortcut around the page or two 
of mathematical tedium which results in the desired expression 


*The subscripts "1" through "5" are used to distinguish the a's and 
b's of the five cases. 
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Figure 4.3: Zero-order-hold 
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, o az 

b 2 (z) " T ~2 


(4.24) 


- 3z + l 


where 


a = exp (-a T) -exp (a T) 
o o 


3 = exp (-a T) + exp (a T) 


Case 3 Part of the radius calculation has the hold and prefilter 
in the same infinite sum [e.g-r^fa)) of Theorem 3.2]. In this case a 
zero-order-hold is cascaded with a first order lag : 


. -sT a 
, . 1-e o 

a (s) = — 

3 s s + a 


(4.24) 


Again the intermediate steps are left to the reader. The final result 
is 

2 


b (z) = — — - + t 

3 2a o z 2 -8* ♦ 1 


(4.25) 


The a and 3 are the same as for the case 2. The denominators of b 2 (z) 
and b -j(z) are the same. 

Example Consider a SISO hybrid compensator with a first order 
lag and a zero-order-hold. The digital computer z-transform is arbit- 
rary. The results of the first three cases of this section can be used 
to place this hybrid compensator inside of a cone. 

Choose the optimal center of the cone: 


c(s) = - h (s)d* (s) f (s) (4.26) 

From Theorem 3.2 the radius r(s) must satisfy the inequality: 

| r ( jco) j > [r 1 (oj) - r 2 «o)l 1//2 (4.27) 
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After a few substitutions it follows that 


r l (a)> ' r 2 (W> = i d(z) i 2 {b i (z) b 2 (z) * t b 3 (z)i jWT 

z=e 


(4.28) 


Case 4 If a. (s) has a state space realization then an analytical 

— 4 

solution can be found for s (co) . This case is useful when the prefilter 
is more complicated than the first order lag of case 2. 

Let the state space realization of a^(s) be 


x = Ax + bu | 
y = cx j 


(4.29) 


The Laplace transform of (4.29) is 


a. (s) = c (si - A) 
4 - 


-1 


(4.30) 


and the Lap? .*oe transform of a 4 (-s) is 


T T T 

at(-s) * -b (si + A ) 
4 - - - 


-1 


(4.31) 


which has the state space realization 


x = - A x + 


y = -b x 


T \ 
c u I 


(4.32) 


The cascade of the. two systems (4.29) and (4.32) has the Laplace 
T 

transform b (s) = a (-s) a. (s). The state space realization of b. (s) is 


x * AX + 
y = cx 


A 

bu 


(4.33) 


1 T ^ 

The transpace a (-s) is used so that A has the form of a Hamiltonian 
matrix. 
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where 


A = 


T T 
-c c -A 


b = 


c = [o b ] 


The causal impulse response of b^(s) is 


v , . . * At £ 

b (t) = c e b 

4 - 


for t > o 

and the z-transform b 4 (z) of the samples b (nT) is 


b ( z) = E b (nT) z 
4 n=o 4 


-1 


-n 


= h (zl - F) g 


where 


(4.34) 


(4.35) 


(4.36) 

(4.37) 


A 

F = exp (AT) 

A 

g = Fb ► 

A 

h = c 


(4.38) 


The result just derived is presented as a Theorem. 

Theorem 4.2 If a^ts) has the state space realization (4.29) then 
2 

T • b (z) I . „ (4.39) 

4 3 COT 

1 z*=e J 

where b 4 (z) is given by (4.37) and (4.38). ■ 

Case 5 In this case a,, (s) is a zero-order-hold cascaded with a, (s) 
5 4 

of the previous case. This case is a generalization of case 3. It is 
useful when a prefilter that is more complicated then first order lag 
is cascaded witn a zero-order-hold. This case is also useful when 
the plant is cascaded with a zero-order-hold. 


sto) = Z l»J - 
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State space realizations of a 4 (s), a 4 (-s), and fc> 4 (z) are given in 
the previous case. Let h(s) be a zero-order-hold and let 


b^s) = h(-s) h(s) 


h 5 (s) = b 4 (s) b^ (s) 


(4.40) 

(4.41) 


The impulse response b^(t) is a convolution of b 4 (t) [the causal impulse 
response of (4.35)] and b^(t) [the non-causal impulse response in Figure 
4.3], Sample b,.(t) and compute the z-transform of these samples. The 
result is 


b 5 (z) = h 5 (zl - F 5 ) 


-1 


+ k r 


(4.42) 


where 


= exp (AT) 


2 5 


exp (AT) £t dT 


+ F 
-o 


exp (At) b (T-T) dT 


(4.43) 


-5 = - 


k 5 = S 


exp (AT) b(T-T) dT 


One more step is needed to avoid the integrations of (4.43) . The 
trick is to use the augmented matrix. 


C = 


A b 0 
0 0 1 
0 0 0 


I 2n 

t 1 
t 1 


(4.44) 


2n 1 1 
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The dironsions of each of the sub-matrices are indicated, where n is the 
order of a^(s). As shown by Van Loan [22], the matrix exponential has the form: 


*1 2i -! 
exp (CT) — 0 f 2 g 2 


0 0 f. 


(4.45) 


where 


F = exp (AT) 


f 2 = 1 


exp (AT) b dx 


(4.46) 


exp (AT) b(T-T) dT 


Hence, by substituting (4.46) into (4.43) 


T 2l + 


(4.47) 


The result just derived for case 5 is summarized in the following 


theorem . 
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Theorem 4,3 Let a^(s) = h(s) a 4 (s), where h(s) is a zero-order-hold 
and a 4 (s) has the state space realization (4.2S). Then 

sl “> - l Kn |2 ‘ T * b 5 (z ’ ( jo® U - 48) 

z=e 

where b^(z) is given by (4.42) and (4.47). ■ 

This completes the section on analytical solutions of (4.1). The 
Fourier transform matrix A(jw) must be a scalar A( joj) =a ( jw) . The matrix 
exponential is used for the analytical solutions of Theorems 4.2 and 4.3. 
The a(jw) of Theorem 4.2 is represented by a state space description, and 
the a(jw) of Theorem 4.3 is represented by a state space description 
cascaded with a zero-order-hold. 
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Appendix to Chapter 4 

This appendix contains the following proof: 

Proof of Theorem 4 . 1 The objective of part (a) is to show the 
validity of (4.9) , which is an upperbound for the remainder r^fco) of the 
truncated series (4.4). An upperbound is found for each term of the 
infinite series for r N (co) , and then an upper bound is found for the sum 
of the infinite series of upperbounds. 

An intermediate step in the proof is to find an upperbound for 


A 



S /i) 2p 

k-N \k/ 


(4.A.1) 


This series is known to converge^ 
is 


if and only if p > — , 


An upper bound 


N-l 




(4. A. 2) 


This inequality is illustrated in Figure 4.A.I. The sum r„ is the area 

of the rectangular boxes, which is less than the area under the curve. 

The main part of the proof now follows. For o < Co < <o and 

- - s 

N > (&)/&)): 

- os 


- (N+l) 

r (co) = Z || A(jco-jco k) |[ + Z II A(-jco+jco k) || (4. A. 3) 

k=-oo S k=N+1 s 

(by definition) 


An analytical solution is known for the case 
2 


Z (M 2 = ~ 

t =i\k; & 
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= E II A( jOH-jW k) J| 2 + E II A(-jCiH-jO) k) II 2 (4. A. 4) 

k«N+l ' * k-N+1 ~ s 

(In the first series replace k by -k. In 
the second series use the identity 
II A(jS)|| - || A(-jto)||). 


< E [b(a>fw k) 2 + b(-w+to k) 2 ] (4. A. 5) 

k«N+l s s 


(for all k both to+to k and -OH-to k are > co ) . 

s so 



[substitute (4.7) for b(w)] 


<2 E 

k=N+l 


0) 


o 

w (k-1) 

s 


2P 


(This step removes the dependence on gj. 
Hence, the convergence is uniform for 

O < 0) < (0 ) 

s 


(4. A. 6) 


(4. A. 7) 



(4. A. 8) 


(4. A. 9) 


This completes the proof of part (a) . 

Part (b) is true if s„(u)) is finite for all to, which is in turn 

N 

true if || A(to) || is bounded for all to. This completes the proof. 
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Figure 4.A.1: An upperbound for an infinite series 
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Remark The upperbound for r N < w ) is only valid for 10 in the funda- 
mental frequency range c < <•» < W . This is a consequence of the fact 
that s N (o>) and r N («) are not periodic. The finite sum s N (w) is centered 
around the fundamental frequency range. To extend the upperbound to other 
frequency ranges then s^(co) must be defined to be centered around the 
other frequency ranges. 

There is no need to go through this trouble, because s(w) is 
periodic. Compute s^(tt) over the fundamental frequency range and then 
shift it as needed by multiples of co . 

3 
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5, ANALYSIS OF HYB RID SYSTEMS AND NUMERICAL EXAMPLES 
5.1 Introduction 

The examples in this chapter show how conic sectors can be used to 
analyze hybrid feedback systems. In the first example the plant is a 
SISO open loop stable system. An analog lead-lag compensator is designed 
by classical control techniques. The analog compensator is discretized, 
a prefilter and hold are chosen, and then the resulting hybrid system 
is analyzed by conventional z-transform techniques. The first new 
technique is the use of Theorem 3.1 to compute the gain of the hybrid 
operator. Theorems 3.2 to 3,7 are then used to analyze closed loop 
stability and robustness. The analysis techniques based on these 
theorems are grouped according to whether the hybrid operator is 
inside of a cone, the loop transfer operator is inside of a cone, or 
the hybrid operator is outside of a cone. The example continues with 
three more uses of conic sectors - selecting the sample rate, comparison 
of discretization techniques, and robustness with respect to extra 
delay. The example then shows how conic sectors can be used to 
analyze command response. 

The second example is not as extensive as the first. The plant 
is SISO and open loop stable, and the hybrid compensator is a discretized 
versron of an analog integrator. An attempt is made to place the hybrid 
compensator inside of a conic sector, but the attempt fails because the 
radius is not finite. The example then goes on to show that Theorems 
3.6 and 3.7 (for which the -K 1 is outside of a cone) can be used to 
analyze the robustness of this integral control problem. 

In the third example the plant is a linearized model of motion 
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in the pitch axis of a high performance aircraft. Both SISO and multi* 
variable analog compensators are designed using Linear Quadratic 
Guassian (LQG) methodology, and then the analog compensators are dis- 
cretized to form hybrid compensators. For the SISO hybrid feedback 
system the compensator is open loop stable and the plant is unstable. 

The hybrid operator can be placed inside of a conic sector but the 
loop transfer operator cannot be. The example shows how to proceed by 
placing a stable version of the loop transfer operator inside of a conic 
sector. For the multivariable version of this example, Theorems 3.6 
and 3.7 are used to analyze robustness. The margins are shown to be 
conservative . 

Before starting with the examples the difference between the 
nominal and actual feedback systems should be made clear. The noroi ~ 
feedback system differs depending on how the conic sectors are applied. 
The actual feedback system is the same (of course) no matter how the 
conic sectors ara applied. 

If the hybrid operator (or loop transfer operator) is inside of 
a cone then the nominal feedback system is an analog LTI feedback 
system. The center of the cone is part of this nominal feedback system. 
Perturbations of the nominal feedback system have two causes: (1) due 
to the use of a hybrid as opposed to an analog ccmpen jator , and (2) due 
to the actual plant being different from the nominal plant. The 
robustness results (Theorems 3.4 and 3.5) assume that the nominal 
feedback system is closed loop stable (as determined by analog 
techniques) and then give sufficient conditions for all perturba- 
tions in a defined set to preserve closed loop stability. 
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On the ether hand, if the hybrid operator is outside of a cone then 
the nomina'. feedback system is a hybrid feedback system. Perturbations 
are due to the actual plant being different from the nominal plant. 

The robustness result (Theorem 3.7) assumes that the nominal feedback 
system is closed loop suable (as determined by digital techniques) 
and then gives sufficient conditions for all perturbations in a defined 
set to be closed loop stable. 


si 
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5. 2 Lead-lag Compensator - 
5.2.1 Classical Control Desig n 

This example begins as a standard classical control problem . The 
open loop plant has the following transfer function: 


g(s) 


150 

(s+1) (s+3) 


(5.1) 


nd 

The plant is single-input-single-output, 2 order, and open loop stable. 
There is no particular system that g(s) is supposed to represent (this 
is an academic example) , but such a transfer function is typical of 
same kinds of rotating mechanical systems. 

The objective of the classical control design is to find a com- 
pensator k(s) that meets the following specifications: 

(i) steady state error < 2% to a step input (dc gai ' > 50) 

(ii) phase margin > 45° 

(iii) maximize crossover frequency, but not above the uncompensated 
value. 


The dc gain of g(s) is 50, so no extra rain is needed to meet the 
steady state error requirement. Phase lead is needed around crossover 
to meet the phase margin specification. A lead compensator provides the 
phase lead, but cannot be used because it would raise the crossover fre- 
quency and violate the third specification. So phase lag is added betow 
crossover, followed by phase lead around crossover. 

Having decided on the type of compensator (phase lag followed by 
phase lead) there are numerous systematic and unsystematic ways to 
select the parameters. This is fairly easily done using a Bode plot. 

The details don't concern us, only the final result: 

(s+3) 2 

(s+ .4) (s+22.5) 


k (s) 


(5.2) 
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The magnitude and phase Bode plots of g(s) and gk(s) are shown in 
Figure 5.1. The closed loop system is stable, with poles at s = -11.7, 
-3.0, and -6.1 ± 1.2 j. The bandwidth is 7 rad/sec (down from the 
uncompensated bandwidth of 12 rad/sec), and the phase margin is 60°, 

The phase margin specification has been exceeded by 15°, which is some- 
what conservative, but is in anticipation of the unavoidable phase lag 
due to the hybrid implementation. 

5.2.2 Hybrid Implementation 

The analog compensator is converted to a hybrid compensator. The 
different parts of the hybrid compensator are: 

f(s) = — -500 . (poles at s = -35± 35. j) 

s + 70s + 2500 

TT 

T = .031416 sec (foldover = — = 100 rad/sec) 

. .80498 (z-,90993) i 2 

lZ) ~ (z-. 98750) (z-. 47744) 

i " sT 

Ms) = (zero-order-hold) 

s 

The sample rate is chosen so that the foldover frequency is about 14 
times the bandwidth (100 versus 7 rad/sec) . 

The prefilter is a 2 n< ^ order Butterworth filter with a break 
poxnt at 50 rad/sec. It contributes 11° of phase lag at 7 rad/sec 
and has a magnitude of .24 at 100 rad/sec. 

The z-transform d(z) is a discretized version of k(s). The Tustin 
with prewarping method was used (see Subsection 2.3.1). The prewarped 
frequency is chosen to be ^ = 3 rad/sec, which is the natural frequency 

of the zeros at s = -3. 



magnitude 


- 170 - 


ORIGINAL PAGE 18 
OF POOR QUALITY 








-i71- 


A zero-order-hold is used to convert the digital sequence at the 
output of the computer to an analog signal. The zero-order-'hold is a 
standard choice for control system design, and in this example it 
contributes 6° of phase lag at 7 rad/sec. 

This completes the choice of the hybrid compensator. There are 
many different ways to choose the hybrid compensator, but it is not the 
intent of this thesis to recommend one way over another. Rather, the 
intent is tc provide tools to analyze a particular choice of a hybrid 
compensator. 

5.2.3 Digital Analysis 

Before going on to the conic sector analysis techniques we do a 
quick digital analysis. The hold and prefilter are grouped with the 
plant and the combination is discretized^: 

Z {h(s)g (s)f (s) } = g (z) 

a 

.0094 (Z + 6.059) (z + .609) (z + .066) ... 

2 T (5 ‘ 4) 

(z - .969) (z - .910) [(z - .145) + (.3) ] 

This combined with the computer z-transform [d(z) of (5.3)] gives the 
discrete loop transfe- function: 

t 3 (z) = g d (z) d (z) (5.5) 

iOJT 

Its discrete Nyquist plot (t^(z) evaluated for z = e ) is shown in 
Figure 5.2. The number of clockwise encirclements (zero) of the 1 
point equals the number of open loop unstable poles, so the digital 
closed loop system is stable. Closed loop stability can also be 

^"The script z indicates the z-transform of the samples of the inverse 
Laplace transform, see [4, Sections 3.4 and 6.2], 
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determined from the location of the digital closed loop poles , which all 
have magnitude < 1: z - .910, .873, ,812± ,204j, and .109± ,271j. 

From the discrete Nyquict plot of Figure 5.2 it is seen that the 
digital closed loop system has a phase margin of 46° and a gain margin 
of lOdB. This implies that the hybrid closed loop system has these 
same phase and gain margins (see Subsection 2.3.4). Therefore, the hybrid 
closed loop system meets the phase margin specification of > 45®. 

5.2.4 Gain of the Hybrid Operator 

The first result of Chapter 3 (Theorem 3.1) is an upper bound on 
the gain of the hybrid operator. For SISO hybrid operators (as in 
this example) the upper bound actually is the gain. 

The gain of the analog compensator is 


max |k(jw)| = |k(jO)| = 1 (5.6) 

CO 


and the gain of the hybrid operator is 


1/2 1/2 


li K ||, = max — 

2 r T 

£ i\l 2 

k K 

* |d*| • 

Z |f | 2 
n 1 n 1 

°<w<T 

- _ 


- 

0=0 

= 1 

= 1 

= 1.004 


The maximum for the hybrid operator occurs at co = o. The hold term 

is =1 at co = o (actually for all k)) , and the computer term is =1 

at co = o (z = 1) . The extra gain is due to aliasing of the prefilter, 

which, in this example, does not amount to much. 

Different prefilters that cutoff less sharply or cutoff at higher 

frequencies will result in a hybrid compensator with a larger gain. 

For example, if f(s) - 50/(s+50) then |j K ||. = 1.094, and if f(s) = 

L 2 

100/ (s+100) then 1 1 K ]j^ = 1.309. As the gain gets higher the 
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compensator is more sensitive to noise. Thus, the higher gains are an 
indication of inadequate prefiltering. 

5.2.5 Hybrid Operator Inside of Cone 

Theorem 3.2 can be used to construct a conic sector that contains 
the hybrid operator. The choice of renter is arbitrary. Here we 
choose the center that minimizes the radius: 

c ( s) = | h(s) d(e ST ) f (s) (5.8) 

The loop transfer function of the nominal feedback system is cg(s). 

At low frequencies c(s) should be a good approximation of k(s) . 

This is checked with the magnitude and phase Bode plots of Figure 5.3. 
The magnitudes are close (within 1 dB) below 30 rad/sec, at which 
point the extra rolloff due to the prefilter and hold causes c(s) to 
diverge from k(s). The phase difference shows up at lower frequencies ~ 
the phases are close (within 5°) below 2 rad/sec. At 7 rad/sec the 
phase lead is reduced from 28° to 12° (which was anticipated) . Above 
100 rad/sec the phase of c(s) rapidly swings over the entire 360° range. 
The equation used to compute the radius is 

1/2 

for o < 0) < (5.9) 

which is a truncated version of (3.11). Because the optimal center is 
used, the r 3 (w) term of (3.6) is equal to zero. A plot of the radius 
is shown in Figure 5.2a. The radius is periodic with period a) g = 700 
rad/sec , 

More important than the magnitude of the radius is the relative 
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magnitude of the center and radius. In Figure 5.4b the magnitude of 
c *r is plotted. This is called the "multiplicative radius." When it 
is «1 then the center can be considered a good approximation of the 
hybrid compensator,* 

When the multiplicative radius is near unity or greater then tne 
center is a poor approximation of hybrid compensator. Over this frequency 
range the analog loop transfer function eg (jw) should be rolled off (have 
magnitude < 1) in order to meet the stability and robustness requirements 
of Theorems 3.3 to 3.5. In this example | eg ( jeo) { < .1 over the fre- 
quency range where |c *r(ju>)| > 1, which should be considered "good". 
Stability 

Theorem 3.3 gives sufficient conditions for the hybrid feedback 
system to be closed loop stable. The three conditions of Theorem 3.3 
will now be checked. 

Condition (i) is that the hybrid operator K is strictly inside 
cone (C,R). The center and radius were constructed so that this is true. 
Condition (ii) is that the nominal system (with the loop transfer function 
eg) is closed loop stable. This is verified by the Nyquist diagram of 
Figure 5.5a.^ The third and last condition is that 

1 rg (1 + eg) * (jw) i < 1 for all w (5.15) 


For example, the center is a good approximation when |c *r(jw)| < .1, 
which is true in Figure 5.4fc for w < 6 rad/sec. 


The number of open loop unstable j*>les is zero. The nu.,,jer of clock- 
wise encirclements of the -1 point is zero. Since these two numbers 
are equal, the analog system is closed loop stable. 


imaginary 


7ft. 
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This is verified in Figure 5.3b. Hence, the hybrid system is closed 
loop stable. Of course, we already know this from the earlier digital 
analysis . 

Robustness - First Attempt 

Beyond mere stability is robustness. Theorem 3.4 gives sufficient 
conditions for an entire set of plants to be closed loop stable. In 
this example the plant uncertainty is characterized as a phase uncertainty 
at the crossover frequency as high as 45°. The first step in applying 
Theorem 3.4 is to model this uncertainty as a multiplicative perturbation. 
As explained in Subsection 2.2.7, a multiplicative perturbation of 

t m (u) = a = .77 (5.16) 

covers phase uncertainties up to ±45° (at all frequencies). 

The robustness condition of Theorem 3.4 is graphically checked in 
Figure 5.6. It is seen that 

> - 5? f ° r au ° < 5 - 171 

This corresponds to a phase uncertainty at all frequencies oi 33° , so 
based on Theorem 3.4 we cannot guarantee that the phase margin specifi- 
cation o £ 45° is actually satisfied. Because Theorem 3.4 giver only 
sufficient conditions, we know that the phase margin is > 33°, but we 
do not have enough information to say whether or not the actual phase 
margin is 45° . 

Robustness - Second Attempt 

Theorem 3.5 can be used to analyze robustness. The procedure 
is to (1) choose a multiplicative perturbation, (2) use part -a) of 


magnitude 
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Theorem 3.5 to conscruct a cone around (I + E ) K, and (3) use part (b) 

in 

of Theorem 3.5 to check if -G* is outside of this cone. If part (b) of 
Theorem 3.5 checks out then the closed loop system is stable for all 

plants defined by G = G (I + E m ) . 

In this example the multiplicative perturbation has already been 

chosen (a = .77). The center of the cone around (I + E )K is given by 

m 

5.8) and is the same as the center of the cone around K. The radius of 

the cone around (I + E )K is computed via (3.23) and is plotted in Figure 

m 

5.7a. It is compared with the radius of the cone around K (Cl = o) . 

The nominal plant will be outside of the cone that contains (I + E m )K 

i f , 

jrg (1 + eg) (joj) j < 1 for all W (5.18) 

This inequality is checked in Figure 5.7b, and it is not satisfied. If 
a is backed down to a = .72 then the inequality iis satisfied. Therefore, 
we can say the guaranteed phase margin is 42°, which doesn't quite make 
it to 45°, but is less conservative than the 33° of Theorem 3.4. 

5.2.6 Loop Transfer Operator Inside of Cone 

The entire loop transfer operator can be placed irside of a cone. 

As discussed in Section 3.3, there are three ways that the plant can be 
included with the compensator: 

Case 1: with the hold (replace h with gh) 

Case 2: with the prefi] „er (replace f with fg) 

Case 3: w^th a combination of the hold and prefilter. (replace 
both h and f with (fgh)^/^) 

For each of the three cases Theorem 3.2 will be used to construct a 
cone. After this is done, Theorem 3.3 will be used to see which of the 
three cases meet the sufficient conditions for closed loop stability. 
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A robustness analysis using Theorems 3,4 and 3.5 is then performed for 
case 3 . 

For each case the center of the cone is 

c(s) = g(s) ^ h<s) d*(s) f(s) (5.19) 

The loop transfer function for the nominal feedback system is c(s). In 
Figure 5.8 the magnitude and phase Bode plots of c(s) and gk(s) are 
compared. The center is a good approximation over the same frequency 
ranges as in Figure 5.3. 

The radii for each of the three cases are shown Figure 5.9a. As 
expected, the radius for case 3 is the smallest for each to. At low 
frequencies the radii are significantly different, with the radius for 
case 1 (plant with hold) being the largest. 

The significant differences came as a surprise, and further analysis 
was conducted to explain the differences. The detailed calculation for 
the radii at to = .01 (the leftmost points in Figure 5.9a) are shown in 
Table 5.1. The infinite series that are part of the equation for the 
radius were truncated at ±100 terms. A large number of significant 
digits are needed in Table 5.1 because accuracy is lost when two large 
numbers are substracted (for case 3, r = r^ - r 2 , where r ~ .001 and both 
r^ and r^ are ~ 2500). ± 

An indication that the radius for case 1 (plant with hold) is 
larger than the radius for case 2 (plant with prefilter) is given by the 


This problem with loss of accuracy does not occur when the double 
infinite sum, r^ of (3.11), is used to compute the radius The double 
infinite sum does not contain the large term (r,, ~ 2500) wl ch is 
added and then subtracted fo convert the double^sum to several single 


sums. 
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frequency (rad/*ec) 

Figure 5.9a: Radii of cones that contain 3 different 

loop transfer operators. 



frequency (rad/sec) 


Figure 5.9b: Center and multiplicative radii. 
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Case 1 


r (j .01) = |d 


• (h i ^ 
( 


E f 


= .99970 I [2499. 72229) [1.00844497] - 


* 4.29 


\l/2 

[2499.722234] J 1/2 


Case 2 


r (j .01) = d « / 


’ (l? I 


l g f 

n n n 1 


- r. 


,) 1/2 


= .99970 ([1] [2499.722254] - [2499.722234]) 1//2 


= .00454 


Case 3 


r (j .01) = d 


= .99970 


= .00156 


• ([? l K 

([49. 


9 k f k ' 


2 - * 2 ) 1/2 


997222284579] 2 - [2499.722233745]\ 1//2 


where 


1 r , , 1 2 
r 2 = ~2 - E h - g v f i 


k k k' 


Table 5.1; Calculations for r ( jco) at (u - .01 
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comparison; 

1 If j 2 - 1.008 > 1 - 4r £ | hj 2 for W = .01 (5.20) 

n 1 n 1 T 2 k K' 

The apparently small difference of .008 is amplified by a multiplication 
of 2500, which is approximately the difference (.008 x 2500 = 4) between 
the two radii at to = .01. 

The intuitive reason for including the plant with the compensator 

is that the extra rolloff of the plant will reduce aliasing and thereby 

make the radius smaller. The comparison (5.20) indicates that including 

the extra rolloff with the prefilter is better than including it with 

hold. The theoretical analysis of Section 3.4 indicates that it is 

1/2 

even better to replace the prefilter and hold with (hgf) 

Moving on now to Figure 5.9b, the multiplicative radii (c 1 r) are 
compared. Despite large differences at low frequencies, they cross 
unity at about the same frequency (oj = 80 rad/sec) . All three are < .1 
over the bandwidth of the analog loop transfer function (to < 7 rad/sec) . 

Theorem 3.3 is used to check sufficient conditions for the hybrid 
system to be closed loop stable. The Nyquist plot of c(jto) is shown in 
Figure 5.10a. The Nyquist plot indicates (for each of the three cases) 
that the nominal analog system is closed locp stable. 

It is in the application of the stability test 

| r (1+c) * ( jto) | < 1 for all to (5.21) 

that the most significant differences in the three cases show up. As 
shown in Figure 5.10b, cases 2 and 3 meet the sufficient conditions 
but case 1 does not. Only one of the cases has to meet the sufficient 
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conditions in order to guarantee closed loop stability. It comes as 
some relief to know that case 3 will always be the least conservative. 

An analysis of robustness is performed just for case 3. Figure 
5.11 shows that 

> .74 for all to (5.22) 

This is the single condition of Theorem 3.4 (appropriately modified for 

case 3) . The bound on the multiplicative perturbation is t (to) = a = .74, 

in 

which corresponds to a guaranteed phase margin of 43®. 

Theorem 3.5 can also be used to compute a guaranteed phase margin, 

A cone is constructed around (I + E^JT, The largest a such that 

| r (1+c) _1 | < 1 for all u (5.23) 

is found to be a = .75, which corresponds to a guaranteed phase margin 
of 44°. This is very close to the specification of 45°. 

5.2.7 Hybri d Operator Outside of Cone 

According to Theorem 3.6, a cone can be constructed such that -K 1 
is outside of the cone if the hybrid operator K 1 + GK ) I is /-2 e _sta ^^ e • 
From the digital analysis of Subsection 5.2.3 we know that this stability 
requirement is satisfied. 

The c.-nter of the cone is the nominal plant g(s). The radius is 
computed via (3.36) of Theorem 3.6. Rather than show a plot of the 
radius, a plot of the multiplicative radius |g ^rfjto) | is shown in 
Figure 5.12. The multiplicative radius is moie useful for graphical 
tests of robustness margins. In this example the plot is almost identi- 
cal to the plot of Figure 5,11. To two significant digits the minima 


1+c - 


c| + |r 


are the same (ot 


.74) 
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Prom Theorem 3.7 we know that the minimum of |g ^r(jo)) | can be used 
to guarantee that the hybrid closed loop system is robust with respect to 
a multiplicative perturbation of the form £^(0)) * a * .74. This corres- 
ponds to a guaranteed phase margin of 43*. 

In this lead-lag example there have been 5 separate tests for ro- 
bustness. They are compared in Table 5.2. The actual phase margin (from 
the digital analysis) is 46 e . The guaranteed phase margins of tests 2 
through 5 are not significantly different. Only test 1 is conservative. 


Table 5.2 

Compar i son oif Guaranteed Phase Margins 


Guaranteed 


Test 

1 

Compensator inside --one 

Theorem 

3.4 

a 

,57 

phase margin 
33° 

2 

11 •« 

3.5 

.72 

42° 

3 

Loop operator inside cone 

3.4 

.74 

43° 

4 

i* «• 

3.5 

.75 

44° 

5 

Compensator outside cone 

3.7 

.74 


The 

major advantage of the outside conic 

sector 

analysis is that 


it can be used with open loop unstable compensators and plants. This 
advantage is not needed in this example. 

5.2.8 Selecting the Sample Rate 

The loop transfer operator can be placed inside of a cone (see 
Subsection 5.2.6). The idea now explored is that the multiplicative 
radius of this cone car. be used to systematically select the sample rate. 

The center is a good approximation of the loop transfer operator 
over the frequency range where the multiplicative radius has magnitude 
« 1. By varying the sample rate, this frequency range can be adjusted 
to correspond to the bandwidth of the analog system (i.e., adjust T so 
that [c ^r ( joj) | <<1 for r such that |gk(jw)[ > 1). 
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This procedure *or selecting the sample rate is demonstrated for 
the lead-lag exampl j , a'.I of the components of the hybrid compensator 
must change as the sample rate is varied. The components are selected 
so that the sample rate is the only variable i 

Prefilter; 2 nd order Butterworth with break frequency ■ ~ 

Computer; Tustin prewarped version of k(s) , with 05, * 3 
Hold; Zero-order-hold. 

For each of five sample rates a cone is constructed around the 
loop transfer operator. The multiplicative radii are shown in Figure 
5.13a, and compared with the center of the cone when T = .031416. 

All but one of the multiplicative radii are < .1 over the bandwidth 
of the analog system (w < 7 rad/sec) . The conclusion based on Figure 
5.13a is that all but the lowest sample rate *T = .31416) is acceptable. 
Further tests should be conducted, however, because the fact that 
|c ^r| << 1 over the analog bandwidth is not by itself a guarantee of 
closed loop stability or adequate robustness margins. 

The magnitude oi the multiplicative radius can be compared with the 
magnitude of other multiplicative uncertainties. This comparison can 
also be used to help select the sample rate. The idea is to select the 
sample rate so that |c ^r(jGJ) | < |^(co) | over the bandwidth of tne ana- 
log system. Lowering the sample rate (increasing T) will cause the 
errors due to sampling to dominate the other uncertainties, which is 
undesirable. On the other hand, increasing the sample rate (lowering T) 

1 The five sample rates are (T = .31416, .062832, .031416, .062832, and 
. 0031416) , which correspond to foldover frequencies of (it/t = 10, 50, 

100, 500, and 1000) Tlie case 3 loop transfer operator is used 
[prefilter and hold replaced by (fgh) 1 ' 2 ]. 



magnitude 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 5.13a: Multiplicative radius of cone that 

loop transfer operator is inside of, 
for each of 5 sample rates. 



frequency (rad/sec) 

Figure 5.13b: Multiplicative radius compared with 

multiplicative perturbation due to 
extra delay 
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will not result in any benefit, because the other errors (t^) are not 

affected and will dominate the errors due to sampling . 

For example, suppose the following delay exists in the feedback loop: 
*«* 

g(s) = g£(s) (5.24) 

where £(s) = e T, = .05 seconds 

d 

The delay can be modelled as the multiplicative perturbation: 

g(s) = g[l + £ (s) ] (5.25) 

m 

where £ (s) = e S ‘ ci **1, T, = .05 seconds 
m a 

In Figure 5.13b the magnitude of £^ (co) is compared with the multiplicative 

radius when T = .031416. The magnitude of £ (co) dominates for <u < 100 

m 

rad/sec. The conclusion is that there is nothing to gain by increasing 
the sample rate (lowering T) . 

5.2.9 Comparison of Discretization Techniques 

Another use of the multiplicative radius is to compare discretiza- 
tion techniques. Tustin with prewarping was used in this lead-lag 
example to transform k(s) into d(z), but this is only one of many 
techniques. Tl»e following discretization techniques are compared: 

Tustin with prewarping: (see Subsection 2.3.1) 

Forward rectangle rule: d(z) = k(s) 

s=(z-l)/T 

Backward rectangle rule: d(z) = k(s) 

s=(z-l)/(zT) 

Pole-zero mapping [4, p. 61] : 

P • T 

If p. is a pole of k(s) then e 1 is a pole of d(z) 

1 z^T 

If is a zero of k(s) then e is a zero of d(z) 

If oo is a zero of k(s) then -1 is a zero of d(z) 

Choose constant so that k(s) 1 = d(z) 

|s=o Z=1 
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The standard by which they are compared is the magnitude of the 
multiplicative radius. Each of the loop transfer operators (case 3) are 
placed in cones. In order to make a comparison the center of the cone 
is chosen to be the same for each of the discretization techniques: 

c(s) = ^ h (s) k (s) f (s) (5.26) 

Note that c(s) is not the optimal center. The center and the different 
multiplicative radii are shown in Figure 5.14. 

The multiplicative radius is smallest when the Tustin with prewarping 
technique is used, following by the pole-zero, backward rectangle, and 
forward rectangle techniques. All of the techniques, however, should 
be considered "good enough" for this example, because all of the multi- 
plicative radii have magnitudes < .1 over the bandwidth of the system 
(u) < 7 rad/sec) . 

There is a great deal of freedom in how to discretize k(s). This 
example demonstrates that to a large extent it does not matter what 
discretization technique is used. 

5.2.10 Extra Delay 

Suppose an extra delay exists in the feedback loop. The delay can 
be modelled as a multiplicative perturbation of the plant, and then any 
of the robustness results (Theorems 3.4, 3.5, and 3.7) can be used to 
determine if the closed loop system is robust with respect to this 
multiplicative perturbation. 

The delay may be a computational delay, in which case it will be a 
fraction of a sample period. On the other hand, the delay may be 
characteristic of the plant, such as a transport delay for mechanical 
systems. In this example the delay is chosen to be T. = 


.05 seconds. 


magnitude 


-196- 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 5.14: 


Comparison of discretization techniques. 
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It is between one and two sample periods, to emphasize that there is no 

need for the delay to be an integer number of sample periods. 

Theorems 3.6 and 3.7 will be used to analyze robustness. The basic 

idea is to construct a cone such that the hybrid operator is outside of 

the cone, and then to show that all perturbations of the plant in a defined 

set are inside of this cone. The delay is modelled as a multiplicative 

perturbation, as shown in equations (5.24) and (5.25). The robustness 

analysis consists of checking the three conditions of Theorem 3.7. 

Condition (i) is that Kd+GK) 1 is L_ -stable. This is true because 

2e 

the discrete closed loop system is stable, as shown in Subsection 5.2.3. 

Condition (\i) is that -K 1 is outside cone (G,R) . The existence of 
such a cone is guaranteed by the fact that condition (i) is true. The 
radius is computed as shown in Theorem 3.6. The multiplicative version 
of this radius, jg 1 r(jO))| , is plotted in Figure 5.15. 

Condition (iii) is that 

£^( 0 )) < |c 1 r(jto) | for all oo (5.27) 

This is shown to be true in Figure 5.15. Hence, all perturbations of 

the plant in a defined set are inside of cone (G,R). 

All three conditions of Theorem 3.7 are satisfied. Hence, the 

hybrid system is robust with respect to an extra delay of T. = .05 

d 


seconds . 



magnitude 
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Figure 5.15: Test for robustness with respect to delay of 

Tj = .05 seconds. 
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5.2.11 Analysis of Command Response 

Conic sectors can be used to analyze the steady state response to 
commands, as discussed in Subsection 2.4.6 and Section 5.6. The res- 
ponse to commands is good if th^ following error signal is small: 

e = (I + GK) x r (5.28) 


Conic sectors are used to find an upperbound for q, the quality measure 
defined by 


lim < q 

II r|| T 


for all r 6 S 


(5.29) 


This inequality only has to be satisfied for input signals in the set 
S C , a r»d by letting the truncation time tend to infinity the tran- 
sient errors in the command response become insignificant. 

Three attempts are made to find an upperbound for qt 

(1) Theorem 3.2 is used to construct a cone that contains GK, 
and then (2.98) to (2.100) are used to compute the 
upperbound for q. 

(2) Theorem 3.2 is used to construct a cone that contains 

( I + GK) 1 , and then (2.102) and (2.105) are used to compute 
the upperbound for q. 

(3) Theorem 3.8 and equations (3.51) to (3.53) are used to compute 
the upperbound for q. 

The first attempt fails due to restrictions on its use. The second 
attempt is conservative, and the third attempt works well. 

The first step for all three attempts is to define a nominal analog 
feedback system and then to compute the nominal quality measure q Q . In 
all of the attempts let the nominal loop transfer function be 

c(s) = ^ g(s) h (s) d* (s) f(s) 


15.30) 
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and let the set of ccamaxvl signals be 

S ■ (sinewaves with frequency < .1 rad/sec} (5.31) 

The nominal quality measure is 

q - . Id+c)* 1 (jtt)| - .02 (5.32) 

The magnitude Bode plot of (1+c) * is shown in Figure 5.16a. 

Earlier in this example* Theorem 3.2 was used to construct a cone 
that contains the loop transfer operator GK. The center and radius of 
this cone can be used via (2.98) to (2.100) to find an upperbound for the 
quality measure q. One of the restrictions of this first attempt is 
that the sufficient conditions for closed loop stability given by Theorem 
3.3 must be satisfied. As shown in Figure 5.10b the sufficient condi- 
tions are not satisfied, and therefore the first attempt at finding an 
upperbound for q fails. 

The second attempt uses Theorem 3.2 to construct a cone that con- 
tains the closed loop operator (I +GK) 1 . The center is (1+c) *, 
where c(s) is given by (5.30). Both the center and radius are shown 
in Figure 5.16a. Equations (2.102) and (2.105) are now used to upper- 
bound the quality measure: 

r 2 q o ' o^l |r<j “ )1 * - 09 (5 ' 33) 

q < q (1+r.) * .11 (5.34) 

- O 2 

This upperbound is significantly more than the nominal quality measure 

of q ■ .02. 
o 


*Case 1 of Subsection 5.2.6. The center and radius are shown in Figures 5.8 
and 5.9. 





V 
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fraqutncy (rad/ mc ) 

Figure 5.16a: Center and radius of a cone that contains (I+GK) , 

when the radius is computed by Theorem 3.2. 


I 

i 



Figure 5.16b: Radius of cone that contains (I+-GK) 1 , when the 

radius is computed by Theorem 3.8 . This cone is 

only valid for input signals that are sinewaves with 
frequency < tt/T = 100 rad/sec. 


L. 
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The second attempt at finding an upperbound for q is conservative 
because the radius computed by Theorem 3.2 does not take into account 
(1) the restricted set of input signals and (2) the truncation time of 
(5.29) which tends to infinity. Addressing the second point, the radius 
computed by Theorem 3.2 is valid for any truncation time [see (2.81)], 
and therefore the upperbound for q which uses (2.105) is likewise valid 
for any truncation time. One interpretation of this conservative 
upperbound is that it is conservative for steady state errors because 
it must also be valid for transient errors. 

The third attempt at finding an upperbound for q makes use of 
Theorem 3.8, which explicitly takes into account the restricted set of 
input signals and the truncation time that tends to infinity. The 
radius computed by (3.49) of Theorem 3.8 is shown in Figure 5.16b, The 
upperbound for q is found by equations (3.51) to (3.53): 

r,q * max Jr(j&>)| ■ 5 x 10 8 (5,35) 

3 ° o<cj< -1 

q < q (1 + r_) - .02 (5.36) 

o J 

The increase from the nominal value of q Q * .02 is insignificant 1 

The use of conic sectors to analyze command response is new with 
[16] and this thesis, and we must be careful not to make hasty genera- 
lizations about the very tight upperbound for the quality measure just 
demonstrated. Nevertheless the following statement appears to be 
justified - the nominal analog feedback system gives a very good 
approximation of the quality measure of a hybrid feedback system. 

Three attempts were made to analyze command response. A straight- 
forward use of the same conic sector used to analyze stability and 
robustness (Theorem 3.2) was shown to give a conservative upperbound 
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for the quality measure. Refinements of this conic sector (Theorem 3.8) 
significantly decrease the upperbound, and indicate that the nominal 
quality measure is a good approximation of the actual quality measure. 

This analysis of command response completes this extensive lead-lag 
compensator example. 


cv * 

5.3 Integral Control Example 

The next example demonstrates that the outside conic sector results 
can be used to analyze control systems with integrators in the forward 
loop. These types of control systems are used to achieve zero steady 
state error to constant inputs. 

In thin simple example the SISO plant is 


9 (.) - ^ 

and the compensator is 


(5.37) 


k(s) - -■ 

8 

The analog system has closed loop poles at s = .5 ± j.87. The bandwidth 
is w * .8 rad/sec, and the phase margin is 52®. 

The following hybrid implementation of k(s) is used: 


f(s) » ~ 2 — (2nd order Butterworth, break frequency * 5 \ 

s +7s+25 rad/sec) J 

T « .31416 (foldover frequency * ^ * 10 rad/sec) \ 

T > (5.39) 

d(z) » .1579 ~~ (Tustin prewarped about w = .8 rc^/tec) I 

Z" X o I 

. -sT ] 

h(s) * { zero-order-hold) ' 

8 

At the crossover frequency (0) ■ .8 rad/sec) the prefilter contributes 
13° of phase lag and the hold contributes 7® of phase lag, so we expect 
the phase margin of the discrete closed loop system (i.e. the actual 
phase margin of the hybrid closed loop system) to drop from 52® to about 
32®. 

Due to the digital integrator there does not exist a conic sector 


with a finite radius that contains the loop transfer operctor. The 
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calculations for the radius are performed to demonstrate this fact. 

The optimal center of the conic sector is chosen: 

c(s) = ~<j(s) h ( s) d*(s) f (s) (5.40) 

The radius is plotted in Figure 5.17a, The problem with the radius is 
immediately apparent - it has peaks of infinite height at frequencies 
that are multiples of 2 tt/t. 

A consequence of the infinite peaks is that the sufficient conditions 

for closed loop stability given by Theorem 3.3 cannot be met. This does 

not necessarily mean that t.h< closed loop system is unstable, it may 

simply mean that the sufficient conditions are "infinitely" conservative. 

Condition (iii) of Theorem 3.3 (for the loop transfer operator inside of 

the cone) is that |r(l+c) ^(jw)| < 1 for all 0 ). Due to the infinite 

peaks of r(jco), this condition cannot be met. 

The optimal center and the multiplicative radius are plotted in 

Figure 5.17b. The multiplicative radius, c *r, is independent of the 

-1 

computer z-transform d(z). Hence, c r does not have the problem with 
the infinite peaks. From Figure 5.17b we see that ic ^r| << ]c| over 
the bandwidth of the analog system (o> < .8 rad/sec) . 

There are no open loop stability restrictions for the outside conic 
sector results of Theorems 3.6 and 3.7. The discrete closed loop system 
is stable, as demonstrated by the discrete Nyquist plot of Figure 5.18a. 
Therefore £ cone can be constructed such that is outside of the cor.e. The 

center is q(juj), and the multiplicative radius is jg *r (j'i)) | , where r ( jfo) is 
computed via equation (3.36). The multiplicative radius, g *r, is shown 


Case 3 is used, which means that the prefilter and hold are each re- 
placed by (fgh)!/^. 
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in Figure 5.18b. Its minimum value is ci ~ .53, which corresponds to a 
guaranteed phase margin of 31°, and which is only 1* less than the 
actual phase margin of 32*. 

In this short example it has been demonstrated for integral control 
problems that (1) there does not exist a conic sector with a finite radius 
that contains the loop transfer operator, (2) Theorem 3.3 cannot be used 
to determine closed loop stability, and (3) Theorems 3.6 and 3.7 can 
be used to determine guaranteed robustness margins. 





-208 


ORIGINAL PAGE IS 
OF POOR QUALITY 







-209- 


5.4 High Performance Aircraft 

The next example contains both a SISO and a 2-input 2-output control 
system design for the pitch axis of a highly maneuverable aircraft. This 
aircraft has served as a test bed for multivariable design techniques 
[23, 24, 25, 26] . The analog compensators used in this example are 
designed using linear quadratic Guassian (LQG) methods, and are similar 
to the designs that appear in [261 . This example goes one step further 


and analyzes a hybrid implementation of these control laws . 

The designs are based on a linearized model of motion in the 

vertical plane. The flight condition is Mach .9 at an altitude of 
25,000 feet. A 4-state, 2-input, 2-output state space model is given 
in Table 5.3. The states and inputs are: 

= forward velocity 
angle-of-attack 
pitch rate 
pitch attitude 


x - 


Sv 

a 

q 

0 


u 



el eva tor/el evon 
canard 


The model leaves out hydraulic actuators with time constants of 70 rad/ 
sec. The uncertainties in the model (data link time delays and bending 
modes) are such that the control loop should be rolled off before 10 
rad/sec [26] . 

5.4.1 L QG Design of SISO Analog Compensator 

In the first design the elevator /elevon input is used to command 
the pitch attitude. The transfer function of the open loop plant is 

-77.81 (s+,0232) (s+1.962) 

2 2 

(s+,2576) (s+5.676) [(s-,6895) + (.2484) ] 


g (s) - 


(5.41) 
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-.02257 


. 00009 


.01233 


-36.62 


-1.90 


11.72 


-18.9 -32.0900 


.9831 -.0073 


-2.6320 


: 


-.9821 


-.4144 


-77.81 


-.76260 


-.00496 


Table 5.3 State space model of highly maneuverable 
aircraft 
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Due to its deliberate relaxed static stability design, the plant is open 
loop unstable, with a pole pair at .6895 ± j.2484. 

The Bode plot of g(s) is shown in Figure 5.19. The bandwidth of 
g(s) is acceptable, hut the magnitude Bode pxOt of g(s) has vo unde- 
sirable features - low gain for w < .1 rad/sec and a resonant peak around 
w = .5 rad/sec. The claim of LQG design methodology is that a k(s) can 
be found that (1) "shapes the loop" to eliminate the undesirable features, 

(2) guarantees closed loop stabxlity, and (3) provides adequate robust- 
ness margins. 

The number of variations of LQG methods is somewhat greater than 
the number of people who know what the abbreviation "I£G" stands for. 

Here we step through a particular LQG design method and provide the 
reader with enough information so that the design can be duplicated. 

The steps are [35] : 

(1) Choose h to shape the magnitude Bode plot of h(sr-A) ^b. 

(2) Choose p to obtain desired bandwidth of k (sl-A) ^b. 

— c 

(3) Solve LQ Riccati equation to obtain k^ . 

(4) Choose process and measurement noise covariance matrices by 
the robustness recovery procedure [27] . 

(5) Solve Kalman filter Riccati equation to obtain k^. 

(6) Use LQ and Kalman filter optimal gains to form analog 
compensator. 

For the linear quadratic (LQ) regulator problem the open loop plant 
is 


x = Ax + bu | 
y = cx ( 


(5.42) 


In this example the plant is controllable and observable. The LQ 
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regulator cost function is 


ORIGINAL PAf?" 1 •• 
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J = 


T T T 

(x h hx + pu u)dt 


( 5 . 43 ) 


The input that minimizes the cost function is 


u = -k x 
-c- 


( 5 . 44 ) 


where 


k 

-c 


1 T 
- b K 
p - - 


( 5 . 45 ) 


and where K > o is the unique solution of the algebraic Riccati equation 


T TIT 

o = AK+KA + hh- - Kbb K 


( 5 . 46 ) 


For the Kalman filter problem the process and measurement noise 
enter the state space system as shown below: 

x = Ax + bu + yC | 

< ( 5 . 47 ) 

y = cx + 0 ; 


The noise sources have zero mean and intensities 


E[C(t)C(s)] = 6 (t-s) 
Ete (t) 0 (s) ] = (t-s) 


( 5 . 48 ) 


The Kalman filter gains are 

k f = ~ c T Z (5.49) 

where E > o is the unique solution of 

o = aZ + Za T + ry T - Zc T cZ (5.50) 

- — — y - - 

After solving the LQ regulator and Kalman filter problems the LQG 
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compensator is formed by 

*k “ { $ “ SE C ' -f 5) *k + *f e ) 

^ (5.51) 

U = -c-k ) 

Tht first step of the design procedure is to choose h to shape the 
magnitude Bode plot of h(sI-A) A b. Let us look ahead a bit and explain 
why we want to do this. The LQ loop transfer function is k c (sI-A) *b, 
and for to such that (k^tjtol-A) *b| >> 1 the following approximation is 
valid: * 

[k^tjtol - A) ^b] * a | h(jtol-A) ^1 for some a (5.52) 


The ultimate objective is to shape the magnitude Bode plot of kg(s) . 

If the robustness recovery procedure is used then for some finite region 
of W: 

kg(jto) 2 k c ( jtol-A) ^b (5.53) 

Therefore, the reason that we want to shape |h(jtol-A) ^bl is that the 
magnitude of the loop transfer function |kg(jto)| will approximate this 
shape . 

The poles of h(sI-A) ^b are the eigenvalues of A. If the zeros are 
at s = -.2 and -2.5 then | h ( jcol— A) *b| has a shape that corrects for 
the undesirable features of the shape of |g(jto)| . An h that places 
zeros at s * -.2 and -2.5 can be computed by the method of Harvey and 
Stein [28] : 


^The constant "a" is used for scaling. It is the shape, not the 
magnitude, of h(sI-A) _1 b that is important. 
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r o —A o 

h = [9.3xlO~ J 1.481x10 -9.061x10 -1.254x10 


The bandwidth of k^CsI-A) depends on the choice of the control 
weight p. A choice of 

p = 1.6 x 10 ^ 


(5.54) 


(5.55) 


results in a bandwidth of 7 rad/sec. 

The Kalman filter is designed by the robustness recovery procedure 
[27] . Choose 


= b 


(5.56) 


and adjust y to select the frequency range over which kg(s) approximates 


k c (sI-A) *b. For this example use 


y = 3.4 x 10 


-7 


(5.57) 


For the choices given above of h, p, y, and y the LQG compensator 
has the transfer function 


. . . -11,750 (s+.195) (s+1.395) (s+6.305) 

k(s) = L 5 2 ~ 

(s+. 02319) (s+1.962) [ (s+261.8) + (261.8) ] 


(5.58) 


The Bode plot of kg(s) is shown in Figure 5.19, The bandwidth is 
7 rad/sec, the phase margin is 64°, and the gain margin is [-14.5dB, 
37.8 dB] . The closed loop system is, of course, stable. 

5.4.2 Conic Sector Analysis of SISO Design 

The analog compensator is converted to a hybrid compensator. This 
is done in a way similar to the earlier examples in this chapter: 


f (s) = 


2500 


(2nd order Butterworth, 


s + 70s + 2500 break frequency = 50 r?d/sec) 


T =« .02 (foJdover frequency = — = 157 rad/sec) 


d (z) = 


-6.233 (Z-. 99611) (z-. 97244) (z- .88119) (z+1) 
(z-. 99954) (z-. 96145) [ (z+. 63763) 2 + (.26229) 2 ] 


(5.59) 


(Tustinized version of k(s), prewarped frequency 
7 rad/sec) 


h(s) = 


1-e 


-sT 


(zero-order-hold) 


The sample rate is 50 samples per second, which is typical for flight 

computers in fighter aircraft. 

Closed Loop Stability by Theorem 3,3 

Conic sector techniques will now be used to determine if the 
hybrid closed loop system is stable. A stable version of the loop 
transfer operator will be placed inside of a cone and then the sufficient 
conditions for closed loop stability given by Theorem 3.3 will be checked. 

The hybrid operator is open loop stable, but the plant is not. The 
loop transfer operator is therefore unstable, and it is not possible to 
place an unstable operator inside of a cone. This problem is circumvented 
by separating the plant into a stable part multiplied by an allpass 
network. The stable part is created by mapping the unstable poles to 
their mirror images about the jto-axis. The allpass network has unit 


magnitude for all to: 
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where 


g (s) 
ap 


[ (s + .6895) 2 + ( . 2484) 2 ] 
t(s - . 6895) 2 + (.2484) 2 ] 


(5.61) 


The hybrid operator and the stable part of the plant are placed 
inside of a cone. The optimal center of the cone is used: 


c (s) = — g hd*f (s) 
T s 


(5.62) 


1/2 


The hold and prefilter are each replaced by (hg s f) > and the radius 
of the cone is computed by 


r(j<0) = | 


20 20 

E E ! (hg f ) (hg f ) 
k=-20 n=-20 y s k n' 

n/k 


1/2 


[d* ( j CO ) 


(5.63) 


The radius is shown in Figure 5.20a. The center and the multiplicative 
radius are in Figure 5.20b. 

The sufficient conditions for closed loop stability given by Theorem 

3.3 are now checked. Condition (i) is that the hybrid operator and the 

stable part of the plant are strictly inside cone (C,R) . This cone has 

just been constructed. Condition (ii) is that the analog system with the 

loop transfer operator c 9 g p^ s ) closed loop stable. This is verified 

by the Nyquist diagram of eg (s) shown in Figure 5.21a. 1 Condition (iii) 

ap 

is that 


| rg (1 + eg ) -1 (jw) | < 1 for all (5.64) 

ap ap 

This condition is verified in Figure 5.21b. Since all three conditions 
of (the appropriately modified) Theorem 3.3 are satisfied, the hybrid 
system is closed loop stable. 

1 The number of open loop unstable poles = 2, which equals the number 

of counter-clockwise encirclements of the -1 point. 
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Figure 5.21b: Stability test of Theorem 3.3. 
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Robust Closed Loop Stability by Theorem 3.5 

The actual plant is g(s) =* g(s) tl + e (s) ] . Assume that the 

m 

multiplicative perturbation is bounded by the constant je^fjco) | < a 
for all to. Theorem 3.5 can be used to find the maximum value of a such 
that g(s) meets the sufficient conditions for closed loop stability. 

The maximum value is Ct = .36, which corresponds to a guaranteed phase 
margin of 21°. 

R obust Closed Loop Stability by Theorem 3.7 

If the discrete closed loop system is stable then a cone can be 

constructed such that -K 1 is outside of the cone. The condition of 

discrete closed loop stability is verified by the discrete Nyquist 

diagram of d*g^*(s) plotted in Figure 5.22a. 

The center of the cone is g(s), and the radius is computed by (3.36) 

of Theorem 3.6. The multiplicative radius |g ^r ( jco) J is plotted in 

Figure 5.22b. By Theorem 3.7 the closed loop system is stable for 

g(s) = g (s) [1 + e (s)] if |e (jco) | < j g 1 r(j 0 )) I . 
in in 

If the multiplicative perturbation is constrained to be constant 
for all co then its magnitude must be less than the minimum value of 
|g r(jco) j, which is Cl = .68. This value of (X corresponds to a guaranteed 
phase margin of 40°, which is less conservative than the 21° obtained by 
Theorem 3.5. The actual phase margin is 46° (obtained from the discrete 
Nyquist diagram of Figure 5.22a). 

This completes the SIE0 version of the high performance aircraft 
example. The analog compensator was designed using LQG methods and then 
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Figure 5.22a: Discrete Nyquist diagram of <3*g d# (jw) 



Figure 5.22b: Multiplicative radius of cone such 

that is outside of the cone. 
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converted to a hybrid con$>ensator . A stable version of the loop transfer 

operator was created by multiplying the unstable plant by an all-pass 

network. Theorems 3.2, 3.3, 3.5, 3.6, and 3.7 were used in various ways 

to analyze closed loop stability and robustness properties. 

5.4.3 LQG Design of Multivariable Analog Compensator 

The high performance aircraft used in this example has multiple 

control surfaces (see [2,26]) that allow independent control of attitude 

and flight path motion. In the pitch axis the basic objective of the 

multivariable design is to achieve independent control of pitch attitude 

and angle-of -attack with approximately equivalent speeds of response. 

This objective is met in the analog design described in [26] . 

The analog design uses LQG techniques. The plant G(s) of Table 5.2 

is augmented to G (s) = — G(s) = C (sI-A ) 1 B by appending integrators 
-a s - -a - -a -a 

to the outputs (pitch attitude ang angle-of-attack) . The Kalman filter 
noise covariance matrices are 


r 

N 


= b (vr 1 ) ) 

I 


where 


V = 


.962 -.274 


-.274 -.962 


Z = diag [1.8 .0077] 


(5.65) 


(5.66) 
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Tbe V and £ matrices are, respectively, the right singular vectors and 

.5 

the singular values of G (ju> ) for w = 10 . They multiply the B 

-a o o -a 

matrix in (5.65) so that the singular values of C (sI-A ) ^ f are 

-a - -a 


approximately equal at (jJ q = 10 


.5 


The LQ regulator weighting matrices Q and R are chosen so that the 
loop transfer function G K(s) approximates the Kalman filter loop trans- 
fer function C (sI-A ) 1 K over a finite frequency range: 

” “3 “I 


T 

c c 


R = 10 


- 6 i i 


(5.67) 


The Kalman filter and Dp regulator problems are solved and used to 
form the LQG compensator K(s) . The analog closed loop system is guaran- 
teed to be closed loop stable. The singular values of G (jw) and G K(jto) 

“3 "“ 3 “ 

are shown in Figure 5.23a, and the "shape" of the loop transfer function 
G K ( joi) is used to analyze various performance measures (contnand response, 

“A* 

bandwidth, disturbance rejection, and so on, see [1]). The singular 
values of I + (G K) ^(jw) are shown in Figure 5.23b and are used to 

*“ *“ 3 ” 

analyze robustness with respect to perturbations inserted at the output 

of the plant [i.e. when G ( joj) is replaced by (I + E )G (ju>)}. In 

-a - -m -a 

Figure 5.23b it is shown that 

a . [I + (G K) _1 ( ju>) ] > a = .54 for all w (5.68) 

mm - -a- — 


itude 
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which indicates that the analog closed loop system will remain stable 
for phase uncertainties up to ±31°, at any frequency, simultaneously in 
any output channel (see Subsection 2.2.7). 

5.4.4 Conic Sector Analysis of Multivariable Design 

The multivariable analog compensator is j K(s) . The integrators 
that were appended to the outputs of the plant are included with K(s) . 
The analog compensator K(s) is now converted to a hybrid compensator. 
The prefilters and holds are chosen to be the same in their respective 
input and output channels: 


F(s) 


T 


D(z) 


H (s) 


2500 


,s +70s+2500> 


I 


.02 seconds 


tustinized version of — K(s) 

s - 

prewarped about 7 rad/sec 



(5.69) 


In this subsection closed loop stability and robustness margins w: 11 
be determined. The conic sector analysis techniques of Theorems 3.6 and 
3.7 will be shown to give conservative robustness margins. 

The discrete version of the multivariable Nyquist criterion is used 
to determine the closed loop stability of the hybrid feedback system. 

The discretized plant is G^*(jw) = (FGH(jU))]*, and the discrete loop 
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transfer function is G d *D*(jw). The Nyquist plot of -1 + det [I + 

G d *D*(ja))] is shown in Figures 5.24a and 5.24b. The first figure (which 

is not to scale) shows the Nyquist plot as the frequency variable to is 
TT IT 

varied from - — to — . The number of counterclockwise encirclements of 

-1 point is 2, which is equal to the number of open loop unstable poles 

1 

of (z) D (z) . Therefore, the hybrid system is closed loop stable . 

An exploded view of the Nyquist plot is shown in Figure 5.24b. 

The nearness to the -1 point is a good robustness margin for SISO sys- 
tems, but not for multivariable systems. Much better is the singular value 

2 

plot of Figure 5.25a, where it is shown that 

0 . [I + (G *D*) _1 (jto) ] > a = .48 for all to (5.70) 

min - -d - — 

Figure 5.25a is used to analyze robustness with respect to multiplicative 
perturbations inserted in the feedback loop just before the digital 
computer, where the physical signal is a discrete sequence. If the 
multiplicative perturbation is a constant diagonal matrix then it can 
be moved in the feedback loop to just after the analog plant. The lower 
bound a = .48 indicates that the hybrid closed loop system will remain 
closed loop stable for phase uncertainties up to ±28°, at any frequency, 
simultaneously in any output channel . 

^The two integrators in D(z) account for the encirclement between to = 0 
and to = 0 + 

2 

Singular values are better than determinants for indicating the nearness 
to singularity of the return difference matrix I + G d *D* ( jto) . 
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Figure 5.24a: Discrete multivariable Nyquist plot. 



Figure 5.24b: Exploded view of 5.24a. 







frequency (rad/MC) 

Figure 5.25b: Multiplicative radius of cone such 

that -K^ is outside of the cone. 
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The conic sector analysis techniques of Theorems 3.6 and 3.7 will 
now be used to analyze robustness. Theorem 3.6 is used to construct a 
cone such that -K 1 is outside of the cone. The center of the cone is 
the nominal plant G, and the radius is the R computed by (3.36). The 
multiplicative radius is r ( jw) ] and is shown in Figure 5.25b. 

According to Theorem 3.7, the hybrid system will remain closed loop 
stable for a multiplicative perturbation of the nominal plant if 


l (to) < 
m 


r (jco) 


a 

max 


[G ( jco) ] 


for all to 


(5.71) 


Unfortunately, this is a conservative result. In Figure 5.25b it is 
shown that 


l (w) < a = .002 for all to (5.72) 

fit 


which corresponds to a guaranteed multi variable phase margin of < 1° . 

More research is needed to remove this conservatism. Some indi- 
cation of why this conservatism occurs is given by the following manipu- 
lation of the corresponding robustness result for multivariable analog 
systems : 


l ( 03 ) 

m 

4=> l (co) 
m 

l (co) 
m 


< o . [I + (GK) 1 ( jco) ] for all to 

min - — 

< 0 _1 [GK (I + GK) _1 (jco) ] for all w 

max — - — 

< O _1 [G { j co) } o [K (I + GK) 1 (jco) ] 

max - max - - — 


(5.73) 

(5.74) 
for all a) 


(5.75) 
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If (5.75) is used instead of (5.73) then the analog robustness test will 
have approximately the same degree of conservatism as the result of 
Theorem 3.7. The problem occurs in (5.75) because the matrix product 
is split up, which is conservative when the minimum and maximum singular 
values of G(jo>) are far apart. * 

This completes the multivariable conic sector example. A multi- 

variable analog compensator was transformed into a hybrid compensator . 

The resulting hybrid feedback system was analyzed with the use of the 

discrete loop transfer function G.D(z) , and then was analyzed with the 

-a- 

use of the conic sectors of Theorems 3.6 and 3.7. The conic sectors 

resulted in a conservative robustness margin, which occured because the 

multiplicative radius is conservative (small) when 0 . (G(jw)) « 

nun - 

a [G ( jco) ] . 

max - J 


Hlote that 1*0 (GG *) < a (G)/a . (G) , which is conservative when 
max — — * max — nun - 


<WS> << W2'- 
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6. MULTIRATE SAMPLING 
6 . 1 Introduction 

The hybrid feedback systems considered up until now have had one 
synchronous sampler. These are the most common types of hybrid feed- 
back systems; however, there exist many applications where the hybrid 
feedback system contains more than one synchronous sampler, each operating 
at a different rate. These are called "multirate sampled-data systems". 

Multirate sampling may occur because the measurements are discrete 
(e.g. radar tracking, sun sensors for dual spin satellites, or image pro- 
cessing algorithms for robot manipulators) and the digital computer used 
to implement the control law is separate and operates at a different sample 
rate. Multirate sampling may also be used if the plant has different time 
scales, in which case multirate sampling can be used to significantly de- 
crease the required computer capacity, i.e. "slow" control loops are 
sampled at a slower rate as compared to "fast" ones. 

A feedback system with two hybrid compensators in the same loop is 
shown in Figure 6.1. This is called a sinjie loop multirate hybrid (SLMRH) 
feedback system. The feedback loop is used for stability augmentation and 
contains a hybrid compensator (modelled by the hybrid operator K^) with a 
sampling period of mT seconds (where m is an integer 1) . The hybrid 
compensator in the forward loop (modelled by has the smaller sampling 
period of T seconds, which allows the output to respond faster to commands 
than to the feedback. The plant has the nominal model g(s) and the multi- 
plicative perturbation e (s) , assumed bounded by |e (jui) I < i ('•>! all id. 

m 1 m 1 m 

All of the components in Figure 6.1 are ft. male input single ou‘put (SISO) . 


1 
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Figure 6.1; Single loop multirate hybrid feedback system 









-233- 


Ccnventional analysis techniques for single rate hybrid feedback 
systems have been extended to the multirate case. A survey with ex- 
tensive references is by Walton [31] . A tutorial treatment of these 
techniques is by Konar and Mahesh [32] . The basic idea is to transform 
the multirate hybrid feedback system into an equivalent single rate 
hybrid feedback system and then apply standard z-transform techniques. 
The only multirate technique used in this chapter is called the "fre- 
quency decomposition" method [33, 34] . Using this method, the SLMRH 
feedback system in Figure 6.1 is broken at some point vhere the physical 
signal is a discrete sequence, and then a z-transform is derived for the 
linear shift invariant loop transfer operator. 

Conic sector analysis techniques also can be extended from single 
rate to multirate hybrid feedback systems. In this chapter they are ex- 
tended to the relatively simple SLMRH feedback systems. Further ex- 
tensions to the much more difficult multiple loop caso will not be 
attempted here.* The material in this chapter is simply a first step 
that demonstrates that further extensions of the conic sector analysis 
technioues are possible to more conplicated multirate systems. 

Conic sectors can be used for multirate hybrid feedback systems 
(either single or multiple loop) just as they can be used for single 
rate hybrid feedback systems and analog feedback systems - to analyze 
closed loop stability, command response, performance, disturbance re- 
jection, and robustness with respect to plant uncertainties. If the 

^Multiple loop multi rate hyorid feedback systems are more common than 
the single loop case because different time scales of the plant 
naturally break down into multiple loops. 
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multirate hybrid operator is placed inside of a conic sector, then the 
center can be used as z. rigorous continuous time linear time invariant 
( CTLTI ) approximation . 

Conic sector analysis techniques are potentially more important for 
multirat^ than for single rate hybrid feedback systems. The feedback 
loop is Broken where the physical signal is analog and complexities due 
to different sampling rates are subsumed by the use of conic sectors. 
Techniques to analyze robustness do not currently exist, whereas conic 
sector techniques are inherently robustness techniques (as discussed 
in Subsection 2.4.5). The z-transform techniques used to convert multi- 
rate to single rate systems have the problem of increased dimensionality 
(proportional to the integer multiple of the sample rates) . The conic 
sector techniques do not suffer from this increased dimensionality. 

In Section 6.2 the SLMRH feedback system is described and frequency 
domain input-output transformations for the multirate hybrid operators 
are presented. Theorems 3.1 to 3.7 for single rate hybrid feedback 
systems are extended to SLMRH feedback systems. These extensions are 
presented in the multiple-part Theorems 6.1 and 6.2. The new conic 
sector analysis techniques can be used to analyze stability and robust- 
ness, but the techniques for analyzing command response (as in Subsection 
2.4.6 and Section 3.5) have not yet been extended. In Section 6.3 a 
multirate version of the lead-lag compensator example of Section 5.2 
is presented. 

6.2 Conic Sectors for Multirate Hybrid Operators 

Conic sector analysis techniques are used by first dividing a 
feedback loop into two subsystems, and second by constructing a cone 
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such that it contains one of the subsystems and the other subsystem is 
outside of it. There are numerous ways to divide up the SLMRH feedback 
system of Figure 6.X. In this chapter the feedback loop is broken at 
points 1 and 2 (at the input and output of the plant) . The two sub- 
systems are (1) the multirate hybrid operator and (2) the actual 

plant G. A conic sector is derived in Subsection 6.2.1 that contains 
and in Subsection 6.2.2 that is outside of. 

6.2.1 Existence of a Cone that Contains K^K 2 

The multirate hybrid operator shown i n Figure 6.2. The 

analysis that follows assumes that the slower sampler is first*. The 
existence of a conic sector that contains K^K 2 follows from the de- 
scription of the input-output transformation from y to u. The assump- 
tion is made that the two samplers are in synchronism, i.e. every mth 

2 

sample the two samplers coincide. . 

The transformation from y to y 2 is described by* 

y 2 (ju>) = h 2 (ju))d 2 (e jamT ) ^ E f 2 ( job- j ~ k)y(jw-j ^ k) (6.1) 

k 


A similar analysis can be performed if the faster sampler is first. 

2 

If the samplers are not in synchronism then an extra delay term must 
be inserted between h„(s) and f. (s) . This does not invalidate what 

1 ~ T d s 

follows. Simply replace h 2 (s) by h^sje 
3 

To avoid confusion oetween sanple times of T and mT seconds the star 
notation (see Subsection 2.3.1), the symbol "w " and the subscripts 

indicating frequency shifts are not used. The notation d 2 (e^ UlnT ) 

indicates d^z) evaluated at 2 = 
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The multi rate hybrid operator K K 

1 2 


Figure 6.2: 
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Similarly, the transformation froir. to u is described by 

u(jo)) = h 1 (ja))d 1 (e : ’ L0T ) ^2 f x ( jto-j ^ n) y 2 ( ju>- j ^ n) (6.2) 

n 

Substitute (6.1) into (6.2) and change the variable of summation to 
obtain 


u(jtu) = £ 1 (jU))d 2 (e jU1Tir ) ~ 2 f 2 ( jw- j k)y(jw-j ^ k) (6.3) 

k n ' 


where 


C x (» = h x ( jw)d 1 (e jU,T ) | 2 f x ( jta-j ^rn)h 2 (ju>-j ^ n) (6.4) 

n 

The multirate hybrid operator described by (6.3) has the same structure 
as the single rate hybrid operator described by (2.51), except that the 
hold h 2 (jw) is replaced by f^(jw). it is precisely this similarity in 
structure that allows the conic sector results to be extended from the 
single to the multirate case.'*' 

The following five part theorem is an extension of Theorems 3.1 
to 3.4. 

Theorem 6,1 Let K^K 2 be the multirate hybrid operator of Figure 
6.2. Assume that is L_ -stable. Let C be any LTX L „ -stable 

operator, and let R and R 1 be LTI l 2e ~stable operators. 

^The same structure ii found if the higher rate sampler is first. The 
conjecture is made that the same structure will be found if the sample 
rates form a ratio that can be expressed as a rational number. 
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(a) ^^^2 strictl y inside cone (C,R) if 


! r ( joj) 1 _> ~~i/2 lr 4^ + 


1/2 


(1-e) 


( 6 . 5 ) 


for all a) and seme £ > 0 


where 


r 4 (w) 


(mT) 


<y e j “' T ) 


z z 

k n?0e 


r / ■ . 2 tt , . . . . . 2tt .2 


( 6 . 6 ) 


r 3 <u) 


1 r , . . 271 .... . iwmT. ... . 277 . , , . . 2TT , . 

S 5f k)d 2 (e k >-'(3“-3 Sf k) 


(b) The optimal center 


( 6 . 7 ) 


C(S) = ^ C l (s)d 2 (eSmT)f 2 (s) 


( 6 . 8 ) 


minimizes the lower bound for |r ( joj) | . 

(c) The gain of the multirate hybrid operator is 


II ¥ 2 ", 


2 0<o)<— 

mT 


jl 1 , l 

V^Sk) 

A 1 ? 

d 2 (e 3 “" r > 

( l (mT) k 


J 



Z l f 2 (3 “- j S n) 

.n- 


,u /2 


(6.9) 
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(d) The SLMRH feedback system of Figure 6.1 is closed loop stable 
for G = G (i.e. for the nominal plant with no perturbation) if 
a C and R exist such that 


(i) 2 is strictly inside cone (C,R) 

(ii) G(I+CG xI is L^-stable (i.e. the nominal analog system 
is closed loop stable) 

(iii) |rg(l+cg) 1 (ju))l < 1 for all u (6.10) 

(e) The SLMRH feedback system is closed loop stable for all pos- 
sible G's if in addition to the three conditions of part (d) 
the following condition is true: 


m 


< U+cg(ju)l - 1 rg ( j(b) | 
cg(jw)| + |rg(ja))| 


for all to. 


( 6 . 11 ) 


The proof of Theorem 6.1 is similar to the proofs of Theorems 3.1 
to 3.4. The following steps prove part (a) of Theorem 6.1 by showing 
that the conic sector inequality (2.81) is satisfied for all input- 
output pairs defined by K^K 2 . 

(1) Define truncated function to convert from truncated to 
function norm. 

(2) Use Parseval's theorem to convert from time to frequency 
domain . 

(3) Use frequency domain inequality of Lemma 3. A. 

(4) Use Parseval's theorem again to convert from frequency to 
time domain. 

(5) Convert from to truncated function norm. 
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Part (b) is proved by showing that the optimal center (6.8) makes r 3 (u)=0 
and therefore minimizes r 4 (w) + r 3 (w) of (6.5). The proof of part (c) 
uses the center c(s)=0 to find an upperbound for the gain of and 

then constructs an input signal that achieves the upperbound. Part (d) 
is an application of Lemma 2.5 to show closed loop stability. Conditions 
(ii) and (iii) of part (d) guarantee that -G 1 is outside cone (C,R) . 

If inequality (6.11) of part (e) is satisfied then all possible -G^s 
are outside cone (C,R). 

If Theorem 6.1 is used to analyze a SLMRH feedback system then the 
nominal feedback system is analog and has the loop transfer function 
c(s)g(s). Stability and robustness properties of the SLMRH feedback 
system depend on the stability of the nominal feedback system. The 
radius of the conic-sector is treated as an additive perturbation of 
c(s) . 

To reiterate. Theorem 6.1 has the following five parts: 

(a) sufficient conditions for the existence of a cone that con- 
tains K . K . 

1 2 

(b) aui optimal center for the cone 

(c) the gain of K^K ^ 

(d) sufficient conditions for closed loop stability 

(e) sufficient conditions for robust closed loop stability 

The analysis techniques laid out in Theorem 6.1 represent a new approach 
to the analysis of multirate hybrid feedback systems . 

6.2.2 Existence of a Cone such that -(K^K^) 1 the Cone 

The alternative to constructing a cone that contains K^K 2 is to 
construct a cone such that - (K^K^) is outside of the cone. The best 
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choice for the center of this cone is G, the nominal plant. It follows 
from Lemma 2.2 that is outside of cone (G,R) if an R can be 

found such that j|R K^K^(1+G The only restriction is 

^ j 

that the multirate hybrid operator K K (I+G K K ) must be L_ -stable. 

11 12 2e 

The multi rate hybrid feedback system that has the closed loop 
operator R K^K^O+G K is shown in Figure 6.3a. By breaking the 

loop at point 1 and manipulating the block diagram the equivalent feed- 
back system shown in Figure 6.3b is obtained. The feedback loop is 
digital and operates at the slower sample period of mT seconds. The 
multirate hybrid operator K^K^il+G K^K J 1 is I- 2 e -sta bl e if this digital 
feedback loop is stable. 

It is possible to obtain the z-transform of the loop transfer func- 
tion of Figure 6.3b by the frequency deconqposition method, see [31] to 
[34] . Define the z-transform of the discretized plant as g d (z) • The 
loop transfer function is d^(z)g^{z) , and stability is checked by the 
discrete Nyquist criterion, which uses a plot of d 2 (e^ tJmT )g { j(e^ amT ) 

TT 

from u) = 0 to — rad/sec. 
mT 

It is not necessary to go to the considerable trouble of finding 
the z-transform g d (z). it is computationally easier to obtain g d (z), 
evaluated at z=e‘* almT , by summing truncated versions of the following 


infinite series: 
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, jumT. 

g d (e ) 


JS E 5 k > f i<5“-3 5F SF *> 

k 


( 6 . 12 ) 


where 


C 2 <j“) = 


f I f 2 (ju-j 
n 


2tt , . . . 2tt 

— n)g(3W-D — n)^ 


. . . 2tt 

(DOJ-d f n) 




(6.131 


Before constructing a cone that ~(K^K^)^ is outside of, an input- 
output description of the multirate operator K^K ^(J+G K^K ^ ) must be 
found. Using Figure 6.36b and equation (6.4), (6.12), and (6.13) we 
obtain 


u(jai) = ? 1 (j a -° j c £(e^ UJmT ) 


5F l E 2»“' 

k 


■ 2lT 

3 inT k y 2 


(jw-j ~ k) (6.14) 


mT 


where 


d ci (e3la ’ T) = 


-1 




(6.15) 


Equation (6.14) has the same form as the single rate hybrid com- 
pensator of (2.51) , which makes possible the following extension of 
Theorems 3 . 6 and 3.7: 

Theorem 6.2 Let K,K,, be the multirate hybrid operator of Figure 

-L 

6.2, let G be any LTI operator such that K^K^I+G K^K^ 1 is L^-stable, 

and let R and R 1 be LTI -stable operators. 

2e 

(a) - (K^K ^) 1 is outside cone (G,R) if 
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lr(ju)) I < for all oj 

where 


(6.16) 


r 9 (w) 


~^~2 [ jj-j k) i 2 l • i d C £( ejwmT > r 

(mT) Ik 


, , , . 2tt , 1 2 

£ 2 (3 “" :i mT n I 


(6.17) 


(b) The multirate hybrid feedback system of Figure 5.1 is closed 
loop stable for all possible G's if an r(jw) exists such that 

(i) K K (1+G K K ) J is L -stable 
12 12 2e 

(ii) is outside cone (G,R) 

(iii) £ m (td) < |g ^r(jco)| for all u). • (6.18) 

If is placed outside of a cone, then the nominal feedback 

system is the SLMRH feedback system with the nominal plant G. The 
stability of the nominal feedback system must be checked, which is done 
by constructing an equivalent single rate hybrid feedback system. It 
is not necessary to explicitly find a z-transform description of the 
single rate system. Instead, stability can be checked with the use of 
equations (6.12) and (6.13). 

The center of the cone that outside of i s the nominal 

plant G. Part (a) of Theorem 6.2 gives a radius for this cone. Part 
(b) is a robustness result which gives sufficient conditions for all 
possible perturbations of the niminal plant to preserve closed loop 


stability. 
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6 . 3 Multirate E xample 

This example is a continuation of the lead-lag compensator example of 
Section 5.2. A second hybrid compensator, operating with twice the sample 
period, is added in the feedback loop. Theorem 6.1 is used to construct 
a cone that contains the multirate hybrid operator and to determine 

closed loop stability. Theorem 6.2 is then used to construct a cone such 
that — * s outs ^^ e cone an< * to determine a robustness margin. 

This multirate example will appear to be similar to the single rate 
example of Section 5.2.^ Though somewhat repetitious, this highlights the 
fact that conic sector analysis techniques are just as useful for multirate 
as for single rate hybrid feedback systems. 

The SLMRH Feedback System 

A block diagram of the SLMRH feedback system is shown in Figure 6.1. 
The transfer function of the nominal plant G is 


g(s) 


150 

(s+1) (s+3) 


(6,19) 


and the hybrid compensator (in the forward loop) is a hybrid implemen- 
tation of the analog compensator 


k(s) 


(s+3) 2 

(s+.4) (s+22.5) 


( 6 . 20 ) 


The prefilter, sample period (T = .031416 seconds) , digital computer 
z-transform, and the hold are given in equation (5.3). 

The components of the hybrid compensator K (in the feedback loop) are 


Specifically, Subsections 5.2.5 and 


5.2.7 
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shown below; 


f 2 (s) 


2T * 
d 2 (z) 

h 2 (3) 


625 


s +35s+625 


(2 nd order But.terworth with break at 
25 rad/sec) 


.062832 seconds (m * 2, ^ * 50 rad/sec) 
1 2z 

■ 2 (digital lead compensator) 

. -2sT 
1— e 

* — (zero-order-hold) 



( 6 . 21 ) 


In the single rate example a continuous estimate of the output signal 
was available. Here the estimate is undated every 2T seconds and held 
constant between updates. The prefilter f 2 (s) smooths the signal (prevents 
aliasing!. The digital computer serves the purpose of adding phase lead 
around the crossover frequency to partially compensate for the phase lag 
due to the prefilter and hold.^ 

A Cone that Contains K , K 

1 2 


The feedback loop is broken before and after the analog plant, where 
the physical signals are analog, and the multirate hybrid operator K^K 
transforms the output signal y to the input signal u. 

Part (a) of Theorem 6.1 is used to construct a cone that contains 
* 1 *, . The optimal center is chosen: 


c(s) =-^-^(8) d 2 (e SmT ) f 2 (s) 


( 6 , 22 ) 


At to = 7 rad/gec the prefilter f (s) decreases the phase by 23°, tne held 
h (s) decreases the phase by 13® f and the computer d^(z) increases the 
phase by 4*, rhe phase margin of the single rate system of Section 5.2 
is 46*, here we expect it to be about 46-23-13+4 = 14* (it turns out to 


be 11*). 



-247- 


where (s) is defined by (6.1). In Figure 6.4 the magnitude and phase 
Bode plots of cg(joo) and kq(jw) are compared. The magnitude of cg(jo)) 
drops off sharply at about 20 rad/sec and has the expected extra phase 
lag at the crossover frequency of 7 rad/sec (phase margin = 14°). 

The radius is computed by (6.5) 1 and is shown in Figure 6,5a, The 
multiplicative radius c ^r(jw) is shown in Figure 6.5b and it> compared to 
cg(joi). Over the bandwidth of the nominal analog system (ou < 7 rad/sec) 
the multiplicative radius is <.12. 

Closed loop stability is determined by the three condi tiens of part 
(d) of Theorem 6,1. Condition (i) is that is strictly inside cone 

(CfR), which is true by the way that the cone was constructed. Condition 
lit) is that the nominal analog system is closed loop stable, which is 
verified l, y the Nyquist plot of cg(joj) in Figure 6,6a. Condition (iii) is 
that 

|rg(l + eg) * ( jw) | < 1 for all CjJ (6.23) 

which is verified in Figure 6.6b. All three conditions are true, hence the 

SIMRH feedback system is closed loop stable. 

A Cone such that -(K K~,> is Outside of the Cone 
I «- 

The nominal feedback system differs depending on how the cone is applied,, 

When is inside of a cone then the nominal fecucack system is analog. 

In this part of the example is outside of a cone, and the 

^Because the optimal center was chosen r (joj) = o. The infinite sum for 
t^(j^) was truncated at ±30 terms, and tne double infinite sum for r^(jw) 
was truncated at ±20 terms. 


phaac angla iwagnituda 




magnitude magnitude 
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Figure 6.5a: 


Radius of cone that 


contains 



Figure 6.5b: Multiplicative radius 
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nominal feedback system is the SI^MRH feedback system with the nominal 
plant G. Its stability is verified by the discrete Nyquist plot of d^g^Cz) 
shown in Figure 6.7a. From this plot we obtain the discrete phase margin 

of 11°. 

Part (a) of Theorem 6.2 is used to construct the cone. The center 
is the nominal plant G, the radius is computed by (6.16), and the multipli- 
cative radius q ^"r(jto) is shown in Figure 6.7b. According to part (b) of 

Theorem 6.2, the SIflRH feedback system remains stable for any g = g(l + e ) 

m 

if 


l (o» < 

m 



r(jco) | 


for all co 


and if L (co) is constrained to be constant then 
m 


(6.24) 


L (co) < a = .19 for all co 
m 


(6.25) 


This constant corresponds to a phase uncertainty up to 11° at any frequency 
(see Subsection 2.2.7). 

This completes the multirate example. The conic sector analysis techni- 
ques of Theorem 6.1 and 6.2 have been used to determine closed loop stability 
and a robustness margin. These conic sector techniques are applied just as 

were the corresponding techniques developed in Chapter 3 for single rate 

2 

hybrid feedback systems. 


X 0f course we already knew the SLMRH feedback system is closed loop stable 
form part (d) of Theorem 6.1. 


2 

The second hybrid compensator has not affected closed loop stability, but has 
signif icantly decreased the robustness margins. A step response will exhibit 
large overshoot. This design is probably not acceptable, but can probably 
be corrected by adjusting d^(z) and/or d^(z) to obtain more phase lead 
around the crossover frequency . 
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7. SUMMARY AND SUGGESTIONS FOR FUTURE RESEARCH 
7.1 Summary 

The major result of this thesis (Theorem 3.2) is the determination 
of a new conic sector which contains a stable hybrid operator. Every- 
thing else in this thesis leads up to , proves , modifies , constructs , 
demonstrates , and extends this basic result. 

The mathematical preliminaries of Chapter 2 lead up to Theorem 3.2. 
Sections 2.2 and 2.3 contain a review of multivariable (singular value) 
analysis techniques for analog (CTLTI) and digital (DTLSI) feedback 
systems. A hybrid system with a single synchronous sampler is defined 
and its properties are discussed. In Section 2.4 a precise mathematical 
framework is set up for a conic sector analysis of nonlinear time varying 
feedback systems. In this general framework, conic sectors are used to 
analyze closed loop stability, robustness properties, and steady state 
response to commands. The conic sector analysis techniques are not useful 
unless a specific conic sector can be found for the feedback system of 
interest. In Section 2.5 specific conic sectors are determined (due to 
Safonov [7,9]) for analog (CTLTI) feedback systems and then used to 
analyze closed loop stability and robustness properties. 

Theorem 3.2 gives sufficient conditions for the existence of a new 
conic sector v/hich is specifically for use in analyzing hybrid feedback 
systems. Part of Chapter 3 is used to prove this basic result. The 
most important part of this proof is Lemma 3. A, which is a frequency 
domain inequality that applies to all possible input-output pairs of 
signals of the hybrid operator. 

The remainder of Chapter 3 (there are a total of 8 theorems) is 
used to modify Theorem 3.2. The entire loop transfer operator (not just 
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the hybrid operator) can be placed inside of a conic sector, but what- 
ever is inside of a conic sector must be open loop stable. Theorem 3.2 
is modified to remove this open loop stability restriction, so that the 
inverse of any hybrid operator (stable or unstable) can be placed outside 
of a conic sector (Theorem 3.6), The conic sectors of Theorems 3.2 and 
3.6 are used for analysis techniques for hybrid feedback systems (just 
as for the general feedback system of Section 2.4) to analyze closed 
loop stability, robustness with respect to modelling uncertainties, 
and steady state response to command inputs. 

The algorithms and numerical results presented in Chapter 4 are 
used to construct the new conic sectors. The center is arbitrary, 
though some are better than others, and there is not any difficulty in 
computing it. The radius, however, can be difficult to compute, and 
various ways to do so are presented and discussed. 

The examples in Chapter 5 demonstrate the usefulness of the new 
conic sector analysis techniques. In addition to the analysis techniques 
mentioned above, the conic sectors are used to select the sample rate 
and compare discretization techniques. In the one multivariable example, 
the robustness results are shown to be conservative. 

The material in Chapter 6 extends the conic sector results of 
Chapter 3 to single loop multirate hybrid (SLMRH) feedback systems. 

This is done by combining the two hybrid operators with different sample 
rates into a multirate hybrid operator which has an input-output trans- 
formation with a structure similar to a single rate hybrid operator with 
the slower of the two sample rates. Once this similai structure is 
established then Theorems 3.1 to 3.7 are extended and then demonstrated 


by an example. 
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The mathematical development leading up to and proving the major 
result (Theorem 3.2) is highlighted as being of interest to control 
thereoticians . The rigorous justification for treating hybrid feedback, 
systems as analog feedback systems, and the specific analysis techniques 
based on Theorem 3.2 and its modifications are highlighted as being of 


interest to control practitioners. 
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*.2 Suggestions for Further Research 

The idea of using conic sectors to analyze hybrid feedback systems 
originated with Gunter Stein and is now about 3 years old. This thesis 
has developed this idea to the point where it can now be said that conic 
sector analysis techniques are constructive and useful for hybrid control 
systems. Time and further developments are needed, however, if these 
conic sector analysis techniques are to become generally accepted. 

Removal of Open Loop Stability Restriction The major restriction 
to the use of conic sectors that contain hybrid operators is that the 
hybrid operator must be open loop stable. This restriction prevents the 
use of some of the conic sector analysis techniques for important control 
systems such as those with integral control action. While it is possible 
to eliminate this open loop stability restriction by placing the inverse 
of the hybrid operator outside of a conic sector,'*' it is still desirable 
to be able to place both open loop stable and unstable hybrid operators 
inside of conic sectors. The advantages of having the hybrid operator 
inside of a conic sector are (1) the center of the conic sector can be 
used as an analog (CTLTI) approximation of the hybrid operator, and 
(2) the nominal feedback system is analog (and therefore one of the 
sufficient conditions for hybrid closed loop stability is that the 
nominal analog system is closed loop stable) . 

It is not clear how (or even if) this open loop stability restriction 


This alternative is significant and leads to the important robustness 
analysis technique of Theorem 3.7. The assumption of open loop stability 
of the hybrid operator is replaced by an assumption that the nominal 
hybrid feedback system (with the nominal plant) is closed loop stable. 
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can be removed, This restriction can be traced back to the Small Gain 
Theorem (see Subsection 2.4.4), which gives sufficient conditions for 
closed loop stability, but requires that the relations that define the 
feedback system be open loop stable. 

One suggestion for removing the open loop stability restriction is 
to define the conic sector using the gap metric [36, Section 4.2, and 37]. 
The basic idea behind the use of the gap metric is that there is some 
exact information which must be known about an unstable system before it 
can be stabilized. For analog (CTLTI) systems this exact information is 
the number of open loop unstable s^plane poles. For hybrid systems this 
exact information may be (it still remains to be shown) the number of 
unstable s-plane poles (of the analog components) plus the number of 
unstable z-plane poles (of the digital components) . 

Generalization of Conic Sector Techniques to Sector Techniques . A 
conic sector is a special case of a sector (see Safonov [7]). it should 
be possible to develop analysis techniques for hybrid systems based on 
the use of sectors. If an operator (of any type) is placed inside of a 
conic sector then the radius of the conic sector is analogous to an 
additive perturbation. Sectors, on the other hand, can be set up so 
that parts of the sector are analogous to addition, multiplication, 
subtraction, and division perturbations (see Lehtomaki [38, p. 86}). 
Sectors can be used to analyze combination operators. If G^ and 

sure inside of sectors then it is possible to find sectors that contain* 

G +G G -G , G G , and G G 1 . These types of operations are needed to 
12 12 12 12 

*Conic sectors can be used for addition, subtraction, and multiplication, 
but the results are conservative. 
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analyze complicated feedback systems with multiple loops (such as multi- 
rate multiloop hybrid feedback systems) . References for combining sectors 
are by Safonov [7, 10]. 

Less Conservative Multivariable Robustness Results . The conic 
sector analysis techniques developed in Chapter 3 are valid for both 
SISO and multivariable hybrid feedback systems, but have been shown by 
example to be conservative for multivariable systems.* This problem 
with conservatism must be corrected. It may be just a matter of breaking 
the feedback loop at a different point in order to analyze robustness, 
but the details remain to be worked out. 

Even for SISO systems the conic sector analysis techniques are 
sensitive to where the feedback loop is broken (see Section 5.2). This 
property is counter-intuitive (at least it is counter to our intuition!) 
and needs to be better understood. 

Synthesis Techniques The emphasis of this thesis has been on the 
development of analysis techniques. What the community of control 
system designers really need, however, are synthesis techniques. 

Therefore, the development of conic sector synthesis techniques is an ill- 
defined but necessary extension of this thesis. 

In some sense the conic sector analysis techniques developed in 
this thesis can be considered synthesis techniques. They can be used 
to compare alternative designs for hybrid compenstators , and they can be 
used as part of an iterative design process (at each step of the 
iteration "engineering judgement" must be used to make changes in the 

1 

In particular, the multivariable robustness technique of Theorem 3.7 
is conservative if the singular values of the nominal plant are far 
apart (see Subsection 5.4<-4) . 
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design) . In the examples it was shown how the gain of the hybrid 
operator can be used to help select the prefilter, and it was shown how 
the magnitude of the multiplicative radius as a function of frequency 
can be used to select the sample rate. 

Ultimately, what conic sectors have to offer is a rigorous justifi- 
cation for keeping the design of hybrid compensator in the analog domain. 
The point of view is that the hybrid compensator is supposed to mimic an 
analog compensator, and the extent to which it does not is a source of 
error (i.e. a perturbation of the nominal analog comepnsator) which is 
treated just as a perturbation of the nominal analog plant. Specific 
algorithms for transforming an analog compensator into a prefilter, 
sampler, digital computer, and hold should have the goal of keeping 
the perturbation of the nominal analog compensator small relative to 
the perturbations of the nominal analog plant. 

Multi rate Sampling Issues Chapter 6 is a preliminary extension 
of conic sector analysis techniques to multirate hybrid feedback systems. 
More work needs to be done in this area, specifically (1) extend the 
results of Chapter 6 to multiloop and multirate hybrid feedback systems 
and (2) relax the assumption that the sample rates form an integer ratio. 

In Chapter 6 two single rate hybrid operators K ^ and were 
combined to form a multirate hybrid operator K ^ K ^ , and then a conic 
sector was constructed that contains ' * different approach that 

should be tried is to place the single rate h/brid operators and K 
into their own respective conic sectors, and then to combine the centers 
and radii to form a conic sector of the multirate hybrid operator. To 
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do so in a nonconservative way may require the use of sectors as opposed 
to conic sectors. 

Asynchronous Sampling Issues One of the assumptioi s made throughout 
this thesis is that the single rate hybrid compensator contains a synchronous 
sampler. This assumption is not valid if the sample periods are 

variable, random, or dependent on the dynamics of the plant. These 
types of asynchronous samplers are rarely used, so there is little need 
to develop conic sector analysis techniques to handle them. It is just 
as well, because we do nor have any suggestions for how to do so. The 
input-output transformation must be defined, and then the conic sector 
inequality (2.81) must be shown to be valid for all possible input- 
output pairs of signals. Unfortunately, there are no systematic ways 
to determine conic sectors. 

More common then asynchronous sampling is synchronous that skewed 

sampling, which occurs when the computational delay is a fraction of a 

sample period. This thesis has not explicitly addressed the issue of 

skewed sampling, but the conic sectors of Chapters 3 and 6 can be 

modified to account for skew by changing the hold transfer function from 
- T d s 

h(s) to h(s)e . This is easier, for instance, then using modified 
z-transforms for digital control systems. 

Asynchronous sampling may occur in multirate hybrid feedback systems, 
especially if the clocks that control the sample intervals are not 
synchronized. Typical behavior is that the sample times slowly drift 
with respect to each other. The z-transform techniques developed for 
multirate systems are difficult to apply to this problem. Conic sectors 
■or sectors) offer the following approach: each hybrid compensator is 
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placed in its own conic sector, and then the conic sectors are combined 
to form a conic sector of the composite operator. By treating each 
hybrid compensator separately it should not matter whether or not they 
are synchronized. 

Finite Wordlenqth Issues If the effects of finite word length of 
the digital computer are included in the model of the hybrid compensator 
then the model is a nonlinear time varying operator. In this thesis 
the effects of finite word length are assumed to be neglicable, thereby 
allowing the model to be a linear time varying operator. With the 36 
and 32 bit microprocessors currently available it is indeed very safe 
to assume that the effects of finite word length are neglicable. 
Nevertheless, this is a problem that has received a lot of attention 
in the digital control literature (see (4, Chapter 7] and [39]). The 
practical effects of finite word length are (1) truncation errors in 
computing, (2) truncation errors in stored parameter values, and (3) 
limit cycles in the feedback loop. 

A motivation for further reserach is to determine if conic sectors 
are useful for the analysis of hybrid feedback systems when the finite 
word length of the digital computer is taken into account. As mentioned 
for asynchronous sampling there are not any systematic ways to determine 
conic sectors, so we cannot offer any specific guidance. If such a conic 
sector is found, however, then its center can be used as an analog 
(CTLTI) approximation of the hybrid compensator; and then the errors 
due to the finite word length can be bounded by the radius and treated 
as a perturbation of the nominal compensator. 

Conic sector analysis techniques were developed for nonlinear 
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systems, so the fact that the effects of finite word length are nonlinear 
does not exclude the use of conic sectors. Difficulties will probably be 
encountered when the conic sector inequality ( 2 . 81 ) is checked for all 
possible input<-output pairs. Also, difficulties may be encountered because 
the definition of stability must be changed so that limit cycles are not 
considered unstable. 

Even it a conic is found it may not be useful due to conservativeness 
(as was the case for the first two conic sectors developed for linear time 
varying hybrid operators) . The particular problem with the assumption of 
finite word length is that input signals of very small amplitude may fall 
below the first truncation level and produce zero output. Linear operators 
used for cone centers and radii cannot distinguish between small and large 
amplitudes, and therefore the radius will treat large amplitude signals 
at any frequency as though their output was zero. 
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